Quantitative Derivation of Effective Evolution Equations for the Dynamics of Bose-Einstein Condensates

This thesis proves certain results concerning an important question in non-equilibrium quantum statistical mechanics which is the derivation of effective evolution equations approximating the dynamics of a system of large number of bosons initially at equilibrium (ground state at very low temperatures). The dynamics of such systems are governed by the time-dependent linear many-body Schroedinger equation from which it is typically difficult to extract useful information due to the number of particles being large. We will study quantitatively (i.e. with explicit bounds on the error) how a suitable one particle non-linear Schroedinger equation arises in the mean field limit as number of particles N → ∞ and how the appropriate corrections to the mean field will provide better approximations of the exact dynamics. In the first part of this thesis we consider the evolution of N bosons, where N is large, with two-body interactions of the form N³ᵝv(Nᵝ⋅), 0≤β≤1. The parameter β measures the strength and the range of interactions. We compare the exact evolution with an approximation which considers the evolution of a mean field coupled with an appropriate description of pair excitations, see [18,19] by Grillakis-Machedon-Margetis. We extend the results for 0 ≤ β < 1/3 in [19, 20] to the case of β < 1/2 and obtain an error bound of the form p(t)/Nᵅ, where α>0 and p(t) is a polynomial, which implies a specific rate of convergence as N → ∞. In the second part, utilizing estimates of the type discussed in the first part, we compare the exact evolution with the mean field approximation in the sense of marginals. We prove that the exact evolution is close to the approximate in trace norm for times of the order o(1)√N compared to log(o(1)N) as obtained in Chen-Lee-Schlein [6] for the Hartree evolution. Estimates of similar type are obtained for stronger interactions as well.

Local Field Potentials in the Gustatory Cortex Carry Taste Information

It has been recently shownthat localfield potentials (LFPs)fromthe auditory and visual cortices carry information about sensory stimuli, but whether this is a universal property of sensory cortices remains to be determined. Moreover, little is known about the temporal dynamics of sensory information contained in LFPs following stimulus onset. Here we investigated the time course of the amount of stimulus information in LFPs and spikes from the gustatory cortex of awake rats subjected to tastants and water delivery on the tongue. We found that the phase and amplitude of multiple LFP frequencies carry information about stimuli, which have specific time courses after stimulus delivery. The information carried by LFP phase and amplitude was independent within frequency bands, since the joint information exhibited neither synergy nor redundancy. Tastant information in LFPs was also independent and had a different time course from the information carried by spikes. These findings support the hypothesis that the brain uses different frequency channels to dynamically code for multiple features of a stimulus.

The unsteady flow field about a right circular cone in unsteady flight.

"These studies were conducted by the General Electric Company, Reentry Systems Department, for the Stability and Control Section of the Flight Dynamics Laboratory of the Air Force Research and Technology Division."

PHYSICAL FACTORS IN B CELL ACTIN DYNAMICS AND ACTIVATION

Cells adapt to their changing world by sensing environmental cues and responding appropriately. This is made possible by complex cascades of biochemical signals that originate at the cell membrane. In the last decade it has become apparent that the origin of these signals can also arise from physical cues in the environment. Our motivation is to investigate the role of physical factors in the cellular response of the B lymphocyte. B cells patrol the body for signs of invading pathogens in the form of antigen on the surface of antigen presenting cells. Binding of antigen with surface proteins initiates biochemical signaling essential to the immune response. Once contact is made, the B cell spreads on the surface of the antigen presenting cell in order to gather as much antigen as possible. The physical mechanisms that govern this process are unexplored. In this research, we examine the role of the physical parameters of antigen mobility and cell surface topography on B cell spreading and activation. Both physical parameters are biologically relevant as immunogens for vaccine design, which can provide laterally mobile and immobile antigens and topographical surfaces. Another physical parameter that influences B cell response and the formation of the cell-cell junction is surface topography. This is biologically relevant as antigen presenting cells have highly convoluted membranes, resulting in variable topography. We found that B cell activation required the formation of antigen-receptor clusters and their translocation within the attachment plane. We showed that cells which failed to achieve these mobile clusters due to prohibited ligand mobility were much less activation competent. To investigate the effect of topography, we use nano- and micro-patterned substrates, on which B cells were allowed to spread and become activated. We found that B cell spreading, actin dynamics, B cell receptor distribution and calcium signaling are dependent on the topographical patterning of the substrate. A quantitative understanding of cellular response to physical parameters is essential to uncover the fundamental mechanisms that drive B cell activation. The results of this research are highly applicable to the field of vaccine development and therapies for autoimmune diseases. Our studies of the physical aspects of lymphocyte activation will reveal the role these factors play in immunity, thus enabling their optimization for biological function and potentially enabling the production of more effective vaccines.

