Here is a daily log for Jakob Grzesik from 2017.6.5 to 2017.8.5 in Saito Lab in Tohoku Univ.
Daily schedule (tentative) †
 09:0009:30 Begin climbing to campus
 09:3009:35 Succumb to burning feeling in calves and take a break
 09:3510:00 Resume climb to lab
 10:0011:00 Finishing/continuing any work from previous day
 11:0012:00 Discussion with Shoufiesan
 12:0013:00 Lunch
 13:0013:30 Prepare for daily report presentation
 13:3014:30 Meeting with Saitosensei
 14:3018:30 Start working on work for tomorrow's presentation and updating Pukiwiki
Goal of the project †
 Synchronization of particle motion
Questions and Answers †
This section is for posting questions from Jakobsan 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.
Q: (Placeholder) †
A: (Placeholder) †
Report †
This part is basically written by Jakobsan. Any other people can add this.
Here the information should be from new to old so that we do not need to scroll.
July 14 †
 Analytically showed how to express synchronized particles as one larger particles with the parameters of the single, larger particle just being the sum of the parameters of the smaller particles in the group.
 Finished coding and began running the synchronized group counter program.
To do:
 Run tests on random initial conditions to see if there is significant difference between the time it takes for a system to synchronize for various initial conditions.
 Come up with 3 hypotheses pertaining to parameters of the system and the time it takes for the system to syncrhonize.
 Find good threshold values for determining when a system is fully synchronized; may depend on number of oscillators.
July 13 †
 Looked into the idea of condensing a group of synchronized oscillators into a single oscillator. Allows for easy desrciption of the fully syncrhonized system.
 Also started running tests to better observe synchronization times based on phase differences to show that two particles with closer phase with synchronize quicker.
 Working on code to count how many synchronized groups are in the system at as a function of time.
To do:
 Finish the code for counting synchronized groups
July 12 †
 Presented to Konosensei and SarahSan along with Ogawasan about my research project for their lab visit.
 Found that the error in the matrix algorithm had to do with the implementation of coupled RungeKutta. The implementation to fix this error will probably make this version of the alogrithm take much longer than it's worth, so will be adjusting older, nonmatrix form of the algorithm.
To do:
 Develop an algorithm based off the old one that allows for more varied parameters and situations.
 Investigate Konosensei's questions from his visit and look at the paper on superradiance he sent, since it's about dipole synchronization.
July 11 †
 Begin retracing steps through both versions of my oscillator program, rederiving carefully the equations and steps used in the first alogrithm.
 Found that the matrix method of the algorithm produced results that seemed strange.
To do:
 Find what is the error in the matrix algorithm and how to fix it.
 Prepare for visit from Konosensei and Sarahsan.
July 10 †
 Successfully transfered data from HDD to SSD and began using new SSD.
 Advertised Nakatani program to students visiting from the University of Washington.
 Noticed that my matrix version of the program gave answers not consistent with the previous version of the algorithm. Began immediately looking through code to find potential errors.
To do:
 Rederive the formulae used in both algorithms to see if they are logically correct.
 Look through code again to see if there are any coding errors.
 Begin working on graphing function that'll let me look at previously collected data and "zoom in" on points of interest.
July 7 †
 Worked on starting tests on the lab computers.
 Met with Saitosensei and attempted to transfer data from old HDD to new SSD. Transfer was unsuccessful and will be attempted again on Monday.
To do:
 Enjoy Matsushima Bay with mother.
July 6 †
 Returned to lab and worked with Shoufiesan, Shirakurasan and Nugrahasensei on getting access to the lab computers and running tests.
 Reported to Saitosensei about the Mid Program meeting and about condition of my computer. Agreed to move data to an SSD.
To do:
 Clear up enough space on HDD to transfer data to new SSD.
 Go over time dependent force in a system and solve equations of motion using Fourier transform.
July 15 †
 Mid Program Meeting in Kyoto.
June 30 †
 Prepared for departure to Mid Program Meeting.
 Successfully setup Ubuntu on a Virtual Box and began following directions from Shirakurasan on getting access to lab computers.
To do:
 Head to Mid Program Meeting
 Work on Stanton problem.
June 29 †
 Linux install was not successful, may be related to my HDD issues.
 Shift towards using Virtual Box installation of Ubuntu.
 Added finishing touches to Mid Program presentation.
