Started by Hasdeo.2013.05.31 This is Moya's logbook for Nano Japan project from 2012.6.2-2012.7.31. #contents * Daily schedule [#daily] - 09:00-10:00 Discussion (1) - 10:00-12:00 Solve the problem and write the progress on pukiwiki - 12:00-13:30 Lunch - 13:30-14:30 Discussion (2) - 14:30-16:30 Solve the problem and write the progress on pukiwiki - 16:30-17:30 e-mail report and problems * Goal of the project. [#goal] - To describe the thermal conductivity of Nanotubes and graphene. The output will be an animation of atomic vibration and its propagation in Mathematica - Keywords: Thermal conductivity, carbon nanotubes, Differential equation, Fourier Transform, Solid-state physics, Mathematica * Schedule for discussion [#discussion] |Time and Day|Mon|Tue|Wed|Thu|Fri| |09:00-10:00|Hasdeo|Tatsumi|Mizuno|Thomas|Nugraha| |13:30-14:30|Thomas|Mizuno|Tatsumi|[group meeting]|Hasdeo| * Division of tutors [#v0d44b28] |Name|Lesson to teach|Misc| |Hasdeo|Phonon in Nanotubes| How to make a presentation, and Japanese class| |Tatsumi|Nanotubes| Printer and Copy machine| |Thomas|Physical Mathematics| Library| |Mizuno|Physical Computation| Computer in Lab| |Nugraha|Phonon in Nanotubes| | - Tutorial takes place in Coffee room, since visitor room will be occupied with many people * Questions and Answers [#QA] - This section is for posting questions from Moya-san and answers from other group members. - Please list here with some simple reasons or details. - For every problem, give a tag double asterisks (**) in the code so that it will appear in the table of contents. - For the answer, give a tag triple asterisks (***) in the code below the problem in order to make a proper alignment. - List from new to old. * Report: [#d6196c2a] ** June 3: [#k93c5c39] Today, Hasdeo-san picked me up from Urban Castle Kawauchi (UCK) with Thomas-san. I learned how to get to the lab by campus bus and have been introduced to the lab, Saito-sensei and the secretaries. Hadeo-san taught me about coupled oscillators, where I learned how to solve for the motion of two masses connected by springs. The goal is to model heat in CNT's as atoms vibrating with springs. Thomas-san also lectured me and taught me some math. He taught me how to solve the 2-dimensional Laplacian using separation of variables. He also provided some insight on the applications solving the laplacian. So far, everyone seems very nice and I am very excited to start working on thermal conductivity, which was a surprise. ** June 4: [#p63feebc] Today, I walked from Urban Castle Kawauchi to the lab and managed to find my way, but I might have strayed a bit! I came early to skype my family back home, but began my lessons with Tatsumi-san at 9:00. He taught me about the graphene and carbon nanotube lattice and how vectors compose the unit lattice to form the Chiral Vector and Translational Vector. He also showed me formulas which describe the nanotube lattice. For homework, I was assigned to build a carbon nanotube out of paper and I have already drawn the Chiral Vector, but have not finished the job as I need to borrow scissors and tape. I also spoke with Saito-sensei and discussed the oscillator between two walls. I have a much better understanding of this situation and Saito-sensei assigned me to work on a similar model including a drag force. I believe I understand this problem, but am not to sure about the graphs which show resonance. After lunch, Mizuno-san introduced me to Mathematica. I do not understand it that well, but for homework I am to research it online and learn different functions. I also spent some time playing with the functions and getting the hang of working with differential equations and plotting. I am also supposed to work on an animation, but this seems rather difficult at the time. Saito-sense also showed me how to calculate different aspects of the CNT today. I learned a lot from our meeting and am very excited to work on the problems he sets up for me. I am to work on a time dependent oscillating function for homework tonight and hopefully I can get this one right and continue to expand. ** June 5: [#gf3952b1] I walked to the lab again today and met with Mizuno-san again to discuss Mathematica. He taught me about Rules and Patterns in programming and showed me several functions in mathematica. For homework I am to study graphics and attempt to make a figure of a particle between springs. I am also to study animation and manipulation, so I can modify this graphic. These are not very easy tasks as I do not know much about programming and I plan on looking for material to help me on mathematica. I have not begun the homework yet as I met with Saito-sensei to discuss more on a particle between two springs with an applied force. I believe I am starting to get the hang of solving these differential equations as I am beginning to recall what I learned last semester in my Differential Equations class. Homework from Saito-sensei involved me recalculating the solutions to the equations and plotting them on mathematica. I think I have successfully plottted them, but plan to ask Saito-sensei if it is a correct graph. I will try to solve the differential equation of the particle with friction and later try to derive the calculations for the properties of the CNT that Saito-sensei showed me yesterday. After lunch, I worked on the plot of the solution on mathematica and went online to learn different functions of mathematica. I learned some shortcuts in typing and how to edit graphics, but I still have no idea how to make good graphics or animations.I met with Tatsumi-san and he explained some basic CNT properties like chiralities and metallic/semiconductor conditions. After this we discussed relationships between energy and temperature and energy and wavelength. We then went on to attempt to solve the Heat Equation. I understood the steps to solve the heat equation up until we had to solve for the constants. Solving for the constants required some ugly integrals with summations and trigonometry tricks that I do not know. I hope to look into this more when I go home and get a better understanding of this. Afterwards I continued to study the motion of the particle and work on its plot in Mathematica. I showed this plot to Saito-sensei and he explained the physical interpretation of the graph which dealt with beats and periods. he also went on about the imaginary part of damped oscillations. I have some homework for tonight on the imaginary solution as the imaginary part is key to understanding the transformation of heat in these situations. ** June 6: [#e24a9280] Today I arrived at the lab early and worked on solving the damped oscillator equation and plotting the solution in mathematica. I successfully plotted the solution and was able manipulate the variables, so I could study the physical event. However, Saito-sensei would later explain this graph in more detail giving me a much better understanding. I meant with Thomas-san and he had me solve the 3-D Laplacian with boundary equations. This problem worked into a real ugly system of equations, but I understood the concept behind it. I continued to reteach myself some differential equations with Thomas-san and he later explained his presentation to me. He taught me about the solving the Laplacian in spherical coordinates, but I couldn't quite understand it and I plan on looking into it again. I then attended to 'group meeting' and I got to learn from Thomas-san's and Nugraha-san's presentations.I learned about the tough standards Saito-sensei has for presentations and a little about CNT's that Nugraha-san talked about. Later on, Saito-sensei explained the damped oscillator and has assigned me to study a coupled oscillator. I think I successfully modeled the non-damped homogeneous solution, but I have had some trouble solving the driven equation and plan to ask Hasdeo-san or Saito-sensei for help later. The goal is still to understand the transfer of energy through the springs. ** June 7: [#wae9b44d] Today I met with Nugraha-san and he helped me solve the driven coupled oscillator that Saito-sensei assigned me. He taught me a better trial solution to use which let me understand the phase difference. However, I think we set up our differential equations wrong and I plan on solving them again. Nugraha-san also introduced infinite coupled oscillators. We have not yet introduced a force to them, but have worked on solving the differential equations. He explained the wavevector and lent me a book to that will help explain what I learned. He assigned me to study an infinite series of different coupled masses. I plan to work on this over the weekend and show it to him when I have time. I met with Saito-sensei as well and he showed me the error in our differential equations for the coupled oscillator and introduced some concepts behind the wavevector and optical/acoustic phonons. ** June 10: [#a77f1b0a] Today,Hasdeo-san helped me finish the homework that Nugraha-san had assigned. This was to solve for the displacement of particles (coupled by springs) in an infinite lattice where the lattice cell is constructed of two particles of different mass as opposed to one. We found the solutions to the system of differential equations and he introduced the concept of zone boundaries and which phonons are acoustic or optical. I understand the math behind these conclusions, but am still rather confused by the physical meaning of frequency as a function of the wavevector. After discussing the acoustic and optical modes, Hasdeo-san introduced the model of thermal conductivity for a crystal lattice. He taught me how the transfer of energy is related to the phase shift in the two particles. I still do not completely understand this concept and plan on studying it more. After lunch I went with Thomas-san and he introduced the laplacian in polar coordinates. We had some trouble converting the initial definition, but I asked Saito-sensei for helped and he helped me work my way through. With some help from the internet, I understand the process of converting it and will begin to studying solving the differential equation. Hasdeo-san also came to my room again and discussed the thermal conductivity again and how it relates to the velocity of the center of mass. For homework, I am to study how heat transfer relates to the phase of different atoms in the atomic lattice. ** June 11: [#ib956120] Today Tatsumi-san was absent, but I managed to work on the problems that were assigned to me yesterday. First I worked on plotting the dispersion plot where frequency is a function of wavevector for a two atom periodic lattice. While I attained the right plot, I am still rather lost behind the meaning of the plot. I am still confused as to the strict definition of a phonon and how the modes are determined at various k and w values. After achieving this graph, I began to work on Hasdeo-san's homework which consisted of calculating a constant which describes heat transfer in the lattice. I managed to calculate this lattice,but am still confused as to how the relationship is attained. I managed to plot the relationship of the constant vs frequency and learned the characteristics of this relationship. However, I do have some questions regarding the imaginary part of the solution and what it means for the conductivity to approach 'negative infinity'. Next I want to attempt to solve for this constant with forces at both walls and then to work on achieving this constant for a periodic lattice. Afterwards, I met with Mizuno-san and he taught me some more about Mathematica. I learned about rules and different functions like Flatten, Transpose, Real, and Refine. These different functions can be really useful for manipulating data in Mathematica. He also helped me work on the periodic lattice problem I had been considering. I also met with Saito-sensei today and he forced me to present what I have learned. He was pretty picky, but I think I did rather well for coming up with it on the spot and barely understanding it myself. I wish I could have talked to Saito-sensei more about science, but I understand the importance of presenting what I know. For tomorrow, I am to work on a presentation of a single atom infinite lattice and learn about boltzmann's constant and heat transfer. I also want to ask Hasdeo-san about the meaning of the velocity of the center of mass and how it relates to energy. **June 12: [#s4419843] Today, I met with Mizuno-san and he continued to show me different functions of Mathematica. I learned different input methods for basic functions and how to handle matrices and vectors.I also learned how to manipulate solutions of differential equations with NDSolve. For homework, I was assigned to look up Graphics and understand them. I looked up different graphics and have began understanding the various how to use lines, polygons, circles and how to edit their style. I first want to make a solid graphic of a single oscillator and then I want to try to animate it. I hope to ask Mizuno-san for help on animating the object. Afterwards, I met with Saito-sensei and he introduced the relationship of energy to our model of springs. He showed me various laws and equations that relate energy, motion, heat, and conductivity. Today, I will study and attempt to understand these various equations. I continued to work on diffusion equations and I have a strong understanding of their derivation and physical implications. I continued to work on understand these and began to research the relationships between phonon modes, energy, and temperature. I stumbled upon the Debye and Einstein models, which Saito-sensei later encouraged me to look at. They involve some difficult math and physics I don't quite understand, but I may need some information that it has to offer. I met with Hasdeo-san to ask him various questions about the wave vector, energy equations and heat relationships. I now have a better understanding of the wave vector and the concept of phonons in solid materials. I also met with Saito-sense again and he has advised to attempt to solve the Heat equation with certain boundary conditions and to study the Debye model. ** June 13: [#zbf07e54] Today I met met with Thomas-san in the morning and we worked on solving the 2D polar laplacian. I understand what we are trying to do and can solve the equation in a general method, but am having trouble with the Fourier series solution. I borrowed Thomas-san's book on Fourier series and plan to study the solution to this differential equation. On my own time, I attempted to study some basic statistical mechanics, but am having a decent amount of trouble with the subject and plan on asking Hasdeo-san and Nugraha-san for help on it tomorrow. I particularly stuck on the derivation of Boltzman's factor and the partition function. I want to understand these topics so that I can understand the Debye Model which provides a relationship between phonon modes, heat, and energy. Hasdeo-san has provided me the assignment of learning to derive the distribution functions, but I am having trouble with this as well. I hope to ask Hasdeo-san for help tomorrow and hopefully I can receive a better understanding of the subject. ** June 14: [#r357f51e] Today I met with Nugraha-san in he morning and he introduced the expansion of coupled oscillators to containing N oscillators. He showed me how to solve for the frequency and amplitudes of movement for the oscillators. In order to solve them, I was instructed to use Mathematica. I have attempted to solve the matrix, but am having some trouble with obtaining the final solution. The program provides several solutions, that I don't quite understand. I also found the relationship of N coupled oscillators with a damping constant in the springs. This was not very difficult, but I still cannot solve the matrix in mathematica, which I will have to learn to do. ** June 17: [#cf8661e7] I spent the first part of my day discussing my project with Hasdeo-san. We talked over the basic ideas behind the project and tried developing a simple model describing energy transfer in a system of oscillators. We left the room confused as we could not understand how to relate the temperature of the wall vibrations and force of the wall. After lunch I met with Hasdeo-san again and we may have created a simple way to observe heat transfer. This method consists of comparing the amplitudes of vibration to energy for specific 'points'. These points are atoms in the lattice which may provide a general structure for a plot fit. We plan to use the general solution to the heat equation to fit the data. I will want to use Mathematica to generate the exact solution. I will begin working on this simple model right now and plan to extend this model to multiple atoms later. I also want to consider using phonons in the model in time. I met with Thomas-san today and we worked on general solutions and fourier series of the heat equation. ** June 18: [#oe7b407b] Tatsumi-san wasn't here to lecture me, so I worked on my project on my own. My initial idea with Hasdeo-san has hit a road block as the 'spring constant' appears in the relationship of the amplitudes. I do not know if there is a way to solve for this constant, or any assumptions I can make to eliminate the constant. I will have to talk with Hasdeo-san later to see if we can find a way to complete the relationship between the amplitudes of movement. I am also interested in determining the frequency of movement of the wall. We are assuming the wall vibrates at a set frequency, but is there anyway to determine frequency? I also studied the heat equation and using Fourier series to find a solution. I also continued to study phonons and have began looking into the quantum nature of phonons, but this topic seems to be to difficult to allow easy application in our model. In the afternoon, Mizuno-san taught me more about mathematica. He showed me how to work with graphics and animations and we created a small animation of the 'balls between springs'. This animation is not that great, but its a start! I think I understand the process very well and he has also taught me several functions as well. I also continued to work on my project, and talked to Hasdeo-san about several problems I had. I think we have a running model, but I will have to input values for the 'spring constant' and the 'frequency of vibration'. I also found errors in some of my derivations, so I will have to find these relationships again. I plan to work on finishing a plot tomorrow and to continue expanding the idea to many atoms. ** June 19: [#f92b7927] This morning I meant with Mizuno-san and we continued to work on the animation of the atoms between springs. We managed to get a realistic animation for the case with several atoms, but are now trying to extend the model to multiple atoms. After showing me how to use Tables, and Block, Mizuno-san assigned me to generate a code that would produce the animation for 'n' atoms. I've been working on this for a while and have hit a problem where I cannot assign the equations of motion to the animated objects. I will need to ask Mizuno-san for help later on. I spent the rest of my day working on my project. I used mathematica to try to solve for 'lambda' which represents the thermal diffusivity of the solid. I managed to solve for constants representing the general solution of the heat equation for my two atom system, but lambda is both negative and imaginary. I do not know if this makes physical sense and will need to ask Hasdeo-san for his opinion on the matter. Hasdeo-san also showed me some software on flex that may help fitting a plot to the points I obtained. ** June 20: [#keffb9a7] This morning I began looking into Fourier transforms as they keep appearing in the quantum definition of phonons. I met with Thomas-san and we worked through some differential equations using Fourier series and we derived the fourier transform for the Gaussian distribution. I understand that Fourier transforms are used to make a time dependent function, only dependent on the frequency of a periodic function, but I need to look into the derivation of fourier transforms and their application physics. However, Thomas-san showed me some interesting applications that Fourier Transforms have in imaging. I spent the afternoon continuing to work on my project. Hasdeo-san and I came to the conclusion that Boltzmann's constant was to small and affecting the calculations on my computer, so we placed Boltzmann's constant=1. This should hopefully keep the values on a larger magnitude which makes it easier to relate. I used this value and obtained a much better list plot of the three points which seems to make physical sense, but I cannot find an acceptable best fit line using the general solution to the heat equation. I keep receiving graphs whose data is way off point. I think I may need may points and I may have to take the time dependence of x along with time to use more points for a better fit. However, I would not know how to fit this into Mathematica and may need to ask for help. I believe that these points are correct, but the general solution of the heat solution is not the correct line that we should be fitting to. Hasdeo-san also began teaching me about Carbon Nanotubes and he reviewed the basic structure of the CNT and began teaching me about energy dispersion and the relationship to the brillioun zone. I hope to learn more about the carbon nanotube and to understand phonons in the CNT so I can later model it. ** June 21: [#t0eb9e69] Today Nugraha-san was busy with his thesis, so I continued to study various things on my own. I continued to study the brillioun zone and phonons in CNT's. I have begun to understand the process of analyzing phonons in CNT's but am a little lost when the math shows up. I will need to ask Hasdeo-san to explain some details to me. Nugraha-san also brought some sheets for me to read that taught me about modes of vibration in solid materials. In the afternoon, I met with Hasdeo-san and we continued to discuss the project. I told him how I had doubt in our heat solution and we checked it to find out the solution was in fact wrong. We found a new solution that I would use for a fit later. When I worked on my own, the fit worked much better and realistic lambda had been found! ** June 24: [#f62b8cf4] Today I met with Hasdeo-san and we continued to talk about my project. I showed him the lambda that I had calculated using our new heat solution. We continued to analyze my mathematica program to make sure all details were taken accounted for and we discussed certain points of the model. These points included talking about 'complex temperature' and 'complex energy' which may not be reasonable. We also noticed that 'lambda' (our constant) changed with time. Since this constant may be dynamic, Hasdeo-san has assigned me to find the 'constant's' dependence on time. I will do this by trying to find multiple values at various t's and plotting for a dependence. I hope that we can find a simple relationship between the values. I have also been assigned to work on a presentation and to begin studying phonons in CNT's and graphene.