Large-scale neural dynamics: simple and complex

We review the use of neural field models for modelling the brain at the large scales necessary for interpreting EEG, fMRI, MEG and optical imaging data. Albeit a framework that is limited to coarse-grained or mean-field activity, neural field models provide a framework for unifying data from different imaging modalities. Starting with a description of neural mass models we build to spatially extended cortical models of layered two-dimensional sheets with long range axonal connections mediating synaptic interactions. Reformulations of the fundamental non-local mathematical model in terms of more familiar local differential (brain wave) equations are described. Techniques for the analysis of such models, including how to determine the onset of spatio-temporal pattern forming instabilities, are reviewed. Extensions of the basic formalism to treat refractoriness, adaptive feedback and inhomogeneous connectivity are described along with open challenges for the development of multi-scale models that can integrate macroscopic models at large spatial scales with models at the microscopic scale.

The dynamics of Brazilian protozoology over the past century

Brazilian scientists have been contributing to the protozoology field for more than 100 years with important discoveries of new species such as Trypanosoma cruzi and Leishmania spp. In this work, we used a Brazilian thesis database (Coordination for the Improvement of Higher Education Personnel) covering the period from 1987-2011 to identify researchers who contributed substantially to protozoology. We selected 248 advisors by filtering to obtain researchers who supervised at least 10 theses. Based on a computational analysis of the thesis databases, we found students who were supervised by these scientists. A computational procedure was developed to determine the advisors’ scientific ancestors using the Lattes Platform. These analyses provided a list of 1,997 researchers who were inspected through Lattes CV examination and allowed the identification of the pioneers of Brazilian protozoology. Moreover, we investigated the areas in which researchers who earned PhDs in protozoology are now working. We found that 68.4% of them are still in protozoology, while 16.7% have migrated to other fields. We observed that support for protozoology by national or international agencies is clearly correlated with the increase of scientists in the field. Finally, we described the academic genealogy of Brazilian protozoology by formalising the “forest” of Brazilian scientists involved in the study of protozoa and their vectors over the past century.

Guided Migration and Collective Behavior: Cell Dynamics during Cancer Progression