To do:
 Submit Mid Prog Presentation.
 Come up with interesting tests to look at.
 Solve the base case presented by Professor Stanton.
June 28 †
 Met with Tokyo Nakatani friends
 Continued working on Ubuntu installation, so didn't do more tests on computer.
 Left install running before I left to see if it would finish by my return.
To do:
 See if Linux installed successfully.
June 27 †
 Traveled to Tokyo
 Tried to do more testing on my laptop using the new matrix algorithm.
 Began attempting the base case of the induction approach Professor Stanton suggesting attempting.
 Began trying to setup a Ubuntu dual boot
To do:
 Come up with interesting cases to look at
 Get Ubuntu installed to begin communicating with lab computers.
June 26 †
 Tried to keep running experiments on laptop, which was slightly fixed after this weekend. I can run python and access the internet through Google Chrome, but cannot access Skype or Microsoft Office programs on my laptop
 Collected data on 64 oscillator case.
 Began talking to Kazusan and Shirakurasan about acessing flex and tube computers for later, more intensive tests. Couldn't set up much due to the state of my laptop.
 Made the matrix version of the N oscillator algorithm.
To do:
 Prepare for departure to Tokyo and mid program meeting.
 Work on fully fixing computer.
 Debug and test matrix algorithm and start using it to do more interesting and variable tests.
June 23 †
 Ran experiments with varying numbers of oscillators, with force being applied on a fraction of the oscillators.
 Found that increase in oscillators in general leads to a larger synchronization oscillation amplitude.
 Also saw that if half the oscillators had an applied force, the time it takes for full synchronization approximately doubles.
 Computer suddenly slowed down and can't open some applications now.
To do:
 Fix computer!!
 Write pseudocode for the matrix implementation.
 Take a look at system on paper to explain the results found today?
June 22 †
 Touched up on presentation
 Presented presentation to Saitosensei, who gave some more useful suggestions on how to touch it up.
 Continued reading through Pine's paper, specifically on the description of collective electron movement and the screening effect.
To do:
 Figure out effect of number of oscillators on the amplitude of position oscillation
 Apply force for a limited time and for a fraction of the total oscillators; this requires that I finish up the matrix version of my Python code so...
 Finish the matrix version of N Oscillators
 Also learn more about screening effect in electrons and its consequences.
June 21 †
 Finished working on draft for powerpoint presentation and showed to Saitosensei
 Using feedback from Saitosensei, began reworking the presentation.
 Started reeding David Pines 1952 paper on plasmons to get better idea of plasmons.
To do:
 Fix up the presentation
 Gain better understanding of Pine's paper
 Work on the matrix version of N oscilaltors RungeKutta appproximation
June 20 †
 Finished code that animated the position vs. velocity graphs for small oscillators.
 Started making presentation for Mid Program Meeting
 Began looking at Kuramoto Model
To do:
 Finish a draft of presentation for Mid Program Meeting
 Try to learn more about plasmons and phonons and the Kuramoto model.
 Begin thinking about questions to ask self and explore consequences of these questions.
June 19 †
 Coded an animation of two waves, one cosine, the other sine, using matplotlib animation on Python.
 Made code that plots position vs. velocity of a small oscillator.
To do:
 Start making presentation for Mid Program Meeting
 Figure out how to animate the position vs velocity graphs over time.
 Think about how to implement Python matrices into program in order to make it easier to come up with and implement more variations and parameter changes.
June 16 †
 Went over phase and its relation to the position and velocity of a particle.
 Discussed fully how to use Fourier Transform in solving differential equations.
 Finished coding and running RungeKutta program for 2,3,4,and 5 oscillators, with varying initial conditions and parameters.
To do:
 Plot phases of particles in phase space.
 Look into how to animate in python or some other software if necessary.
 Figure out how to implement matrix math in python to make future customization of N oscillators problem easier to implement.
 Take a look at WXMaxima computer algebra system
 Weekend: Report on the quality of Ramen Jiro, and check out Sendai Castle, while exploring other parts of the city.
June 15 †
 Learned about different ways to derive and interpret Q value, in particular, Saitosensei showed me a method using the Fourier Transform, which worked well in deriving values for overdamped and underdamped oscillators.
 Formalized and derived motions of equations for a system of N small oscillators on top of one larger oscillator
 Wrote up algorithm for doing a RungeKutta approximation for N oscillators.
To do:
 Use Fourier Transform to solve differential equations of 2nd order involving constant and nonlinear forces.
 Code the algorithm that I've written up for RungeKutta for N osciallations.
June 14 †
 Finished getting Fast Fourier Transform to work on data points for the system described by two oscillators and plotted the dominant frequencies of the system. Also played around with parameters to see how that would affect the steady state frequency and synchronization.
To do:
 Derive Q value for a critically damped system, overdamped system, underdamped system.
 Extend Python function to evaluate system for an arbitrary N oscillators.
 Think about how to measure synchronization of a system and find the phase difference between two objects in a system.
June 13 †
 Got Fast Fourier Transform to work in Python, plotting a line graph that shows the dominant frequencies of the system, whose position is approximated by the previously written RungeKutta Approximation program.
 Attended group meeting and watched Tatsumisan's presentation of his Doctorate Thesis
To do:
 Find the engineering applications of synchronization.
 Learn about Q value
To do:
June 12 †
 Finished writing a Python program to approximate and graph solutions for coupled linear differential equations using the Runge Kutta method.
 Met with Saitosensei and Shoufiesan and discussed Fourier transform
 Found an implementation in Python of Fast Fourier Transform with Shoufiesan and used it to plot already known functions in the frequency domain.
To do:
 Implement Python program on example from Thursday
 Figure out how to implement Fast Fourier Transform on discrete data set in Python.
June 9 †
 Discussed with Shoufie how to extend the RungeKutta algorithm I made for a system of a single differential equation to a system of many coupled differential equations. Learning how to better implement matrices in python looks like it'll be helpful.
 Found out that my derivation for EulerLagrange equation involving nonconservative forces had a mistake that led to different results than expected. Will be looking over and correcting it this weekend.
To do:
 Enjoy the weekend and check out the Kizuna Matsuri, which will be in town.
 Fix derivation of Euler Lagrange equation and apply it a simple situation to verify it works, and then apply it to the more complicated example I discussed with Saitosensei on Thursday.
 Find damping of a system as a function of frequency.
June 8 †
 Solved problem for HW and did some more exercises with Lagrangian mechanics and ODEs.
 Learned about an example of synchronization phenomenon with Saitosensei. Involved first getting equations of motion using Lagrangian mechanics, then applying Linear Algebra and differential equations. Ended up with system that was not analytically solvable and needed to be numerically approximated.
 Worked on deriving EulerLagrange equation in systems with nonconservative forces.
To do:
 Begin learning how to plot equations effectively. JSXGraph looks very promising, but difficult to learn.
 Design programs to numerically solve differential equation systems. Start off with Euler's method, than build upon that to do the Runge Kutta method.
 Solidify understanding of EL in nonconservative situation to the extent that I can follow and present the derivation.
June 7 †
 Walked to campus, took the right path, so it only took about 25 mins to make it to the lab.
 Plotted results for last night's HW with varying parameters.
 Practiced getting equations of motion for more complicated systems using Lagrangian Mechanics.
 Considered systems with damping forces and modelled motion using differential equations.
To do:
 Apply linear algebra and ODEs to solve a general second order differential equations problem.
 Learn how to solve ODEs and linear algebra questions using python(Start with Euler's approximation method).
June 6 †
 Walked to campus, got lost at campus. Took about 45 mins.
 Finished working on HW about action and the Euler Lagrange equation
 Worked on learning about differential equations.
 Worked on applying Lagrangian mechanics to double pendulum.
 Had lunch with Nullisan at Espace Ouvert
 First "experiment": synchronization of vegetable can movement within a box due to forces exerted while carrying it and walking at a certain pace.
To do:
 Solve equations of motion for a box on a spring with some initial external force over a fixed period of time, both before and after external force is removed.
 Numerically solve ODEs in Python(perhaps using numpy and scipy packages).
June 5 †
 Shoufiesan picked me up from Urban Castle Kawauchi, took me to campus by subway (International Center > Aobayama, 250 yen ~20 mins total)
 Nugrahasensei helped me set up lab server access, mail client, etc.
 Got a quick bento lunch on campus and ate with Saitosensei who played his ukulele
 Learned about Lagrangian mechanics, specifically in the context of a simple pendulum.
 Learned a bit about the general solutions to secondorder differential equations.
To do:
 Learn what "action" in physics is
 Derive the EulerLagrange Equation