This dissertation focuses on gaining understanding of cell migration and collective behavior through a combination of experiment, analysis, and modeling techniques. Cell migration is a ubiquitous process that plays an important role during embryonic development and wound healing as well as in diseases like cancer, which is a particular focus of this work. As cancer cells become increasingly malignant, they acquire the ability to migrate away from the primary tumor and spread throughout the body to form metastatic tumors. During this process, changes in gene expression and the surrounding tumor environment can lead to changes in cell migration characteristics. In this thesis, I analyze how cells are guided by the texture of their environment and how cells cooperate with their neighbors to move collectively. The emergent properties of collectively moving groups are a particular focus of this work as collective cell dynamics are known to change in diseases such as cancer. The internal machinery for cell migration involves polymerization of the actin cytoskeleton to create protrusions that---in coordination with retraction of the rear of the cell---lead to cell motion. This actin machinery has been previously shown to respond to the topography of the surrounding surface, leading to guided migration of amoeboid cells. Here we show that epithelial cells on nanoscale ridge structures also show changes in the morphology of their cytoskeletons; actin is found to align with the ridge structures. The migration of the cells is also guided preferentially along the ridge length. These ridge structures are on length scales similar to those found in tumor microenvironments and as such provide a system for studying the response of the cells' internal migration machinery to physiologically relevant topographical cues. In addition to sensing surface topography, individual cells can also be influenced by the pushing and pulling of neighboring cells. The emergent properties of collectively migrating cells show interesting dynamics and are relevant for cancer progression, but have been less studied than the motion of individual cells. We use Particle Image Velocimetry (PIV) to extract the motion of a collectively migrating cell sheet from time lapse images. The resulting flow fields allow us to analyze collective behavior over multiple length and time scales. To analyze the connection between individual cell properties and collective migration behavior, we compare experimental flow fields with the migration of simulated cell groups. Our collective migration metrics allow for a quantitative comparison between experimental and simulated results. This comparison shows that tissue-scale decreases in collective behavior can result from changes in individual cell activity without the need to postulate the existence of subpopulations of leader cells or global gradients. In addition to tissue-scale trends in collective behavior, the migration of cell groups includes localized dynamic features such as cell rearrangements. An individual cell may smoothly follow the motion of its neighbors (affine motion) or move in a more individualistic manner (non-affine motion). By decomposing individual motion into both affine and non-affine components, we measure cell rearrangements within a collective sheet. Finally, finite-time Lyapunov exponent (FTLE) values capture the stretching of the flow field and reflect its chaotic character. Applying collective migration analysis techniques to experimental data on both malignant and non-malignant human breast epithelial cells reveals differences in collective behavior that are not found from analyzing migration speeds alone. Non-malignant cells show increased cooperative motion on long time scales whereas malignant cells remain uncooperative as time progresses. Combining multiple analysis techniques also shows that these two cell types differ in their response to a perturbation of cell-cell adhesion through the molecule E-cadherin. Non-malignant MCF10A cells use E-cadherin for short time coordination of collective motion, yet even with decreased E-cadherin expression, the cells remain coordinated over long time scales. In contrast, the migration behavior of malignant and invasive MCF10CA1a cells, which already shows decreased collective dynamics on both time scales, is insensitive to the change in E-cadherin expression.

PHYTOPLANKTON AND NUTRIENT DYNAMICS WITH A FOCUS ON NITROGEN FORM IN THE ANACOSTIA RIVER, IN WASHINGTON, D.C. AND WEST LAKE, IN HANGZHOU, CHINA

Nutrient loading has been linked with severe water quality impairment, ranging from hypoxia to increased frequency of harmful algal blooms (HABs), loss of fisheries, and changes in biodiversity. Waters around the globe are experiencing deleterious effects of eutrophication; however, the relative amount of nitrogen (N) and phosphorus (P) reaching these waters is not changing proportionately, with high N loads increasingly enriched in chemically-reduced N forms. Research involving two urban freshwater and nutrient enriched systems, the Anacostia River, USA, a tributary of the Potomac River feeding into the Chesapeake Bay, and West Lake, Hangzhou, Zhejiang Province, China, was conducted to assess the response of phytoplankton communities to changing N-form and N/P-ratios. Field observations involving the characterization of ambient phytoplankton communities and N-forms, as well as experimental (nutrient enrichment) manipulations were used to understand shifts in phytoplankton community composition with increasing NH4+ loads. In both locations, a >2-fold increase in ambient NH4+:NO3- ratios was followed by a shift in the phytoplankton community, with diatoms giving way to chlorophytes and cyanobacteria. Enrichment experiments mirrored this, in that samples enriched with NH4+ lead to increased abundance of chlorophytes and cyanobacteria. This work shows that in both of these systems experiencing nutrient enrichment that NH4+ supports communities dominated by more chlorophytes and cyanobacteria than other phytoplankton groups.

Guided Migration and Collective Behavior: Cell Dynamics during Cancer Progression

This dissertation focuses on gaining understanding of cell migration and collective behavior through a combination of experiment, analysis, and modeling techniques. Cell migration is a ubiquitous process that plays an important role during embryonic development and wound healing as well as in diseases like cancer, which is a particular focus of this work. As cancer cells become increasingly malignant, they acquire the ability to migrate away from the primary tumor and spread throughout the body to form metastatic tumors. During this process, changes in gene expression and the surrounding tumor environment can lead to changes in cell migration characteristics. In this thesis, I analyze how cells are guided by the texture of their environment and how cells cooperate with their neighbors to move collectively. The emergent properties of collectively moving groups are a particular focus of this work as collective cell dynamics are known to change in diseases such as cancer. The internal machinery for cell migration involves polymerization of the actin cytoskeleton to create protrusions that---in coordination with retraction of the rear of the cell---lead to cell motion. This actin machinery has been previously shown to respond to the topography of the surrounding surface, leading to guided migration of amoeboid cells. Here we show that epithelial cells on nanoscale ridge structures also show changes in the morphology of their cytoskeletons; actin is found to align with the ridge structures. The migration of the cells is also guided preferentially along the ridge length. These ridge structures are on length scales similar to those found in tumor microenvironments and as such provide a system for studying the response of the cells' internal migration machinery to physiologically relevant topographical cues. In addition to sensing surface topography, individual cells can also be influenced by the pushing and pulling of neighboring cells. The emergent properties of collectively migrating cells show interesting dynamics and are relevant for cancer progression, but have been less studied than the motion of individual cells. We use Particle Image Velocimetry (PIV) to extract the motion of a collectively migrating cell sheet from time lapse images. The resulting flow fields allow us to analyze collective behavior over multiple length and time scales. To analyze the connection between individual cell properties and collective migration behavior, we compare experimental flow fields with the migration of simulated cell groups. Our collective migration metrics allow for a quantitative comparison between experimental and simulated results. This comparison shows that tissue-scale decreases in collective behavior can result from changes in individual cell activity without the need to postulate the existence of subpopulations of leader cells or global gradients. In addition to tissue-scale trends in collective behavior, the migration of cell groups includes localized dynamic features such as cell rearrangements. An individual cell may smoothly follow the motion of its neighbors (affine motion) or move in a more individualistic manner (non-affine motion). By decomposing individual motion into both affine and non-affine components, we measure cell rearrangements within a collective sheet. Finally, finite-time Lyapunov exponent (FTLE) values capture the stretching of the flow field and reflect its chaotic character. Applying collective migration analysis techniques to experimental data on both malignant and non-malignant human breast epithelial cells reveals differences in collective behavior that are not found from analyzing migration speeds alone. Non-malignant cells show increased cooperative motion on long time scales whereas malignant cells remain uncooperative as time progresses. Combining multiple analysis techniques also shows that these two cell types differ in their response to a perturbation of cell-cell adhesion through the molecule E-cadherin. Non-malignant MCF10A cells use E-cadherin for short time coordination of collective motion, yet even with decreased E-cadherin expression, the cells remain coordinated over long time scales. In contrast, the migration behavior of malignant and invasive MCF10CA1a cells, which already shows decreased collective dynamics on both time scales, is insensitive to the change in E-cadherin expression.

Whole-cell-analysis of live cardiomyocytes using wide-field interferometric phase microscopy.

We apply wide-field interferometric microscopy techniques to acquire quantitative phase profiles of ventricular cardiomyocytes in vitro during their rapid contraction with high temporal and spatial resolution. The whole-cell phase profiles are analyzed to yield valuable quantitative parameters characterizing the cell dynamics, without the need to decouple thickness from refractive index differences. Our experimental results verify that these new parameters can be used with wide field interferometric microscopy to discriminate the modulation of cardiomyocyte contraction dynamics due to temperature variation. To demonstrate the necessity of the proposed numerical analysis for cardiomyocytes, we present confocal dual-fluorescence-channel microscopy results which show that the rapid motion of the cell organelles during contraction preclude assuming a homogenous refractive index over the entire cell contents, or using multiple-exposure or scanning microscopy.

Chemical Reaction Dynamics of the Liquid/Vapor Interface Studied by Mass Spectrometry

This thesis presents investigations of chemical reactions occurring at the liquid/vapor interface studied using novel sampling methodologies coupled with detection by mass spectrometry. Chapters 2 and 3 utilize the recently developed technique of field-induced droplet ionization mass spectrometry (FIDI-MS), in which the application of a strong electric field to a pendant microliter droplet results in the ejection of highly charged progeny droplets from the liquid surface. In Chapter 2, this method is employed to study the base-catalyzed dissociation of a surfactant molecule at the liquid/vapor interface upon uptake of ammonia from the gas phase. This process is observed to occur without significant modulation of the bulk solution pH, suggesting a transient increase in surface pH following the uptake of gaseous ammonia. Chapter 3 presents real-time studies of the oxidation of the model tropospheric organic compound glycolaldehyde by photodissociation of iron (III) oxalate complexes. The oxidation products of glycolaldehyde formed in this process are identified, and experiments in a deoxygenated environment identify the role of oxygen in the oxidation pathway and in the regeneration of iron (III) following photo-initiated reduction. Chapter 4 explores alternative methods for the study of heterogeneous reaction processes by mass spectrometric sampling from liquid surfaces. Bursting bubble ionization (BBI) and interfacial sampling with an acoustic transducer (ISAT) generate nanoliter droplets from a liquid surface that can be sampled via the atmospheric pressure interface of a mass spectrometer. Experiments on the oxidation of oleic acid by ozone using ISAT are also presented. Chapters 5 and 6 detail mechanistic studies and applications of free-radical-initiated peptide sequencing (FRIPS), a technique employing gas-phase free radical chemistry to the sequencing of peptides and proteins by mass spectrometry. Chapter 5 presents experimental and theoretical studies on the anomalous mechanism of dissociation observed in the presence of serine and threonine residues in peptides. Chapter 6 demonstrates the combination of FRIPS with ion mobility-mass spectrometry (IM-MS) for the separation of isomeric peptides.

Disorder driven transitions in non-equilibrium quantum systems

This thesis presents studies of the role of disorder in non-equilibrium quantum systems. The quantum states relevant to dynamics in these systems are very different from the ground state of the Hamiltonian. Two distinct systems are studied, (i) periodically driven Hamiltonians in two dimensions, and (ii) electrons in a one-dimensional lattice with power-law decaying hopping amplitudes. In the first system, the novel phases that are induced from the interplay of periodic driving, topology and disorder are studied. In the second system, the Anderson transition in all the eigenstates of the Hamiltonian are studied, as a function of the power-law exponent of the hopping amplitude.

In periodically driven systems the study focuses on the effect of disorder in the nature of the topology of the steady states. First, we investigate the robustness to disorder of Floquet topological insulators (FTIs) occurring in semiconductor quantum wells. Such FTIs are generated by resonantly driving a transition between the valence and conduction band. We show that when disorder is added, the topological nature of such FTIs persists as long as there is a gap at the resonant quasienergy. For strong enough disorder, this gap closes and all the states become localized as the system undergoes a transition to a trivial insulator.

Interestingly, the effects of disorder are not necessarily adverse, disorder can also induce a transition from a trivial to a topological system, thereby establishing a Floquet Topological Anderson Insulator (FTAI). Such a state would be a dynamical realization of the topological Anderson insulator. We identify the conditions on the driving field necessary for observing such a transition. We realize such a disorder induced topological Floquet spectrum in the driven honeycomb lattice and quantum well models.

Finally, we show that two-dimensional periodically driven quantum systems with spatial disorder admit a unique topological phase, which we call the anomalous Floquet-Anderson insulator (AFAI). The AFAI is characterized by a quasienergy spectrum featuring chiral edge modes coexisting with a fully localized bulk. Such a spectrum is impossible for a time-independent, local Hamiltonian. These unique characteristics of the AFAI give rise to a new topologically protected nonequilibrium transport phenomenon: quantized, yet nonadiabatic, charge pumping. We identify the topological invariants that distinguish the AFAI from a trivial, fully localized phase, and show that the two phases are separated by a phase transition.

The thesis also present the study of disordered systems using Wegner's Flow equations. The Flow Equation Method was proposed as a technique for studying excited states in an interacting system in one dimension. We apply this method to a one-dimensional tight binding problem with power-law decaying hoppings. This model presents a transition as a function of the exponent of the decay. It is shown that the the entire phase diagram, i.e. the delocalized, critical and localized phases in these systems can be studied using this technique. Based on this technique, we develop a strong-bond renormalization group that procedure where we solve the Flow Equations iteratively. This renormalization group approach provides a new framework to study the transition in this system.

MEASURING HEAD IMPACT CONTACT PRESSURE IN COLLEGIATE FOOTBALL GAMES TO CORRELATE HEAD KINEMATICS TO BRAIN KINETICS ELUCIDATING BRAIN INJURY DYNAMICS

Does a brain store thoughts and memories the way a computer saves its files? How can a single hit or a fall erase all those memories? Brain Mapping and traumatic brain injuries (TBIs) have become widely researched fields today. Many researchers have been studying TBIs caused to adult American football players however youth athletes have been rarely considered for these studies, contradicting to the fact that American football enrolls highest number of collegiate and high-school children than adults. This research is an attempt to contribute to the field of youth TBIs. Earlier studies have related head kinematics (linear and angular accelerations) to TBIs. However, fewer studies have dealt with brain kinetics (impact pressures and stresses) occurring during head-on collisions. The National Operating Committee on Standards for Athletic Equipment (NOCSAE) drop tests were conducted for linear impact accelerations and the Head Impact Contact Pressures (HICP) calculated from them were applied to a validated FE model. The results showed lateral region of the head as the most vulnerable region to damage from any drop height or impact distance followed by posterior region. The TBI tolerance levels in terms of Von-Mises and Maximum Principal Stresses deduced for lateral impact were 30 MPa and 18 MPa respectively. These levels were corresponding to 2.625 feet drop height. The drop heights beyond this value will result in TBI causing stress concentrations in human head without any detectable structural damage to the brain tissue. This data can be utilized for designing helmets that provide cushioning to brain along with providing a resistance to shear.

PARAMETRIC STUDY OF REAXFF SIMULATION PARAMETERS FOR MOLECULAR DYNAMICS MODELING OF REACTIVE CARBON GASES

Abstract The development of innovative carbon-based materials can be greatly facilitated by molecular modeling techniques. Although the Reax Force Field (ReaxFF) can be used to simulate the chemical behavior of carbon-based systems, the simulation settings required for accurate predictions have not been fully explored. Using the ReaxFF, molecular dynamics (MD) simulations are used to simulate the chemical behavior of pure carbon and hydrocarbon reactive gases that are involved in the formation of carbon structures such as graphite, buckyballs, amorphous carbon, and carbon nanotubes. It is determined that the maximum simulation time step that can be used in MD simulations with the ReaxFF is dependent on the simulated temperature and selected parameter set, as are the predicted reaction rates. It is also determined that different carbon-based reactive gases react at different rates, and that the predicted equilibrium structures are generally the same for the different ReaxFF parameter sets, except in the case of the predicted formation of large graphitic structures with the Chenoweth parameter set under specific conditions.

Vortex Shedding Dynamics in Long Aspect-Ratio Aerodynamics Bodies

The focus of the current dissertation is to study qualitatively the underlying physics of vortex-shedding and wake dynamics in long aspect-ratio aerodynamics in incompressible viscous flow through the use of the KLE method. We carried out a long series of numerical experiments in the cases of flow around the cylinder at low Reynolds numbers. The study of flow at low Reynolds numbers provides an insight in the fluid physics and also plays a critical role when applying to stalled turbine rotors. Many of the conclusions about the qualitative nature of the physical mechanisms characterizing vortex formation, shedding and further interaction analyzed here at low Re could be extended to other Re regimes and help to understand the separation of the boundary layers in airfoils and other aerodynamic surfaces. In the long run, it aims to provide a better understanding of the complex multi-physics problems involving fluid-structure-control interaction through improved mathematical computational models of the multi-physics process. Besides the scientific conclusions produced, the research work on streamlined and bluff-body condition will also serve as a valuable guide for the future design of blade aerodynamics and the placement of wind turbines and hydrakinetic turbines, increasing the efficiency in the use of expensive workforce, supplies, and infrastructure. After the introductory section describing the main fields of application of wind power and hydrokinetic turbines, we describe the main features and theoretical background of the numerical method used here. Then, we present the analysis of the numerical experimentation results for the oscillatory regime right before the onset of vortex shedding for circular cylinders. We verified the wake length of the closed near-wake behind the cylinder and analysed the decay of the wake at the wake formation region, and then studied the St-Re relationship at the Reynolds numbers before the wake sheds compared to the experimental data. We found a theoretical model that describes the time evolution of the amplitude of fluctuations in the vorticity field on the twin vortex wake, which accurately matches the numerical results in terms of the frequency of the oscillation and rate of decay. We also proposed a model based on an analog circuit that is able to interpret the concerning flow by reducing the number of degrees of freedom. It follows the idea of the non-linear oscillator and resembles the dynamics mechanism of the closed near-wake with a common configured sine wave oscillator. This low-dimensional circuital model may also help to understand the underlying physical mechanisms, related to vorticity transport, that give origin to those oscillations.