Rotational Brownian Dynamics simulations of clathrin cage formation
The self-assembly of nearly rigid proteins into ordered aggregates is well suited for modeling by the patchy particle approach. Patchy particles are traditionally simulated using Monte Carlo methods, to study the phase diagram, while Brownian Dynamics simulations would reveal insights into the assembly dynamics. However, Brownian Dynamics of rotating anisotropic particles gives rise to a number of complications not encountered in translational Brownian Dynamics. We thoroughly test the Rotational Brownian Dynamics scheme proposed by Naess and Elsgaeter [Macromol. Theory Simul. 13, 419 (2004); Naess and Elsgaeter Macromol. Theory Simul. 14, 300 (2005)], confirming its validity. We then apply the algorithm to simulate a patchy particle model of clathrin, a three-legged protein involved in vesicle production from lipid membranes during endocytosis. Using this algorithm we recover time scales for cage assembly comparable to those from experiments. We also briefly discuss the undulatory dynamics of the polyhedral cage
Fast simulation of Brownian dynamics in a crowded environment
Smith, Stephen; Grima, Ramon
2016-01-01
Brownian dynamics simulations are an increasingly popular tool for understanding spatially-distributed biochemical reaction systems. Recent improvements in our understanding of the cellular environment show that volume exclusion effects are fundamental to reaction networks inside cells. These systems are frequently studied by incorporating inert hard spheres (crowders) into three-dimensional Brownian dynamics simulations, however these methods are extremely slow owing to the sheer number of p...
Diffusion of Particle in Hyaluronan Solution, a Brownian Dynamics Simulation
Takasu, Masako; Tomita, Jungo
2004-04-01
Diffusion of a particle in hyaluronan solution is investigated using Brownian dynamics simulation. The slowing down of diffusion is observed, in accordance with the experimental results. The temperature dependence of the diffusion is calculated, and a turnover is obtained when the temperature is increased.
Fast simulation of Brownian dynamics in a crowded environment
Smith, Stephen
2016-01-01
Brownian dynamics simulations are an increasingly popular tool for understanding spatially-distributed biochemical reaction systems. Recent improvements in our understanding of the cellular environment show that volume exclusion effects are fundamental to reaction networks inside cells. These systems are frequently studied by incorporating inert hard spheres (crowders) into three-dimensional Brownian dynamics simulations, however these methods are extremely slow owing to the sheer number of possible collisions between particles. Here we propose a rigorous "crowder-free" method to dramatically increase simulation speed for crowded biochemical reaction systems by eliminating the need to explicitly simulate the crowders. We consider both the case where the reactive particles are point particles, and where they themselves occupy a volume. We use simulations of simple chemical reaction networks to confirm that our simplification is just as accurate as the original algorithm, and that it corresponds to a large spee...
Brownian dynamics simulations with hard-body interactions: Spherical particles
Behringer, Hans; 10.1063/1.4761827
2012-01-01
A novel approach to account for hard-body interactions in (overdamped) Brownian dynamics simulations is proposed for systems with non-vanishing force fields. The scheme exploits the analytically known transition probability for a Brownian particle on a one-dimensional half-line. The motion of a Brownian particle is decomposed into a component that is affected by hard-body interactions and into components that are unaffected. The hard-body interactions are incorporated by replacing the affected component of motion by the evolution on a half-line. It is discussed under which circumstances this approach is justified. In particular, the algorithm is developed and formulated for systems with space-fixed obstacles and for systems comprising spherical particles. The validity and justification of the algorithm is investigated numerically by looking at exemplary model systems of soft matter, namely at colloids in flow fields and at protein interactions. Furthermore, a thorough discussion of properties of other heurist...
Brownian dynamics simulation for modeling ion permeation across bionanotubes.
Krishnamurthy, Vikram; Chung, Shin-Ho
2005-03-01
The principles underlying Brownian dynamics (BD), its statistical consistency, and algorithms for practical implementation are outlined here. The ability to compute current flow across ion channels confers a distinct advantage to BD simulations compared to other simulation techniques. Thus, two obvious applications of BD ion channels are in calculation of the current-voltage and current-concentration curves, which can be directly compared to the physiological measurements to assess the reliability of the model and predictive power of the method. We illustrate how BD simulations are used to unravel the permeation dynamics in two biological ion channels-the KcsA K+ channel and CIC Cl- channel. PMID:15816176
Active microrheology of Brownian suspensions via Accelerated Stokesian Dynamics simulations
Chu, Henry; Su, Yu; Gu, Kevin; Hoh, Nicholas; Zia, Roseanna
2015-11-01
The non-equilibrium rheological response of colloidal suspensions is studied via active microrheology utilizing Accelerated Stokesian Dynamics simulations. In our recent work, we derived the theory for micro-diffusivity and suspension stress in dilute suspensions of hydrodynamically interacting colloids. This work revealed that force-induced diffusion is anisotropic, with qualitative differences between diffusion along the line of the external force and that transverse to it, and connected these effects to the role of hydrodynamic, interparticle, and Brownian forces. This work also revealed that these forces play a similar qualitative role in the anisotropy of the stress and in the evolution of the non-equilibrium osmotic pressure. Here, we show that theoretical predictions hold for suspensions ranging from dilute to near maximum packing, and for a range of flow strengths from near-equilibrium to the pure-hydrodynamic limit.
Brownian dynamics simulations of ellipsoidal magnetizable particle suspensions
Torres-Díaz, I.; Rinaldi, C.
2014-06-01
The rotational motion of soft magnetic tri-axial ellipsoidal particles suspended in a Newtonian fluid has been studied using rotational Brownian dynamics simulations by solving numerically the stochastic angular momentum equation in an orientational space described by the quaternion parameters. The model is applicable to particles where the effect of shape anisotropy is dominant. The algorithm quantifies the magnetization of a monodisperse suspension of tri-axial ellipsoids in dilute limit conditions under applied constant and time-varying magnetic fields. The variation of the relative permeability with the applied magnetic field of the particle's bulk material was included in the simulations. The results show that the equilibrium magnetization of a suspension of magnetizable tri-axial ellipsoids saturates at high magnetic field amplitudes. Additionally, the dynamic susceptibility at low magnetic field intensity presents a peak in the out-of-phase component, which is significantly smaller than the in-phase component and depends on the Langevin parameter. The dynamic magnetization of the particle suspension is in phase with the magnetic field at low and high frequencies far from the peak of the out-of-phase component.
Magnetoviscosity in dilute ferrofluids from rotational brownian dynamics simulations.
Soto-Aquino, D; Rinaldi, C
2010-10-01
Ferrofluids are suspensions of magnetic nanoparticles which respond to imposed magnetic fields by changing their viscosity without losing their fluidity. Prior work on modeling the behavior of ferrofluids has focused on using phenomenological suspension-scale continuum equations. A disadvantage of this approach is the controversy surrounding the equation describing the rate of change of the ferrofluid magnetization, the so-called magnetization relaxation equation. In this contribution the viscosity of dilute suspensions of spherical magnetic nanoparticles suspended in a Newtonian fluid and under applied shear and constant magnetic fields is studied through rotational brownian dynamics simulations. Simulation results are compared with the predictions of suspension-scale models based on three magnetization relaxation equations. Excellent agreement is observed between simulation results and the predictions of an equation due to Martsenyuk, Raikher, and Shliomis. Good qualitative agreement is observed with predictions of other equations, although these models fail to accurately predict the magnitude and shear rate dependence of the magnetic-field-dependent effective viscosity. Finally, simulation results over a wide range of conditions are collapsed into master curves using a Mason number defined based on the balance of hydrodynamic and magnetic torques. PMID:21230393
Simultaneous orthokinetic and perikinetic coagulations (SOPCs) are studied for small and large Peclet numbers (Pe) using Brownian dynamics simulation. The results demonstrate that the contributions of the Brownian motion and the shear flow to the overall coagulation rate are basically not additive. At the early stages of coagulation with small Peclet numbers, the ratio of overall coagulation rate to the rate of pure perikinetic coagulation is proportional to Pe1/2, while with high Peclet numbers, the ratio of overall coagulation rate to the rate of pure orthokinetic coagulation is proportional to Pe−1/2. Moreover, our results show that the aggregation rate generally changes with time for the SOPC, which is different from that for pure perikinetic and pure orthokinetic coagulations. By comparing the SOPC with pure perikinetic and pure orthokinetic coagulations, we show that the redistribution of particles due to Brownian motion can play a very important role in the SOPC. In addition, the effects of redistribution in the directions perpendicular and parallel to the shear flow direction are different. This perspective explains the behavior of coagulation due to the joint effects of the Brownian motion (perikinetic) and the fluid motion (orthokinetic). (electromagnetism, optics, acoustics, heat transfer, classical mechanics, and fluid dynamics)
Analysis of Brownian Dynamics Simulations of Reversible Bimolecular Reactions
Lipková, Jana
2011-01-01
A class of Brownian dynamics algorithms for stochastic reaction-diffusion models which include reversible bimolecular reactions is presented and analyzed. The method is a generalization of the λ-bcȳ model for irreversible bimolecular reactions which was introduced in [R. Erban and S. J. Chapman, Phys. Biol., 6(2009), 046001]. The formulae relating the experimentally measurable quantities (reaction rate constants and diffusion constants) with the algorithm parameters are derived. The probability of geminate recombination is also investigated. © 2011 Society for Industrial and Applied Mathematics.
Brownian dynamics in a confined geometry. Experiments and numerical simulations
Garnier, Nicolas; Ostrowsky, N.
1991-01-01
The Brownian dynamics of a colloidal suspension is measured in the immediate vicinity of a rigid surface by the Evanescent Quasielastic Light Scattering Technique. A net decrease of the measured diffusion coefficient is observed, due to the hydrodynamic slowing down of the particles very close to the wall. This effect is all the more important when the particles are allowed to get closer to the wall, i.e. when the range of the static wall/particle repulsive interaction decreases. It thus prov...
Coupling all-atom molecular dynamics simulations of ions in water with Brownian dynamics
Erban, Radek
2015-01-01
Molecular dynamics (MD) simulations of ions (K$^+$, Na$^+$, Ca$^{2+}$ and Cl$^-$) in aqueous solutions are investigated. Water is described using the SPC/E model. A stochastic coarse-grained description for ion behaviour is presented and parameterized using MD simulations. It is given as a system of coupled stochastic and ordinary differential equations, describing the ion position, velocity and acceleration. The stochastic coarse-grained model provides an intermediate description between all-atom MD simulations and Brownian dynamics (BD) models. It is used to develop a multiscale method which uses all-atom MD simulations in parts of the computational domain and (less detailed) BD simulations in the remainder of the domain.
Jung, Jiyun; Lee, Jumin; Kim, Jun Soo
2015-03-01
We present a simulation study on the mechanisms of a phase separation in dilute fluids of Lennard-Jones (LJ) particles as a model of self-interacting molecules. Molecular dynamics (MD) and Brownian dynamics (BD) simulations of the LJ fluids are employed to model the condensation of a liquid droplet in the vapor phase and the mesoscopic aggregation in the solution phase, respectively. With emphasis on the cluster growth at late times well beyond the nucleation stage, we find that the growth mechanisms can be qualitatively different: cluster diffusion and coalescence in the MD simulations and Ostwald ripening in the BD simulations. We also show that the rates of the cluster growth have distinct scaling behaviors during cluster growth. This work suggests that in the solution phase the random Brownian nature of the solute dynamics may lead to the Ostwald ripening that is qualitatively different from the cluster coalescence in the vapor phase.
A comparison of lattice-Boltzmann and Brownian dynamics simulations of dilute polymer solutions
Ladd, Tony; Kekre, Rahul; Butler, Jason
2008-11-01
We have compared lattice-Boltzmann and Brownian dynamics simulations of a single flexible polymer, in isolation and in confined geometries. In the case of the isolated chain we find agreement to within 1% in the diffusion coefficient and the Rouse mode relaxation times. We have obtained good agreement for the concentration profiles in a bounded shear flow, but the Brownian dynamics simulations currently use a superposition of the hydrodynamic fields generated by the walls. We expect to know the effects of the inter-wall correction by the time of the meeting. We have gone to some lengths to match the conditions of both simulations as closely as possible. We use identical potential parameters and correct for the differences between the periodic boundaries used in the LB simulations and the unbounded domains used in the BD simulations. We use very long runs, of the order of 10000 times the longest relaxation time, to reduce the statistical uncertainties to less than 0.1%. We find excellent agreement in the relaxation times over a wide range of temperatures and fluid viscosity. The most quantitative agreement is achieved in the weak coupling limit, where the hydrodynamic radius of the monomers is less than one quarter of the lattice spacing.
Matrix-free Brownian dynamics simulation technique for semidilute polymeric solutions
Saadat, Amir; Khomami, Bamin
2015-09-01
Evaluating the concentration dependence of static and dynamic properties of macromolecules in semidilute polymer solutions requires accurate calculation of long-range hydrodynamic interactions (HI) and short range excluded volume (EV) forces. In conventional Brownian dynamics simulations (BDS), computation of HI necessitates construction of a dense diffusion tensor commonly performed via Ewald summation. Krylov subspace techniques allow efficient decomposition of this tensor [computational cost scales as O (N2) , where N is the total number of beads in bead-spring representation of macromolecules in a simulation box] and computation of Brownian displacements in the box. In this paper, a matrix-free approach for calculation of HI is implemented which leads to O (N logN ) scaling of computational expense. The fidelity of the algorithm is demonstrated by evaluating the asymptotic value of center-of-mass diffusivity of polymer molecules at very low concentrations and their radius of gyration scaling as a function of number of beads for dilute and semidilute solutions (with concentrations up to 5 times the overlap concentration). In turn, a favorable comparison between our results and the blob theory is shown.
Chun, Myung-Suk; Kim, Chongyoup; Lee, Duck E.
2009-05-01
In our recent Brownian dynamics (BD) simulation study, the structure and dynamics of anionic polyelectrolyte xanthan in bulk solution as well as confined spaces of slitlike channel were examined by applying a coarse-grained model with nonlinear bead-spring discretization of a whole chain [J. Jeon and M.-S. Chun, J. Chem. Phys. 126, 154904 (2007)]. This model goes beyond other simulations as they did not consider both long-range electrostatic and hydrodynamic interactions between pairs of beads. Simulation parameters are obtained from the viscometric method of rheology data on the native and sonicated xanthan polysaccharides, which have a contour length less than 1μm . The size of the semiflexible polyelectrolyte can be well described by the wormlike chain model once the electrostatic effects are taken into account by the persistence length measured at a long length scale. For experimental verifications, single molecule visualization was performed on fluorescein-labeled xanthan using an inverted fluorescence microscope, and the motion of an individual molecule was quantified. Experimental results on the conformational changes in xanthan chain in the electrolyte solution have a reasonable trend to agree with the prediction by BD simulations. In the translational diffusion induced by the Debye screening effect, the simulation prediction reveals slightly higher values compared to those of our measurements, although it agrees with the literature data. Considering the experimental restrictions, our BD simulations are verified to model the single polyelectrolyte well.
DNA electrophoresis in tri-block copolymer gels--experiments and Brownian dynamics simulation
Wei, Ling; van Winkle, David H.
2015-03-01
The mobility of double-stranded DNA ladders in Pluronics®P105, P123 and F127, was measured by two-dimensional gel electrophoresis. Pluronics®are triblock copolymers which form gel-like phases of micelles arranged with cubic order at room temperature. A 10 base pair and a 25 base pair DNA ladder were used as samples in gel electrophoresis. The monotonically decreasing mobility with increasing length observed in the agarose separations is not observed in separations in Pluronics®. Rather, a complicated dependence of mobility on DNA length is observed, where mobility vs. length increases for short DNA molecules then decreases for longer molecules. There is also a variation of mobility with length correlated to the micelle diameter. Brownian dynamics simulations of a discrete wormlike chain model were performed to simulate short DNA molecules migrating in free solution and in a face-centered cubic matrix. By incorporating hydrodynamic interactions, the trend of simulated length-dependent mobility qualitatively agrees with experimental measurements.
Cook, Sara Iliafar
T, respectively). In addition to the binding strength of ssDNA nucleotide to surfaces, it is equally as important to understand the dynamics of these interactions. The force response of a simple chain-like polymeric molecule (representative of single stranded DNA) was studied using Brownian dynamics to shed light on these dynamics and the features that may be masked in SMFS experiments. Through simulations at slow peeling rates, our Brownian dynamics model confirmed the predictions of an equilibrium statistical thermodynamic model. Faster removal rates resulted in deviations from equilibrium which were dominated by a combination of Stokes (viscous) drag and a finite desorption rate of the monomeric units. Furthermore, the force probe's thermal fluctuations were shown to be affected by the spring constant of the contact mode AFM cantilever Consequently, this effect provided evidence on the source of disappearance for certain key features such as force spikes, associated with the desorption of individual links and predicted by the statistical thermodynamic model under displacement control, from SMFS experiments. In studying the elastic response of a freely jointed chain stretched in 2D and 3D, we obtained analytical expressions for two modes of stretching: i) when force is applied only to one end of the chain, and ii) when the applied force is distributed uniformly throughout the chain. By comparing, we confirmed that these expressions correctly predict the results obtained from our Brownian dynamics simulations as well as experimental results from the literature.
LI Minghai; Liu, Yongsheng; Bansil, Rama
2010-01-01
The kinetics of the transformation from the hexagonal packed cylinder (HEX) phase to the face-centered-cubic (FCC) phase was simulated using Brownian Dynamics for an ABA triblock copolymer in a selective solvent for the A block. The kinetics was obtained by instantaneously changing either the temperature of the system or the well-depth of the Lennard-Jones potential. Detailed analysis showed that the transformation occurred via a rippling mechanism. The simulation results indicated that the o...
Soto-Aquino, D; Rosso, D; Rinaldi, C
2011-11-01
Ferrofluids are colloidal suspensions of magnetic nanoparticles that exhibit normal liquid behavior in the absence of magnetic fields but respond to imposed magnetic fields by changing their viscosity without loss of fluidity. The response of ferrofluids to constant shear and magnetic fields has received a lot of attention, but the response of ferrofluids to oscillatory shear remains largely unexplored. In the present work we used rotational Brownian dynamics to study the dynamic properties of ferrofluids with thermally blocked nanoparticles under oscillatory shear and constant magnetic fields. Comparisons between simulations and modeling using the ferrohydrodynamics equations were also made. Simulation results show that, for small rotational Péclet number, the in-phase and out-of-phase components of the complex viscosity depend on the magnitude of the magnetic field and frequency of the shear, following a Maxwell-like model with field-dependent viscosity and characteristic time equal to the field-dependent transverse magnetic relaxation time of the nanoparticles. Comparison between simulations and the numerical solution of the ferrohydrodynamic equations shows that the oscillatory rotational magnetoviscosity for an oscillating shear field obtained using the kinetic magnetization relaxation equation quantitatively agrees with simulations for a wide range of Péclet number and Langevin parameter but has quantitative deviations from the simulations at high values of the Langevin parameter. These predictions indicate an apparent elastic character to the rheology of these suspensions, even though we are considering the infinitely dilute limit in which there are negligible particle-particle interactions and, as such, chains do not form. Additionally, an asymptotic analytical solution of the ferrohydrodynamics equations, valid for Peferrofluids under conditions of small shear rates. At higher shear rates the Cox-Merz rule ceases to apply. PMID:22181497
Sánchez, J H; Rinaldi, C
2009-03-15
The magnetic field dependent viscosity (magnetoviscosity) of dilute suspensions of magnetic tri-axial ellipsoidal particles suspended in a Newtonian fluid and under applied shear and magnetic fields was studied numerically. Brownian dynamics simulations were performed to compute the intrinsic magnetoviscosity of the suspension. Results are presented for the response of dilute suspensions of ellipsoidal particles to constant magnetic and shear flow fields. Suspensions of ellipsoidal particles show a significant effect of aspect ratio on the intrinsic magnetoviscosity of the suspension, and this effect is more pronounced as the aspect ratio becomes more extreme. The use of an effective rotational diffusion coefficient D(r,eff) collapses the normalized intrinsic magnetoviscosity of all suspensions to a master curve as a function of Péclet number with the Langevin parameter alpha=(mu(0)muH)/(k(B)T) as parameter, up to a critical value of alpha for which the results for suspensions of spherical particles deviate from those of suspensions of ellipsoids. This discrepancy is attributed to the action of the shear-torque on the ellipsoidal particles, which tends to orient these particles in the direction of maximum deformation of the simple shear flow, and which does not act on spherical particles. PMID:19100560
Markutsya, Sergiy; Fox, Rodney; Vigil, Dennis; Subramaniam, Shankar
2009-11-01
Nanoparticle synthesis in turbulent reactors subjects anoparticle aggregates to a homogeneous, time-varying shear flow. The shear flow results in anisotropic clusters and it is of interest to characterize the structural properties of these clusters and their effects on initiation and acceleration of aggregation, the restructuring of clusters, and their breakage. The anisotropic structure of a sheared cluster is characterized by the ratio of the major to minor axis length of the approximating ellipsoid oriented along the cluster moment of inertia tensor's principal axes. Brownian dynamics simulations show that shear flow dramatically changes the structure of aggregates by initiating the formation of more compact structures at smaller length scales perpendicular to the shear direction, and anisotropic, cigar--like structures along the shear direction. More compact clusters correspond to higher local volumetric potential energy density. Therefore, we classify the compactness and anisotropy of sheared clusters on a map of local volumetric potential energy density versus ratio of the principal values of the cluster's moment of inertia tensor. The effect of shear on breakage of clusters is characterized by the radius of gyration Rg^cr of the largest stable aggregate for a given value of the imposed steady shear rate (P'eclet number).
Brownian Dynamics Simulation of two-dimensional nanosheets under extensional flow
Xu, Yueyi; Green, Micah
2014-11-01
We investigated the morphology change of two-dimensional nanosheets under extensional flow using a coarse-grained model. Nanosheets such as graphene are promising materials for a variety of materials and electronics applications; extensional flow fields are used to cast or process liquid nanosheet dispersions in several processing techniques, including spin coating and compression molding. Process parameters, including bending stiffness and Weissenberg numbers can have a significant impact on the nanosheet morphology and the physical properties of the finished products. We use Brownian Dynamics simulations to study the impact of external flow field on a two-dimensional bead-rod lattice model. Our model was previously demonstrated for steady shear flow. Here we studied the change of morphology of graphene over time and varied the sheet size, bending stiffness and Weissenberg number. Our results showed a flattening behavior that increases with Weissenberg number. Our results also showed significant differences between nanosheets as a function of bending stiffness, with contrasting ``plate'' and ``washrag'' results under extension. The intrinsic viscosity first experiences a drop with Weissenberg number followed by a plateau associated with maximum extension.
From Molecular Dynamics to Brownian Dynamics
Erban, Radek
2014-01-01
Three coarse-grained molecular dynamics (MD) models are investigated with the aim of developing and analyzing multiscale methods which use MD simulations in parts of the computational domain and (less detailed) Brownian dynamics (BD) simulations in the remainder of the domain. The first MD model is formulated in one spatial dimension. It is based on elastic collisions of heavy molecules (e.g. proteins) with light point particles (e.g. water molecules). Two three-dimensional MD models are then investigated. The obtained results are applied to a simplified model of protein binding to receptors on the cellular membrane. It is shown that modern BD simulators of intracellular processes can be used in the bulk and accurately coupled with a (more detailed) MD model of protein binding which is used close to the membrane.
GPU accelerated Monte Carlo simulation of Brownian motors dynamics with CUDA
Spiechowicz, J; Machura, L
2014-01-01
This work presents an updated and extended guide on methods of a proper acceleration of the Monte Carlo integration of stochastic differential equations with the commonly available NVIDIA Graphics Processing Units using the CUDA programming environment. We outline the general aspects of the scientific computing on graphics cards and demonstrate them with two models of a well known phenomenon of the noise induced transport of Brownian motors in periodic structures. As a source of fluctuations in the considered systems we selected the three most commonly occurring noises: the Gaussian white noise, the white Poissonian noise and the dichotomous process also known as a random telegraph signal. The detailed discussion on various aspects of the applied numerical schemes is also presented. The measured speedup can be of the astonishing order of 2000 when compared to a typical CPU. This number significantly expands the range of problems solvable by use of stochastic simulations, allowing even an interactive research ...
An understanding of particle transport is necessary to reduce contamination of semiconductor wafers during low-pressure processing. The trajectories of particles in these reactors are determined by external forces (the most important being neutral fluid drag, thermophoresis, electrostatic, viscous ion drag, and gravitational), by Brownian motion (due to neutral and charged gas molecule collisions), and by particle inertia. Gas velocity and temperature fields are also needed for particle transport calculations, but conventional continuum fluid approximations break down at low pressures when the gas mean free path becomes comparable to chamber dimensions. Thus, in this work we use a massively parallel direct simulation Monte Carlo method to calculate low-pressure internal gas flow fields which show temperature jump and velocity slip at the reactor boundaries. Because particle residence times can be short compared to particle response times in these low-pressure systems (for which continuum diffusion theory fails), we solve the Langevin equation using a numerical Lagrangian particle tracking model which includes a fluctuating Brownian force. Because of the need for large numbers of particle trajectories to ensure statistical accuracy, the particle tracking model is also implemented on a massively parallel computer. The particle transport model is validated by comparison to the Ornstein endash Furth theoretical result for the mean square displacement of a cloud of particles. For long times, the particles tend toward a Maxwellian spatial distribution, while at short times, particle spread is controlled by their initial (Maxwellian) velocity distribution. Several simulations using these techniques are presented for particle transport and deposition in a low pressure, parallel-plate reactor geometry. The corresponding particle collection efficiencies on a wafer for different particle sizes, gas temperature gradients, and gas pressures are evaluated
Li, Minghai; Bansil, Rama
2010-01-01
The kinetics of the transformation from the hexagonal packed cylinder (HEX) phase to the face-centered-cubic (FCC) phase was simulated using Brownian Dynamics for an ABA triblock copolymer in a selective solvent for the A block. The kinetics was obtained by instantaneously changing either the temperature of the system or the well-depth of the Lennard-Jones potential. Detailed analysis showed that the transformation occurred via a rippling mechanism. The simulation results indicated that the order-order transformation (OOT) was a nucleation and growth process when the temperature of the system instantly jumped from 0.8 to 0.5. The time evolution of the structure factor obtained by Fourier Transformation showed that the peak intensities of the HEX and FCC phases could be fit well by an Avrami equation.
Mereghetti, Paolo; Wade, Rebecca C.
2012-07-26
High macromolecular concentrations are a distinguishing feature of living organisms. Understanding how the high concentration of solutes affects the dynamic properties of biological macromolecules is fundamental for the comprehension of biological processes in living systems. In this paper, we describe the implementation of mean field models of translational and rotational hydrodynamic interactions into an atomically detailed many-protein brownian dynamics simulation method. Concentrated solutions (30-40% volume fraction) of myoglobin, hemoglobin A, and sickle cell hemoglobin S were simulated, and static structure factors, oligomer formation, and translational and rotational self-diffusion coefficients were computed. Good agreement of computed properties with available experimental data was obtained. The results show the importance of both solvent mediated interactions and weak protein-protein interactions for accurately describing the dynamics and the association properties of concentrated protein solutions. Specifically, they show a qualitative difference in the translational and rotational dynamics of the systems studied. Although the translational diffusion coefficient is controlled by macromolecular shape and hydrodynamic interactions, the rotational diffusion coefficient is affected by macromolecular shape, direct intermolecular interactions, and both translational and rotational hydrodynamic interactions.
Theers, Mario; Westphal, Elmar; Gompper, Gerhard; Winkler, Roland G.
2016-03-01
The friction and diffusion coefficients of rigid spherical colloidal particles dissolved in a fluid are determined from velocity and force autocorrelation functions by mesoscale hydrodynamic simulations. Colloids with both slip and no-slip boundary conditions are considered, which are embedded in fluids modeled by multiparticle collision dynamics with and without angular momentum conservation. For no-slip boundary conditions, hydrodynamics yields the well-known Stokes law, while for slip boundary conditions the lack of angular momentum conservation leads to a reduction of the hydrodynamic friction coefficient compared to the classical result. The colloid diffusion coefficient is determined by integration of the velocity autocorrelation function, where the numerical result at shorter times is combined with the theoretical hydrodynamic expression for longer times. The suitability of this approach is confirmed by simulations of sedimenting colloids. In general, we find only minor deviations from the Stokes-Einstein relation, which even disappear for larger colloids. Importantly, for colloids with slip boundary conditions, our simulation results contradict the frequently assumed additivity of local and hydrodynamic diffusion coefficients.
Mohammadi, Maziar; Larson, Eric D.; Liu, Jun; Larson, Ronald G.
2015-01-01
Brownian dynamics simulations are performed to study the binding kinetics in the dilute-sphere limit by considering interactions of two spheres under shear flow across the entire range of Peclet numbers, spanning both perikinetic (diffusion-controlled) and orthokinetic (flow-controlled) coagulation regimes. The dilute regime is attained by carrying out two-sphere simulations in periodic boxes of different sizes and aspect ratios and extrapolating toward the infinite box limit. Effects of particle type (Janus and isotropic particles), shear rate, hydrodynamic interactions, and inter-particle potential are explored. We find that rectangular boxes with appropriate aspect ratios overcome a particle "shadow effect" that cannot be overcome with cubic boxes unless huge boxes are used. With rectangular boxes, we obtain converged binding kinetics for the whole Peclet number range, while cubic boxes of increasing size allow converged results only in the absence of flow. We consider the effect of binding both in a secondary minimum controlled by a combination of electrostatic repulsion and depletion attraction, as well as in a primary minimum governed by induced-dipole attraction. Results are computed using both realistic interaction potentials and by replacing the potential with a simple cutoff gap distance at which binding is deemed to occur. Results agree with several existing reports including Smoluchowski predictions in the zero- and infinite-shear-rate limits, and high-Pe perturbation results of Feke and Schowalter [J. Fluid Mech. 133, 17-35 (1983)] at Peclet numbers (Pe) above 100. Finally, we compute binding times for anisotropic Janus particles which have both repulsive and attractive faces, for a wide range of Pe number.
Brownian Dynamics of charged particles in a constant magnetic field
Hou, L J; Piel, A; Shukla, P K
2009-01-01
Numerical algorithms are proposed for simulating the Brownian dynamics of charged particles in an external magnetic field, taking into account the Brownian motion of charged particles, damping effect and the effect of magnetic field self-consistently. Performance of these algorithms is tested in terms of their accuracy and long-time stability by using a three-dimensional Brownian oscillator model with constant magnetic field. Step-by-step recipes for implementing these algorithms are given in detail. It is expected that these algorithms can be directly used to study particle dynamics in various dispersed systems in the presence of a magnetic field, including polymer solutions, colloidal suspensions and, particularly complex (dusty) plasmas. The proposed algorithms can also be used as thermostat in the usual molecular dynamics simulation in the presence of magnetic field.
Brownian shape dynamics in fission
Randrup Jørgen; Möller Peter
2013-01-01
It was recently shown that remarkably accurate fission-fragment mass distributions are obtained by treating the nuclear shape evolution as a Brownian walk on previously calculated five-dimensional potentialenergy surfaces; the current status of this novel method is described here.
Brownian shape dynamics in fission
Randrup Jørgen
2013-12-01
Full Text Available It was recently shown that remarkably accurate fission-fragment mass distributions are obtained by treating the nuclear shape evolution as a Brownian walk on previously calculated five-dimensional potentialenergy surfaces; the current status of this novel method is described here.
Babaei, Hasan; Keblinski, Pawel; Khodadadi, J. M.
2013-02-01
It has been recently demonstrated through experiments that the observed high enhancements in thermal conductivity of nanofluids are due to aggregation of nanoparticles rather than the previously stated mechanism of the Brownian motion-induced micro-convection. In this paper, we use equilibrium molecular dynamics simulations to investigate the role of micro-convection on the thermal conductivity of well-dispersed nanofluids. We show that while the individual terms in the heat current autocorrelation function associated with nanoparticle diffusion achieve significant values, these terms essentially cancel each other if correctly defined average enthalpy expressions are subtracted. Otherwise, erroneous thermal conductivity enhancements will be predicted, which are attributed to Brownian motion-induced micro-convection. Consequently, micro-convection does not contribute noticeably to the thermal conductivity and the predicted thermal conductivity enhancements are consistent with the effective medium theory.
3-d Brownian dynamics simulations of the smallest units of an active biological material
Luettmer-Strathmann, Jutta; Paudyal, Nabina; Adeli Koudehi, Maral
Motor proteins generate stress in a cytoskeletal network by walking on one strand of the network while being attached to another one. A protein walker in contact with two elements of the network may be considered the smallest unit of an active biological material. In vitro experiments, mathematical modeling and computer simulations have provided important insights into active matter on large and on very small length and time scales. However, it is still difficult to model the effects of local environment and interactions at intermediate scales. Recently, we developed a coarse-grained, three-dimensional model for a motor protein transporting cargo by walking on a substrate. In this work, we simulate a tethered motor protein pulling a substrate with elastic response. As the walker progresses, the retarding force due to the substrate tension increases until contact fails. We present simulation results for the effect of motor-protein activity on the tension in the substrate and the effect of the retarding force on the processivity of the molecular motor.
Sanchez, Jorge H. [Department of Chemical Engineering, University of Puerto Rico, Mayaguez campus, P.O. Box 9046, Mayaguez, PR 00681 (Puerto Rico); Facultad de Ingenieria Quimica, Universidad Pontificia Bolivariana, Medellin (Colombia); Rinaldi, Carlos [Department of Chemical Engineering, University of Puerto Rico, Mayaguez campus, P.O. Box 9046, Mayaguez, PR 00681 (Puerto Rico)], E-mail: crinaldi@uprm.edu
2009-10-15
The rotational Brownian motion of magnetized tri-axial ellipsoidal particles (orthotropic particles) suspended in a Newtonian fluid, in the dilute suspension limit, under applied d.c. and a.c. magnetic fields was studied using rotational Brownian dynamics simulations. The algorithm describing the change in the suspension magnetization was obtained from the stochastic angular momentum equation using the fluctuation-dissipation theorem and a quaternion formulation of orientation space. Simulation results are in agreement with the Langevin function for equilibrium magnetization and with single-exponential relaxation from equilibrium at small fields using Perrin's effective relaxation time. Dynamic susceptibilities for ellipsoidal particles of different aspect ratios were obtained from the response to oscillating magnetic fields of different frequencies and described by Debye's model for the complex susceptibility using Perrin's effective relaxation time. Simulations at high equilibrium and probe fields indicate that Perrin's effective relaxation time continues to describe relaxation from equilibrium and response to oscillating fields even beyond the small field limit.
Brownian motion on random dynamical landscapes
Suñé Simon, Marc; Sancho, José María; Lindenberg, Katja
2016-03-01
We present a study of overdamped Brownian particles moving on a random landscape of dynamic and deformable obstacles (spatio-temporal disorder). The obstacles move randomly, assemble, and dissociate following their own dynamics. This landscape may account for a soft matter or liquid environment in which large obstacles, such as macromolecules and organelles in the cytoplasm of a living cell, or colloids or polymers in a liquid, move slowly leading to crowding effects. This representation also constitutes a novel approach to the macroscopic dynamics exhibited by active matter media. We present numerical results on the transport and diffusion properties of Brownian particles under this disorder biased by a constant external force. The landscape dynamics are characterized by a Gaussian spatio-temporal correlation, with fixed time and spatial scales, and controlled obstacle concentrations.
Stochastic description of quantum Brownian dynamics
Yan, Yun-An; Shao, Jiushu
2016-08-01
Classical Brownian motion has well been investigated since the pioneering work of Einstein, which inspired mathematicians to lay the theoretical foundation of stochastic processes. A stochastic formulation for quantum dynamics of dissipative systems described by the system-plus-bath model has been developed and found many applications in chemical dynamics, spectroscopy, quantum transport, and other fields. This article provides a tutorial review of the stochastic formulation for quantum dissipative dynamics. The key idea is to decouple the interaction between the system and the bath by virtue of the Hubbard-Stratonovich transformation or Itô calculus so that the system and the bath are not directly entangled during evolution, rather they are correlated due to the complex white noises introduced. The influence of the bath on the system is thereby defined by an induced stochastic field, which leads to the stochastic Liouville equation for the system. The exact reduced density matrix can be calculated as the stochastic average in the presence of bath-induced fields. In general, the plain implementation of the stochastic formulation is only useful for short-time dynamics, but not efficient for long-time dynamics as the statistical errors go very fast. For linear and other specific systems, the stochastic Liouville equation is a good starting point to derive the master equation. For general systems with decomposable bath-induced processes, the hierarchical approach in the form of a set of deterministic equations of motion is derived based on the stochastic formulation and provides an effective means for simulating the dissipative dynamics. A combination of the stochastic simulation and the hierarchical approach is suggested to solve the zero-temperature dynamics of the spin-boson model. This scheme correctly describes the coherent-incoherent transition (Toulouse limit) at moderate dissipation and predicts a rate dynamics in the overdamped regime. Challenging problems
Multiscale Reaction-Diffusion Algorithms: PDE-Assisted Brownian Dynamics
Franz, Benjamin
2013-06-19
Two algorithms that combine Brownian dynami cs (BD) simulations with mean-field partial differential equations (PDEs) are presented. This PDE-assisted Brownian dynamics (PBD) methodology provides exact particle tracking data in parts of the domain, whilst making use of a mean-field reaction-diffusion PDE description elsewhere. The first PBD algorithm couples BD simulations with PDEs by randomly creating new particles close to the interface, which partitions the domain, and by reincorporating particles into the continuum PDE-description when they cross the interface. The second PBD algorithm introduces an overlap region, where both descriptions exist in parallel. It is shown that the overlap region is required to accurately compute variances using PBD simulations. Advantages of both PBD approaches are discussed and illustrative numerical examples are presented. © 2013 Society for Industrial and Applied Mathematics.
An exact solution to Brownian dynamics of a reversible bimolecular reaction in one dimension
Smith, Stephen; Grima, Ramon
2016-01-01
Brownian dynamics is a popular fine-grained method for simulating systems of interacting particles, such as chemical reactions. Though the method is simple to simulate, it is generally assumed that the dynamics is impossible to solve exactly and analytically, aside from some trivial systems. We here give the first exact analytical solution to a non-trivial Brownian dynamics system: the reaction $A+B\\xrightleftharpoons[]{}C$ in equilibrium in one-dimensional periodic space. The solution is a f...
Static/structural characteristics of non-covalent complexes, formed by terminally charged hyperbranched polymers and oppositely charged neutralizing linear polyelectrolytes, are examined by means of Brownian dynamics computer simulations. Excluded-volume, electrostatic and hydrodynamic interactions are taken into account in implicit solvent. Three pairs of complexes consisting of linear chains and hyperbranched molecules each bearing different molecular weight and distinctly diverse topologies are examined under conditions of varying electrostatic interactions. The findings from the present work demonstrate that through an appropriate modification of internal structure and external stimuli, key attributes of such complexes like size, shape and local density distribution, can be tuned at desired levels, rendering them promising candidates for a wide range of pertinent nanoscale applications
Dynamics and Efficiency of Brownian Rotors
Bauer, Wolfgang R
2008-01-01
Brownian rotors play an important role in biological systems and in future nano-technological applications. However the mechanisms determining their dynamics, efficiency and performance remain to be characterized. Here the F0 portion of the F-ATP synthase is considered as a paradigm of a Brownian rotor. In a generic analytical model we analyze the stochastic rotation of F0-like motors as a function of the driving free energy difference and of the free energy profile the rotor is subjected to. The latter is composed of the rotor interaction with its surroundings, of the free energy of chemical transitions, and of the workload. The dynamics and mechanical efficiency of the rotor depends on the magnitude of its stochastic motion driven by the free energy energy difference and its rectification on the reaction-diffusion path. We analyze which free energy profiles provide maximum flow and how their arrangement on the underlying reaction-diffusion path affects rectification and -- by this -- the efficiency.
Dynamical objectivity in quantum Brownian motion
Tuziemski, J.; Korbicz, J. K.
2015-11-01
Classical objectivity as a property of quantum states —a view proposed to explain the observer-independent character of our world from quantum theory, is an important step in bridging the quantum-classical gap. It was recently derived in terms of spectrum broadcast structures for small objects embedded in noisy photon-like environments. However, two fundamental problems have arisen: a description of objective motion and applicability to other types of environments. Here we derive an example of objective states of motion in quantum mechanics by showing the formation of dynamical spectrum broadcast structures in the celebrated, realistic model of decoherence —Quantum Brownian Motion. We do it for realistic, thermal environments and show their noise-robustness. This opens a potentially new method of studying the quantum-to-classical transition.
An elementary singularity-free Rotational Brownian Dynamics algorithm for anisotropic particles
Ilie, Ioana M.; Briels, Wim J. [Computational Biophysics, Faculty of Science and Technology, University of Twente, P.O. Box 217, 7500 AE Enschede (Netherlands); MESA+ Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE Enschede (Netherlands); Otter, Wouter K. den, E-mail: w.k.denotter@utwente.nl [Computational Biophysics, Faculty of Science and Technology, University of Twente, P.O. Box 217, 7500 AE Enschede (Netherlands); MESA+ Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE Enschede (Netherlands); Multi Scale Mechanics, Faculty of Engineering Technology, University of Twente, P.O. Box 217, 7500 AE Enschede (Netherlands)
2015-03-21
Brownian Dynamics is the designated technique to simulate the collective dynamics of colloidal particles suspended in a solution, e.g., the self-assembly of patchy particles. Simulating the rotational dynamics of anisotropic particles by a first-order Langevin equation, however, gives rise to a number of complications, ranging from singularities when using a set of three rotational coordinates to subtle metric and drift corrections. Here, we derive and numerically validate a quaternion-based Rotational Brownian Dynamics algorithm that handles these complications in a simple and elegant way. The extension to hydrodynamic interactions is also discussed.
An elementary singularity-free Rotational Brownian Dynamics algorithm for anisotropic particles
Brownian Dynamics is the designated technique to simulate the collective dynamics of colloidal particles suspended in a solution, e.g., the self-assembly of patchy particles. Simulating the rotational dynamics of anisotropic particles by a first-order Langevin equation, however, gives rise to a number of complications, ranging from singularities when using a set of three rotational coordinates to subtle metric and drift corrections. Here, we derive and numerically validate a quaternion-based Rotational Brownian Dynamics algorithm that handles these complications in a simple and elegant way. The extension to hydrodynamic interactions is also discussed
Greives, Nicholas; Zhou, Huan-Xiang
2012-10-01
A method developed by Northrup et al. [J. Chem. Phys. 80, 1517 (1984)], 10.1063/1.446900 for calculating protein-ligand binding rate constants (ka) from Brownian dynamics (BD) simulations has been widely used for rigid molecules. Application to flexible molecules is limited by the formidable computational cost to treat conformational fluctuations during the long BD simulations necessary for ka calculation. Here, we propose a new method called BDflex for ka calculation that circumvents this problem. The basic idea is to separate the whole space into an outer region and an inner region, and formulate ka as the product of kE and bar η _d, which are obtained by separately solving exterior and interior problems. kE is the diffusion-controlled rate constant for the ligand in the outer region to reach the dividing surface between the outer and inner regions; in this exterior problem conformational fluctuations can be neglected. bar η _d is the probability that the ligand, starting from the dividing surface, will react at the binding site rather than escape to infinity. The crucial step in reducing the determination of bar η _d to a problem confined to the inner region is a radiation boundary condition imposed on the dividing surface; the reactivity on this boundary is proportional to kE. By confining the ligand to the inner region and imposing the radiation boundary condition, we avoid multiple-crossing of the dividing surface before reaction at the binding site and hence dramatically cut down the total simulation time, making the treatment of conformational fluctuations affordable. BDflex is expected to have wide applications in problems where conformational fluctuations of the molecules are crucial for productive ligand binding, such as in cases where transient widening of a bottleneck allows the ligand to access the binding pocket, or the binding site is properly formed only after ligand entrance induces the closure of a lid.
Dynamical 3-Space: Anisotropic Brownian Motion Experiment
Cahill R. T.
2015-01-01
In 2014 Jiapei Dai reported evidence of anisotropic Brownian motion of a toluidine blue colloid solution in water. In 2015 Felix Scholkmann analysed the Dai data and detected a sidereal time dependence, indicative of a process driving the preferred Brownian mo- tion diffusion direction to a star-based preferred direction. Here we further analyse the Dai data and extract the RA and Dec of that preferred direction, and relate the data to previous determinations from NASA Spacecr...
Momentum conserving Brownian dynamics propagator for complex soft matter fluids.
Padding, J T; Briels, W J
2014-12-28
We present a Galilean invariant, momentum conserving first order Brownian dynamics scheme for coarse-grained simulations of highly frictional soft matter systems. Friction forces are taken to be with respect to moving background material. The motion of the background material is described by locally averaged velocities in the neighborhood of the dissolved coarse coordinates. The velocity variables are updated by a momentum conserving scheme. The properties of the stochastic updates are derived through the Chapman-Kolmogorov and Fokker-Planck equations for the evolution of the probability distribution of coarse-grained position and velocity variables, by requiring the equilibrium distribution to be a stationary solution. We test our new scheme on concentrated star polymer solutions and find that the transverse current and velocity time auto-correlation functions behave as expected from hydrodynamics. In particular, the velocity auto-correlation functions display a long time tail in complete agreement with hydrodynamics. PMID:25554134
From Brownian Dynamics to Markov Chain: An Ion Channel Example
Chen, Wan
2014-02-27
A discrete rate theory for multi-ion channels is presented, in which the continuous dynamics of ion diffusion is reduced to transitions between Markovian discrete states. In an open channel, the ion permeation process involves three types of events: an ion entering the channel, an ion escaping from the channel, or an ion hopping between different energy minima in the channel. The continuous dynamics leads to a hierarchy of Fokker-Planck equations, indexed by channel occupancy. From these the mean escape times and splitting probabilities (denoting from which side an ion has escaped) can be calculated. By equating these with the corresponding expressions from the Markov model, one can determine the Markovian transition rates. The theory is illustrated with a two-ion one-well channel. The stationary probability of states is compared with that from both Brownian dynamics simulation and the hierarchical Fokker-Planck equations. The conductivity of the channel is also studied, and the optimal geometry maximizing ion flux is computed. © 2014 Society for Industrial and Applied Mathematics.
Van den Broeck, C; Kawai, R
2006-06-01
Onsager symmetry implies that a Brownian motor, driven by a temperature gradient, will also perform a refrigerator function upon loading. We analytically calculate the corresponding heat flow for an exactly solvable microscopic model and compare it with molecular dynamics simulations. PMID:16803223
Studying protein assembly with reversible Brownian dynamics of patchy particles
Assembly of protein complexes like virus shells, the centriole, the nuclear pore complex, or the actin cytoskeleton is strongly determined by their spatial structure. Moreover, it is becoming increasingly clear that the reversible nature of protein assembly is also an essential element for their biological function. Here we introduce a computational approach for the Brownian dynamics of patchy particles with anisotropic assemblies and fully reversible reactions. Different particles stochastically associate and dissociate with microscopic reaction rates depending on their relative spatial positions. The translational and rotational diffusive properties of all protein complexes are evaluated on-the-fly. Because we focus on reversible assembly, we introduce a scheme which ensures detailed balance for patchy particles. We then show how the macroscopic rates follow from the microscopic ones. As an instructive example, we study the assembly of a pentameric ring structure, for which we find excellent agreement between simulation results and a macroscopic kinetic description without any adjustable parameters. This demonstrates that our approach correctly accounts for both the diffusive and reactive processes involved in protein assembly
Long, Hai; Chang, Christopher H; King, Paul W; Ghirardi, Maria L; Kim, Kwiseon
2008-10-01
The [FeFe] hydrogenase from the green alga Chlamydomonas reinhardtii can catalyze the reduction of protons to hydrogen gas using electrons supplied from photosystem I and transferred via ferredoxin. To better understand the association of the hydrogenase and the ferredoxin, we have simulated the process over multiple timescales. A Brownian dynamics simulation method gave an initial thorough sampling of the rigid-body translational and rotational phase spaces, and the resulting trajectories were used to compute the occupancy and free-energy landscapes. Several important hydrogenase-ferredoxin encounter complexes were identified from this analysis, which were then individually simulated using atomistic molecular dynamics to provide more details of the hydrogenase and ferredoxin interaction. The ferredoxin appeared to form reasonable complexes with the hydrogenase in multiple orientations, some of which were good candidates for inclusion in a transition state ensemble of configurations for electron transfer. PMID:18621810
Brownian dynamics of confined rigid bodies
Delong, Steven; Balboa Usabiaga, Florencio; Donev, Aleksandar, E-mail: donev@courant.nyu.edu [Courant Institute of Mathematical Sciences, New York University, New York, New York 10012 (United States)
2015-10-14
We introduce numerical methods for simulating the diffusive motion of rigid bodies of arbitrary shape immersed in a viscous fluid. We parameterize the orientation of the bodies using normalized quaternions, which are numerically robust, space efficient, and easy to accumulate. We construct a system of overdamped Langevin equations in the quaternion representation that accounts for hydrodynamic effects, preserves the unit-norm constraint on the quaternion, and is time reversible with respect to the Gibbs-Boltzmann distribution at equilibrium. We introduce two schemes for temporal integration of the overdamped Langevin equations of motion, one based on the Fixman midpoint method and the other based on a random finite difference approach, both of which ensure that the correct stochastic drift term is captured in a computationally efficient way. We study several examples of rigid colloidal particles diffusing near a no-slip boundary and demonstrate the importance of the choice of tracking point on the measured translational mean square displacement (MSD). We examine the average short-time as well as the long-time quasi-two-dimensional diffusion coefficient of a rigid particle sedimented near a bottom wall due to gravity. For several particle shapes, we find a choice of tracking point that makes the MSD essentially linear with time, allowing us to estimate the long-time diffusion coefficient efficiently using a Monte Carlo method. However, in general, such a special choice of tracking point does not exist, and numerical techniques for simulating long trajectories, such as the ones we introduce here, are necessary to study diffusion on long time scales.
Brownian dynamics of confined rigid bodies
Delong, Steven; Balboa Usabiaga, Florencio; Donev, Aleksandar
2015-10-01
We introduce numerical methods for simulating the diffusive motion of rigid bodies of arbitrary shape immersed in a viscous fluid. We parameterize the orientation of the bodies using normalized quaternions, which are numerically robust, space efficient, and easy to accumulate. We construct a system of overdamped Langevin equations in the quaternion representation that accounts for hydrodynamic effects, preserves the unit-norm constraint on the quaternion, and is time reversible with respect to the Gibbs-Boltzmann distribution at equilibrium. We introduce two schemes for temporal integration of the overdamped Langevin equations of motion, one based on the Fixman midpoint method and the other based on a random finite difference approach, both of which ensure that the correct stochastic drift term is captured in a computationally efficient way. We study several examples of rigid colloidal particles diffusing near a no-slip boundary and demonstrate the importance of the choice of tracking point on the measured translational mean square displacement (MSD). We examine the average short-time as well as the long-time quasi-two-dimensional diffusion coefficient of a rigid particle sedimented near a bottom wall due to gravity. For several particle shapes, we find a choice of tracking point that makes the MSD essentially linear with time, allowing us to estimate the long-time diffusion coefficient efficiently using a Monte Carlo method. However, in general, such a special choice of tracking point does not exist, and numerical techniques for simulating long trajectories, such as the ones we introduce here, are necessary to study diffusion on long time scales.
Solano, Carlos J F; Pothula, Karunakar R; Prajapati, Jigneshkumar D; De Biase, Pablo M; Noskov, Sergei Yu; Kleinekathöfer, Ulrich
2016-05-10
All-atom molecular dynamics simulations have a long history of applications studying ion and substrate permeation across biological and artificial pores. While offering unprecedented insights into the underpinning transport processes, MD simulations are limited in time-scales and ability to simulate physiological membrane potentials or asymmetric salt solutions and require substantial computational power. While several approaches to circumvent all of these limitations were developed, Brownian dynamics simulations remain an attractive option to the field. The main limitation, however, is an apparent lack of protein flexibility important for the accurate description of permeation events. In the present contribution, we report an extension of the Brownian dynamics scheme which includes conformational dynamics. To achieve this goal, the dynamics of amino-acid residues was incorporated into the many-body potential of mean force and into the Langevin equations of motion. The developed software solution, called BROMOCEA, was applied to ion transport through OmpC as a test case. Compared to fully atomistic simulations, the results show a clear improvement in the ratio of permeating anions and cations. The present tests strongly indicate that pore flexibility can enhance permeation properties which will become even more important in future applications to substrate translocation. PMID:27088446
Non-Brownian diffusion in lipid membranes: Experiments and simulations.
Metzler, R; Jeon, J-H; Cherstvy, A G
2016-10-01
The dynamics of constituents and the surface response of cellular membranes-also in connection to the binding of various particles and macromolecules to the membrane-are still a matter of controversy in the membrane biophysics community, particularly with respect to crowded membranes of living biological cells. We here put into perspective recent single particle tracking experiments in the plasma membranes of living cells and supercomputing studies of lipid bilayer model membranes with and without protein crowding. Special emphasis is put on the observation of anomalous, non-Brownian diffusion of both lipid molecules and proteins embedded in the lipid bilayer. While single component, pure lipid bilayers in simulations exhibit only transient anomalous diffusion of lipid molecules on nanosecond time scales, the persistence of anomalous diffusion becomes significantly longer ranged on the addition of disorder-through the addition of cholesterol or proteins-and on passing of the membrane lipids to the gel phase. Concurrently, experiments demonstrate the anomalous diffusion of membrane embedded proteins up to macroscopic time scales in the minute time range. Particular emphasis will be put on the physical character of the anomalous diffusion, in particular, the occurrence of ageing observed in the experiments-the effective diffusivity of the measured particles is a decreasing function of time. Moreover, we present results for the time dependent local scaling exponent of the mean squared displacement of the monitored particles. Recent results finding deviations from the commonly assumed Gaussian diffusion patterns in protein crowded membranes are reported. The properties of the displacement autocorrelation function of the lipid molecules are discussed in the light of their appropriate physical anomalous diffusion models, both for non-crowded and crowded membranes. In the last part of this review we address the upcoming field of membrane distortion by elongated membrane
Differential dynamic microscopy to characterize Brownian motion and bacteria motility
Germain, David; Leocmach, Mathieu; Gibaud, Thomas
2016-03-01
We have developed a lab module for undergraduate students, which involves the process of quantifying the dynamics of a suspension of microscopic particles using Differential Dynamic Microscopy (DDM). DDM is a relatively new technique that constitutes an alternative method to more classical techniques such as dynamic light scattering (DLS) or video particle tracking (VPT). The technique consists of imaging a particle dispersion with a standard light microscope and a camera and analyzing the images using a digital Fourier transform to obtain the intermediate scattering function, an autocorrelation function that characterizes the dynamics of the dispersion. We first illustrate DDM in the textbook case of colloids under Brownian motion, where we measure the diffusion coefficient. Then we show that DDM is a pertinent tool to characterize biological systems such as motile bacteria.
Brownian Dynamics of Colloidal Particles in Lyotropic Chromonic Liquid Crystals
Martinez, Angel; Collings, Peter J.; Yodh, Arjun G.
We employ video microscopy to study the Brownian dynamics of colloidal particles in the nematic phase of lyotropic chromonic liquid crystals (LCLCs). These LCLCs (in this case, DSCG) are water soluble, and their nematic phases are characterized by an unusually large elastic anisotropy. Our preliminary measurements of particle mean-square displacement for polystyrene colloidal particles (~5 micron-diameter) show diffusive and sub-diffusive behaviors moving parallel and perpendicular to the nematic director, respectively. In order to understand these motions, we are developing models that incorporate the relaxation of elastic distortions of the surrounding nematic field. Further experiments to confirm these preliminary results and to determine the origin of these deviations compared to simple diffusion theory are ongoing; our results will also be compared to previous diffusion experiments in nematic liquid crystals. We gratefully acknowledge financial support through NSF DMR12-05463, MRSEC DMR11-20901, and NASA NNX08AO0G.
Mériguet, G; Jardat, M; Turq, P
2004-09-22
We present Brownian dynamics simulations of real charge-stabilized ferrofluids, which are stable colloidal dispersions of magnetic nanoparticles, with and without the presence of an external magnetic field. The colloidal suspensions are treated as collections of monodisperse spherical particles, bearing point dipoles at their centers and undergoing translational and rotational Brownian motions. The overall repulsive isotropic interactions between particles, governed by electrostatic repulsions, are taken into account by a one-component effective pair interaction potential. The potential parameters are fitted in order that computed structure factors are close to the experimental ones. Two samples of ferrofluid differing by the particle diameter and consequently by the intensity of the magnetic interaction are considered here. The magnetization and birefringence curves are computed: a deviation from the ideal Langevin behaviors is observed if the dipolar moment of particles is sufficiently large. Structure factors are also computed from simulations with and without an applied magnetic field H: the microstructure of the repulsive ferrofluid becomes anisotropic under H. Even our simple modeling of the suspension allows us to account for the main experimental features: an increase of the peak intensity is observed in the direction perpendicular to the field whereas the peak intensity decreases in the direction parallel to the field. PMID:15367036
Chung, Shin-Ho; Corry, Ben
2007-01-01
In the narrow segment of an ion conducting pathway, it is likely that a permeating ion influences the positions of the nearby atoms that carry partial or full electronic charges. Here we introduce a method of incorporating the motion of charged atoms lining the pore into Brownian dynamics simulations of ion conduction. The movements of the carbonyl groups in the selectivity filter of the KcsA channel are calculated explicitly, allowing their bond lengths, bond angles, and dihedral angels to c...
Chen, Shing Bor
2015-12-01
Brownian dynamics simulation has been employed to study the dynamic behavior of particles in three-dimensional ordered porous media subject to a sinusoidal force field. The media comprises interconnected spherical cavities arranged in a simple cubic lattice. The thermal noise assists the particles to undergo cavity hopping, leading to a displacement behavior analogous to stochastic resonance, when the imposed field is strong enough but not aligned with the aperture lines, and the oscillation frequency is not too high. The periodic mean trajectory depends on the strength, frequency, and orientation of the imposed field. At sufficiently large field strength, the periodic particle displacement can become nonsinusoidal due to the strong hindrance and pinning effect of the cavity wall. PMID:26764630
Dubina, Sean Hyun; Wedgewood, Lewis Edward
2016-07-01
Ferrofluids are often favored for their ability to be remotely positioned via external magnetic fields. The behavior of particles in ferromagnetic clusters under uniformly applied magnetic fields has been computationally simulated using the Brownian dynamics, Stokesian dynamics, and Monte Carlo methods. However, few methods have been established that effectively handle the basic principles of magnetic materials, namely, Maxwell's equations. An iterative constraint method was developed to satisfy Maxwell's equations when a uniform magnetic field is imposed on ferrofluids in a heterogeneous Brownian dynamics simulation that examines the impact of ferromagnetic clusters in a mesoscale particle collection. This was accomplished by allowing a particulate system in a simple shear flow to advance by a time step under a uniformly applied magnetic field, then adjusting the ferroparticles via an iterative constraint method applied over sub-volume length scales until Maxwell's equations were satisfied. The resultant ferrofluid model with constraints demonstrates that the magnetoviscosity contribution is not as substantial when compared to homogeneous simulations that assume the material's magnetism is a direct response to the external magnetic field. This was detected across varying intensities of particle-particle interaction, Brownian motion, and shear flow. Ferroparticle aggregation was still extensively present but less so than typically observed.
Dynamics of non-Brownian fiber suspensions under periodic shear.
Franceschini, Alexandre; Filippidi, Emmanouela; Guazzelli, Elisabeth; Pine, David J
2014-09-21
We report experiments studying the dynamics of dense non-Brownian fiber suspensions subjected to periodic oscillatory shear. We find that periodic shear initially causes fibers to collide and to undergo irreversible diffusion. As time progresses, the fibers tend to orient in the vorticity direction while the number of collisions decreases. Ultimately, the system goes to one of two steady states: an absorbing steady state, where collisions cease and the fibers undergo reversible trajectories; an active state, where fibers continue to collide causing them to diffuse and undergo irreversible trajectories. Collisions between fibers can be characterized by an effective volume fraction Φ with a critical volume fraction Φc that separates absorbing from active (diffusing) steady states. The effective volume fraction Φ depends on the mean fiber orientation and thus decreases in time as fibers progressively orient under periodic shear. In the limit that the temporal evolution of Φ is slow compared to the activity relaxation time τ, all the data for all strain amplitudes and all concentrations can be scaled onto a single master curve with a functional dependence well-described by t(-β/ν)R(e(-t)R), where tR is the rescaled time. As Φ → Φc, τ diverges. Therefore, for experiments in which Φ(t) starts above Φc but goes to a steady state below Φc, departures from scaling are observed for Φ very near Φc. The critical exponents are measured to be β = 0.84 ± 0.04 and ν = 1.1 ± 0.1, which is consistent with the Manna universality class for directed percolation. PMID:25068577
Nonlinear Brownian motion - mean square displacement
W.Ebeling
2004-01-01
Full Text Available The stochastic dynamics of self-propelled Brownian particles is studied by means of the Langevin and the Fokker-Planck approach. We model the driving by a nonlinear friction function which has a negative part at small velocities, leading to active Brownian motion of the particles. The mean square displacement is estimated analytically and compared with numerical simulations.
Maldonado-Camargo, L.; Torres-Díaz, I.; Chiu-Lam, A.; Hernández, M.; Rinaldi, C.
2016-08-01
We demonstrate how dynamic magnetic susceptibility measurements (DMS) can be used to estimate the relative contributions of Brownian and Néel relaxation to the dynamic magnetic response of a magnetic fluid, a suspension of magnetic nanoparticles. The method applies to suspensions with particles that respond through Brownian or Néel relaxation and for which the characteristic Brownian and Néel relaxation times are widely separated. First, we illustrate this using magnetic fluids consisting of mixtures of particles that relax solely by the Brownian or Néel mechanisms. Then, it is shown how the same approach can be applied to estimate the relative contributions of Brownian and Néel relaxation in a suspension consisting of particles obtained from a single synthesis and whose size distribution straddles the transition from Néel to Brownian relaxation.
Molecular dynamics test of the Brownian description of Na(+) motion in water
Wilson, M. A.; Pohorille, A.; Pratt, L. R.
1985-01-01
The present paper provides the results of molecular dynamics calculations on a Na(+) ion in aqueous solution. Attention is given to the sodium-oxygen and sodium-hydrogen radial distribution functions, the velocity autocorrelation function for the Na(+) ion, the autocorrelation function of the force on the stationary ion, and the accuracy of Brownian motion assumptions which are basic to hydrodynamic models of ion dyanmics in solution. It is pointed out that the presented calculations provide accurate data for testing theories of ion dynamics in solution. The conducted tests show that it is feasible to calculate Brownian friction constants for ions in aqueous solutions. It is found that for Na(+) under the considered conditions the Brownian mobility is in error by only 60 percent.
Cosseddu, Salvatore M; Allen, Michael P; Rodger, P M; Luchinsky, Dmitry G; McClintock, Peter V E
2013-01-01
The statistical and dynamical properties of ions in the selectivity filter of the KcsA ion channel are considered on the basis of molecular dynamics (MD) simulations of the KcsA protein embedded in a lipid membrane surrounded by an ionic solution. A new approach to the derivation of a Brownian dynamics (BD) model of ion permeation through the filter is discussed, based on unbiased MD simulations. It is shown that depending on additional assumptions, ion's dynamics can be described either by under-damped Langevin equation with constant damping and white noise or by Langevin equation with a fractional memory kernel. A comparison of the potential of the mean force derived from unbiased MD simulations with the potential produced by the umbrella sampling method demonstrates significant differences in these potentials. The origin of these differences is an open question that requires further clarifications.
Cosseddu, S. M.; Khovanov, I. A.; Allen, M. P.; Rodger, P. M.; Luchinsky, D. G.; McClintock, P. V. E.
2013-10-01
The statistical and dynamical properties of ions in the selectivity filter of the KcsA ion channel are considered on the basis of molecular dynamics (MD) simulations of the KcsA protein embedded in a lipid membrane surrounded by an ionic solution. A new approach to the derivation of a Brownian dynamics (BD) model of ion permeation through the filter is discussed, based on unbiased MD simulations. It is shown that depending on additional assumptions, ion's dynamics can be described either by under-damped Langevin equation with constant damping and white noise or by Langevin equation with a fractional memory kernel. A comparison of the potential of the mean force derived from unbiased MD simulations with the potential produced by the umbrella sampling method demonstrates significant differences in these potentials. The origin of these differences is an open question that requires further clarifications.
Khan, Siddique J.
We carry out Brownian Dynamics Simulations to study the self-assembly of ligated gold nanoparticles for various ligand chain lengths. First, we develop a phenomenological model for an effective nanoparticle-nanoparticle pair potential by treating the ligands as flexible polymer chains. Besides van der Waals interactions, we incorporate both the free energy of mixing and elastic contributions from compression of the ligands in our effective pair potentials. The separation of the nanoparticles at the potential minimum compares well with experimental results of gold nanoparticle superlattice constants for various ligand lengths. Next, we use the calculated pair potentials as input to Brownian dynamics simulations for studying the formation of nanoparticle assembly in three dimensions. For dodecanethiol ligated nanoparticles in toluene, our model gives a relatively shallower well depth and the clusters formed after a temperature quench are compact in morphology. Simulation results for the kinetics of cluster growth in this case are compared with phase separations in binary mixtures. For decanethiol ligated nanoparticles, the model well depth is found to be deeper, and simulations show hybrid, fractal-like morphology for the clusters. Cluster morphology in this case shows a compact structure at short length scales and a fractal structure at large length scales. Growth kinetics for this deeper potential depth is compared with the diffusion-limited cluster-cluster aggregation (DLCA) model. We also did simulation studies of nanoparticle supercluster (NPSC) nucleation from a temperature quenched system. Induction periods are observed with times that yield a reasonable supercluster interfacial tension via classical nucleation theory (CNT). However, only the largest pre-nucleating clusters are dense and the cluster size can occasionally range greater than the critical size in the pre-nucleation regime until a cluster with low enough energy occurs, then nucleation ensues. Late
Brownian agents and active particles collective dynamics in the natural and social sciences
Schweitzer, Frank
2007-01-01
""This book lays out a vision for a coherent framework for understanding complex systems"" (from the foreword by J. Doyne Farmer). By developing the genuine idea of Brownian agents, the author combines concepts from informatics, such as multiagent systems, with approaches of statistical many-particle physics. This way, an efficient method for computer simulations of complex systems is developed which is also accessible to analytical investigations and quantitative predictions. The book demonstrates that Brownian agent models can be successfully applied in many different contexts, ranging from
Generalized Scaling and the Master Variable for Brownian Magnetic Nanoparticle Dynamics
Reeves, Daniel B.; Yipeng Shi; Weaver, John B.
2016-01-01
Understanding the dynamics of magnetic particles can help to advance several biomedical nanotechnologies. Previously, scaling relationships have been used in magnetic spectroscopy of nanoparticle Brownian motion (MSB) to measure biologically relevant properties (e.g., temperature, viscosity, bound state) surrounding nanoparticles in vivo. Those scaling relationships can be generalized with the introduction of a master variable found from non-dimensionalizing the dynamical Langevin equation. T...
Brownian Agents and Active Particles: Collective Dynamics in the Natural and Social Sciences
This is a book about the modelling of complex systems and, unlike many books on this subject, concentrates on the discussion of specific systems and gives practical methods for modelling and simulating them. This is not to say that the author does not devote space to the general philosophy and definition of complex systems and agent-based modelling, but the emphasis is definitely on the development of concrete methods for analysing them. This is, in my view, to be welcomed and I thoroughly recommend the book, especially to those with a theoretical physics background who will be very much at home with the language and techniques which are used. The author has developed a formalism for understanding complex systems which is based on the Langevin approach to the study of Brownian motion. This is a mesoscopic description; details of the interactions between the Brownian particle and the molecules of the surrounding fluid are replaced by a randomly fluctuating force. Thus all microscopic detail is replaced by a coarse-grained description which encapsulates the essence of the interactions at the finer level of description. In a similar way, the influences on Brownian agents in a multi-agent system are replaced by stochastic influences which sum up the effects of these interactions on a finer scale. Unlike Brownian particles, Brownian agents are not structureless particles, but instead have some internal states so that, for instance, they may react to changes in the environment or to the presence of other agents. Most of the book is concerned with developing the idea of Brownian agents using the techniques of statistical physics. This development parallels that for Brownian particles in physics, but the author then goes on to apply the technique to problems in biology, economics and the social sciences. This is a clear and well-written book which is a useful addition to the literature on complex systems. It will be interesting to see if the use of Brownian agents becomes
Long-time diffusivity of DNA chains in nanochannels: A Brownian dynamics study
Jain, Aashish; Dorfman, Kevin
2015-03-01
The simplest approach to calculate the diffusivity of any polymer chain is to use the double sum Kirkwood formula, which is based on preaveraging approximation of diffusion tensor. The error due to the preaveraging approximation has been reported by a number of researchers in the context of free solution by computing both Kirkwood diffusivity D (K) (also known as short-time diffusivity) and long-time diffusivity DL. In nanochannels, the main approach to compute the diffusivity is the Kirkwood formula. However, the error due to the preaveraging approximation is not known in a confined system. We use Brownian dynamics simulation algorithm with excluded volume and hydrodynamic interactions to calculate both short-time and long-time diffusivities of DNA chains in nanochannels, and compare them for a range of channel sizes and DNA chain sizes. Our results indicate that the long-time diffusivity is always smaller than the short-time diffusivity, which is consistent with the result obtained in free solution using linear response theory DL
Hybrid finite element and Brownian dynamics method for diffusion-controlled reactions
Bauler, Patricia; Huber, Gary A.; McCammon, J. Andrew
2012-01-01
Diffusion is often the rate determining step in many biological processes. Currently, the two main computational methods for studying diffusion are stochastic methods, such as Brownian dynamics, and continuum methods, such as the finite element method. This paper proposes a new hybrid diffusion method that couples the strengths of each of these two methods. The method is derived for a general multidimensional system, and is presented using a basic test case for 1D linear and radially symmetri...
Application of GPU processing for Brownian particle simulation
Cheng, Way Lee; Sheharyar, Ali; Sadr, Reza; Bouhali, Othmane
2015-01-01
Reports on the anomalous thermal-fluid properties of nanofluids (dilute suspension of nano-particles in a base fluid) have been the subject of attention for 15 years. The underlying physics that govern nanofluid behavior, however, is not fully understood and is a subject of much dispute. The interactions between the suspended particles and the base fluid have been cited as a major contributor to the improvement in heat transfer reported in the literature. Numerical simulations are instrumental in studying the behavior of nanofluids. However, such simulations can be computationally intensive due to the small dimensions and complexity of these problems. In this study, a simplified computational approach for isothermal nanofluid simulations was applied, and simulations were conducted using both conventional CPU and parallel GPU implementations. The GPU implementations significantly improved the computational performance, in terms of the simulation time, by a factor of 1000-2500. The results of this investigation show that, as the computational load increases, the simulation efficiency approaches a constant. At a very high computational load, the amount of improvement may even decrease due to limited system memory.
Differential Dynamic Microscopy to characterize Brownian motion and bacteria motility
Germain, David; Leocmach, Mathieu; Gibaud, Thomas
2015-01-01
We have developed a lab work module where we teach undergraduate students how to quantify the dynamics of a suspension of microscopic particles, measuring and analyzing the motion of those particles at the individual level or as a group. Differential Dynamic Microscopy (DDM) is a relatively recent technique that precisely does that and constitutes an alternative method to more classical techniques such as dynamics light scattering (DLS) or video particle tracking (VPT). DDM consists in imagin...
Speck, Thomas [Institut für Physik, Johannes Gutenberg-Universität Mainz, Staudingerweg 7-9, 55128 Mainz (Germany); Menzel, Andreas M.; Bialké, Julian; Löwen, Hartmut [Institut für Theoretische Physik II, Heinrich-Heine-Universität, D-40225 Düsseldorf (Germany)
2015-06-14
Recently, we have derived an effective Cahn-Hilliard equation for the phase separation dynamics of active Brownian particles by performing a weakly non-linear analysis of the effective hydrodynamic equations for density and polarization [Speck et al., Phys. Rev. Lett. 112, 218304 (2014)]. Here, we develop and explore this strategy in more detail and show explicitly how to get to such a large-scale, mean-field description starting from the microscopic dynamics. The effective free energy emerging from this approach has the form of a conventional Ginzburg-Landau function. On the coarsest scale, our results thus agree with the mapping of active phase separation onto that of passive fluids with attractive interactions through a global effective free energy (motility-induced phase transition). Particular attention is paid to the square-gradient term necessary for the phase separation kinetics. We finally discuss results from numerical simulations corroborating the analytical results.
Generalized Langevin Theory Of The Brownian Motion And The Dynamics Of Polymers In Solution
The review deals with a generalization of the Rouse and Zimm bead-spring models of the dynamics of flexible polymers in dilute solutions. As distinct from these popular theories, the memory in the polymer motion is taken into account. The memory naturally arises as a consequence of the fluid and bead inertia within the linearized Navier-Stokes hydrodynamics. We begin with a generalization of the classical theory of the Brownian motion, which forms the basis of any theory of the polymer dynamics. The random force driving the Brownian particles is not the white one as in the Langevin theory, but “colored”, i.e., statistically correlated in time, and the friction force on the particles depends on the history of their motion. An efficient method of solving the resulting generalized Langevin equations is presented and applied to the solution of the equations of motion of polymer beads. The memory effects lead to several peculiarities in the time correlation functions used to describe the dynamics of polymer chains. So, the mean square displacement of the polymer coils contains algebraic long-time tails and at short times it is ballistic. It is shown how these features reveal in the experimentally observable quantities, such as the dynamic structure factors of the scattering or the viscosity of polymer solutions. A phenomenological theory is also presented that describes the dependence of these quantities on the polymer concentration in solution. (author)
A Brownian motion technique to simulate gasification and its application to C/C composite ablation
Ablation of carbon-carbon composites (C/C) results in a heterogeneous surface recession mainly due to some gasification processes (oxidation, sublimation) possibly coupled to bulk mass transfer. In order to simulate and analyse the material/environment interactions during ablation, a Brownian motion simulation method featuring special Random Walk rules close to the wall has been implemented to efficiently simulate mass transfer in the low Peclet number regime. A sticking probability law adapted to this kind of Random Walk has been obtained for first-order heterogeneous reactions. In order to simulate the onset of surface roughness, the interface recession is simultaneously handled in 3D using a Simplified Marching Cube discretization. This tool is validated by comparison to analytical models. Then, its ability to provide reliable and accurate solutions of ablation phenomena in 3D is illustrated. (authors)
Generalized Scaling and the Master Variable for Brownian Magnetic Nanoparticle Dynamics
Reeves, Daniel B.; Shi, Yipeng; Weaver, John B.
2016-01-01
Understanding the dynamics of magnetic particles can help to advance several biomedical nanotechnologies. Previously, scaling relationships have been used in magnetic spectroscopy of nanoparticle Brownian motion (MSB) to measure biologically relevant properties (e.g., temperature, viscosity, bound state) surrounding nanoparticles in vivo. Those scaling relationships can be generalized with the introduction of a master variable found from non-dimensionalizing the dynamical Langevin equation. The variable encapsulates the dynamical variables of the surroundings and additionally includes the particles’ size distribution and moment and the applied field’s amplitude and frequency. From an applied perspective, the master variable allows tuning to an optimal MSB biosensing sensitivity range by manipulating both frequency and field amplitude. Calculation of magnetization harmonics in an oscillating applied field is also possible with an approximate closed-form solution in terms of the master variable and a single free parameter. PMID:26959493
Chernov, N.; Dolgopyat, D.
2008-01-01
A classical model of Brownian motion consists of a heavy molecule submerged into a gas of light atoms in a closed container. In this work we study a 2D version of this model, where the molecule is a heavy disk of mass M and the gas is represented by just one point particle of mass m = 1, which interacts with the disk and the walls of the container via elastic collisions. Chaotic behavior of the particles is ensured by convex (scattering) walls of the container. We prove that the position and ...
Shit, Anindita; Chattopadhyay, Sudip; Chaudhuri, Jyotipratim Ray
2011-06-01
We explore the Brownian dynamics in the quantum regime (by investigating the quantum Langevin and Smoluchowski equations) in terms of an effective time-independent Hamiltonian in the presence of a rapidly oscillating field. We achieve this by systematically expanding the time-dependent system-reservoir Hamiltonian in the inverse of driving frequency with a systematic time-scale separation and invoking a quantum gauge transformation within the framework of Floquet theorem. PMID:21797289
Elamin, Khalid; Swenson, Jan
2015-03-01
Aqueous solutions of glycerol are investigated by dynamic light scattering (DLS) over the whole concentration range (10-98 wt.% water) and in the temperature range 283-303 K. The measurements reveal one slow relaxation process in the geometry of polarized light scattering. This process is present in the whole concentration range, although it is very weak at the highest and lowest water concentrations and is considerably slower than the structural α relaxation, which is too fast to be observed on the experimental time scale in the measured temperature range. The relaxation time of the observed process exhibits a 1/q2 dependence, proving that it is due to long-range translational diffusion. The Stokes-Einstein relation is used to estimate the hydrodynamic radius of the diffusing particles and from these calculations it is evident that the observed relaxation process is due to the Brownian motion of single or a few glycerol molecules. The fact that it is possible to study the self-diffusion of such small molecules may stimulate a broadening of the research field used to be covered by the DLS technique. PMID:25871109
Chavanis, Pierre-Henri; Sire, Clément
2006-06-01
We propose a general kinetic and hydrodynamic description of self-gravitating Brownian particles in d dimensions. We go beyond the usual approximations by considering inertial effects and finite-N effects while previous works use a mean-field approximation valid in a proper thermodynamic limit (N --> +infinity) and consider an overdamped regime (xi --> +infinity). We recover known models in some particular cases of our general description. We derive the expression of the virial theorem for self-gravitating Brownian particles and study the linear dynamical stability of isolated clusters of particles and uniform systems by using techniques introduced in astrophysics. We investigate the influence of the equation of state, of the dimension of space, and of the friction coefficient on the dynamical stability of the system. We obtain the exact expression of the critical temperature Tc for a multicomponents self-gravitating Brownian gas in d = 2. We also consider the limit of weak frictions, xi --> 0, and derive the orbit-averaged Kramers equation. PMID:16906911
Brownian Dynamics of a Suspension of Particles with Constrained Voronoi Cell Volumes
Singh, John P.
2015-06-23
© 2015 American Chemical Society. Solvent-free polymer-grafted nanoparticle fluids consist of inorganic core particles fluidized by polymers tethered to their surfaces. The attachment of the suspending fluid to the particle surface creates a strong penalty for local variations in the fluid volume surrounding the particles. As a model of such a suspension we perform Brownian dynamics of an equilibrium system consisting of hard spheres which experience a many-particle potential proportional to the variance of the Voronoi volumes surrounding each particle (E = α(V
De Biase, Pablo M; Markosyan, Suren; Noskov, Sergei
2015-02-01
The transport of ions and solutes by biological pores is central for cellular processes and has a variety of applications in modern biotechnology. The time scale involved in the polymer transport across a nanopore is beyond the accessibility of conventional MD simulations. Moreover, experimental studies lack sufficient resolution to provide details on the molecular underpinning of the transport mechanisms. BROMOC, the code presented herein, performs Brownian dynamics simulations, both serial and parallel, up to several milliseconds long. BROMOC can be used to model large biological systems. IMC-MACRO software allows for the development of effective potentials for solute-ion interactions based on radial distribution function from all-atom MD. BROMOC Suite also provides a versatile set of tools to do a wide variety of preprocessing and postsimulation analysis. We illustrate a potential application with ion and ssDNA transport in MspA nanopore. PMID:25503688
Millen, J; Deesuwan, T; Barker, P; Anders, J
2014-06-01
Einstein realized that the fluctuations of a Brownian particle can be used to ascertain the properties of its environment. A large number of experiments have since exploited the Brownian motion of colloidal particles for studies of dissipative processes, providing insight into soft matter physics and leading to applications from energy harvesting to medical imaging. Here, we use heated optically levitated nanospheres to investigate the non-equilibrium properties of the gas surrounding them. Analysing the sphere's Brownian motion allows us to determine the temperature of the centre-of-mass motion of the sphere, its surface temperature and the heated gas temperature in two spatial dimensions. We observe asymmetric heating of the sphere and gas, with temperatures reaching the melting point of the material. This method offers opportunities for accurate temperature measurements with spatial resolution on the nanoscale, and provides a means for testing non-equilibrium thermodynamics. PMID:24793558
Mezzasalma, Stefano A
2007-03-15
The theoretical basis of a recent theory of Brownian relativity for polymer solutions is deepened and reexamined. After the problem of relative diffusion in polymer solutions is addressed, its two postulates are formulated in all generality. The former builds a statistical equivalence between (uncorrelated) timelike and shapelike reference frames, that is, among dynamical trajectories of liquid molecules and static configurations of polymer chains. The latter defines the "diffusive horizon" as the invariant quantity to work with in the special version of the theory. Particularly, the concept of universality in polymer physics corresponds in Brownian relativity to that of covariance in the Einstein formulation. Here, a "universal" law consists of a privileged observation, performed from the laboratory rest frame and agreeing with any diffusive reference system. From the joint lack of covariance and simultaneity implied by the Brownian Lorentz-Poincaré transforms, a relative uncertainty arises, in a certain analogy with quantum mechanics. It is driven by the difference between local diffusion coefficients in the liquid solution. The same transformation class can be used to infer Fick's second law of diffusion, playing here the role of a gauge invariance preserving covariance of the spacetime increments. An overall, noteworthy conclusion emerging from this view concerns the statistics of (i) static macromolecular configurations and (ii) the motion of liquid molecules, which would be much more related than expected. PMID:17223124
Quantal Brownian Motion from RPA dynamics: The master and Fokker-Planck equations
From the purely quantal RPA description of the damped harmonic oscillator and of the corresponding Brownian Motion within the full space (phonon subspace plus reservoir), a master equation (as well as a Fokker-Planck equation) for the reduced density matrix (for the reduced Wigner function, respectively) within the phonon subspace is extracted. The RPA master equation agrees with the master equation derived by the time-dependent perturbative approaches which utilize Tamm-Dancoff Hilbert spaces and invoke the rotating wave approximation. Since the RPA yields a full, as well as a contracted description, it can account for both the kinetic and the unperturbed oscillator momenta. The RPA description of the quantal Brownian Motion contrasts with the descriptions provided by the time perturbative approaches whether they invoke or not the rotating wave approximation. The RPA description also contrasts with the phenomenological phase space quantization. (orig.)
Dynamics of 2D Stochastic non-Newtonian fluids driven by fractional Brownian motion
Li, Jin; Huang, Jianhua
2011-01-01
A 2D Stochastic incompressible non-Newtonian fluids driven by fractional Bronwnian motion with Hurst parameter $H \\in (1/2,1)$ is studied. The Wiener-type stochastic integrals are introduced for infinite-dimensional fractional Brownian motion. Four groups of assumptions, including the requirement of Nuclear operator or Hilbert-Schmidt operator, are discussed. The existence and regularity of stochastic convolution for the corresponding additive linear stochastic equation are obtained under eac...
Chavanis, Pierre-Henri; Sire, Clément
2006-06-01
We derive the virial theorem appropriate to the generalized Smoluchowski-Poisson (GSP) system describing self-gravitating Brownian particles in an overdamped limit. We extend previous works by considering the case of an unbounded domain and an arbitrary equation of state. We use the virial theorem to study the diffusion (evaporation) of an isothermal Brownian gas above the critical temperature Tc in dimension d = 2 and show how the effective diffusion coefficient and the Einstein relation are modified by self-gravity. We also study the collapse at T = Tc and show that the central density increases logarithmically with time instead of exponentially in a bounded domain. Finally, for d > 2, we show that the evaporation of the system is essentially a pure diffusion slightly slowed down by self-gravity. We also study the linear dynamical stability of stationary solutions of the GSP system representing isolated clusters of particles and investigate the influence of the equation of state and of the dimension of space on the dynamical stability of the system. PMID:16906910
Lu, Zhumin; Dou, Hua-Shu; Phan-Thien, Nhan; Yeo, Khoon Seng
2008-01-01
In this paper, the finite element method is combined with the Brownian Configuration Field (BFC) method to simulate the fibre suspension flow in axisymmetric contraction and expansion passages. In order to solve for the high stress at high concentration, the Discrete Adaptive Viscoelastic Stress Splitting (DAVSS) method is employed. For the axisymmetric contraction and expansion passages with different geometry ratios, the results obtained are compared to available constitutive models and experiments. The predicted vortex length for dilute suspensions agrees well with experimental data in literature. Our numerical results show clearly the effect on vortex enhancement with increase of the volume fractions and the aspect ratios. Effect of aspect ratio of fibres on the vortex length is also studied. It is found that for the lower expansion ratio flows the vortex dimension in the corner region is fairly independent of fibre concentration and aspect ratio of fibres while the said vortex dimension increases with th...
Ryo Kanada
Full Text Available Kinesin is a family of molecular motors that move unidirectionally along microtubules (MT using ATP hydrolysis free energy. In the family, the conventional two-headed kinesin was experimentally characterized to move unidirectionally through "walking" in a hand-over-hand fashion by coordinated motions of the two heads. Interestingly a single-headed kinesin, a truncated KIF1A, still can generate a biased Brownian movement along MT, as observed by in vitro single molecule experiments. Thus, KIF1A must use a different mechanism from the conventional kinesin to achieve the unidirectional motions. Based on the energy landscape view of proteins, for the first time, we conducted a set of molecular simulations of the truncated KIF1A movements over an ATP hydrolysis cycle and found a mechanism exhibiting and enhancing stochastic forward-biased movements in a similar way to those in experiments. First, simulating stand-alone KIF1A, we did not find any biased movements, while we found that KIF1A with a large friction cargo-analog attached to the C-terminus can generate clearly biased Brownian movements upon an ATP hydrolysis cycle. The linked cargo-analog enhanced the detachment of the KIF1A from MT. Once detached, diffusion of the KIF1A head was restricted around the large cargo which was located in front of the head at the time of detachment, thus generating a forward bias of the diffusion. The cargo plays the role of a diffusional anchor, or cane, in KIF1A "walking."
Dynamics simulations for engineering macromolecular interactions
Robinson-Mosher, Avi; Shinar, Tamar; Silver, Pamela A.; Way, Jeffrey
2013-06-01
The predictable engineering of well-behaved transcriptional circuits is a central goal of synthetic biology. The artificial attachment of promoters to transcription factor genes usually results in noisy or chaotic behaviors, and such systems are unlikely to be useful in practical applications. Natural transcriptional regulation relies extensively on protein-protein interactions to insure tightly controlled behavior, but such tight control has been elusive in engineered systems. To help engineer protein-protein interactions, we have developed a molecular dynamics simulation framework that simplifies features of proteins moving by constrained Brownian motion, with the goal of performing long simulations. The behavior of a simulated protein system is determined by summation of forces that include a Brownian force, a drag force, excluded volume constraints, relative position constraints, and binding constraints that relate to experimentally determined on-rates and off-rates for chosen protein elements in a system. Proteins are abstracted as spheres. Binding surfaces are defined radially within a protein. Peptide linkers are abstracted as small protein-like spheres with rigid connections. To address whether our framework could generate useful predictions, we simulated the behavior of an engineered fusion protein consisting of two 20 000 Da proteins attached by flexible glycine/serine-type linkers. The two protein elements remained closely associated, as if constrained by a random walk in three dimensions of the peptide linker, as opposed to showing a distribution of distances expected if movement were dominated by Brownian motion of the protein domains only. We also simulated the behavior of fluorescent proteins tethered by a linker of varying length, compared the predicted Förster resonance energy transfer with previous experimental observations, and obtained a good correspondence. Finally, we simulated the binding behavior of a fusion of two ligands that could
Stokesian Dynamic Simulations of Colloid Assembly at a Fluid Interface
Dani, Archit; Maldarelli, Charles
2015-11-01
The collective dynamics and self-assembly of colloids floating at a gas/liquid or a liquid/liquid interface is a balance between deterministic lateral interaction forces, e.g. capillary attraction and dipolar electrostatic repulsion if the particles are charged, viscous resistance to colloid motion along the surface and thermal fluctuations. As the colloid size decreases, thermal (Brownian) forces become important and can affect the self assembly into ordered patterns and crystal structures that are the starting point for materials applications. Stokesian dynamics simulations are presented to describe the lateral organization of particles along the surface in Brownian dominated regimes that includes (using a pairwise approximation) capillary attraction and the hydrodynamic interaction between particles (incorporating the effect of the particle immersion depth) and thermal fluctuations. Clustering, fractal growth and particle ordering are observed at critically large values of the Peclet numbers, while smaller values yield states in which particles remain uncorrelated in space and more widely separated.
Role of Brownian motion on the thermal conductivity enhancement of nanofluids
Gupta, Amit; Kumar, Ranganathan
2007-11-01
This study involves Brownian dynamics simulations of a real nanofluid system in which the interparticle potential is determined based on Debye length and surface interaction of the fluid and the solid. This paper shows that Brownian motion can increase the thermal conductivity of the nanofluid by 6% primarily due to "random walk" motion and not only through diffusion. This increase is limited by the maximum concentration for each particle size and is below that predicted by the effective medium theory. Beyond the maximum limit, particle aggregates begin to form. Brownian motion contribution stays as a constant beyond a certain particle diameter.
Gomez-Marin, A.; Sancho, J. M.
2004-01-01
In this paper we present a model of a symmetric Brownian motor (SBM) which changes the sign of its velocity when the temperature gradient is inverted. The velocity, external work and efficiency are studied as a function of the temperatures of the baths and other relevant parameters. The motor shows a current reversal when another parameter (a phase shift) is varied. Analytical predictions and results from numerical simulations are performed and agree very well. Generic properties of this type...
Jian-Zhong Lin; Fu-Jun Gan
2012-01-01
The Brownian coagulation of nanoparticles with initial bimodal size distribution,i.e.,mode i and j,is numerically studied using the moment method.Evolutions of particle number concentration,geometric average diameter and geometric standard deviation are given in the free molecular regime,the continuum regime,the free molecular regimand transition regime,the free molecular regime and contin-uum regime,respectively.The results show that,both in the free molecular regime and the continuum regime,the number concentration of mode i and j decreases with increasing time.The evolutions of particle geometric average diameter with different initial size distribution are quite different.Both intra-modal and inter-modal coagulation finally make the polydispersed size distribution become monodispersed.As time goes by,the size distribution with initial bimodal turns to be unimodal and shifts to a larger particle size range.In the free molecular regime and transition regime,the intermodal coagulation becomes dominant when the number concentrations of mode i and j are of the same order.The effects of the number concentration of mode i and mode j on the evolution of geometric average diameter of mode j are negligible,while the effects of the number concentration of mode j on the evolution of geometric average diameter of mode j is distinct.In the free molecular regime and continuum regime,the higher the initial number concentration of mode j,the more obvious the variation of the number concentration of mode i.
Effect of interfaces on the nearby Brownian motion
Huang, Kai
2016-01-01
Near-boundary Brownian motion is a classic hydrodynamic problem of great importance in a variety of fields, from biophysics to micro-/nanofluidics. However, due to challenges in experimental measurements of near-boundary dynamics, the effect of interfaces on Brownian motion has remained elusive. Here, we report a computational study of this effect using microsecond-long large-scale molecular dynamics simulations and our newly developed Green-Kubo relation for friction at the liquid-solid interface. Our computer experiment unambiguously reveals that the t^(-3/2) long-time decay of the velocity autocorrelation function of a Brownian particle in bulk liquid is replaced by a t^(-5/2) decay near a boundary. We discover a general breakdown of traditional no-slip boundary condition at short time scales and we show that this breakdown has a profound impact on the near-boundary Brownian motion. Our results demonstrate the potential of Brownian-particle based micro-/nano-sonar to probe the local wettability of liquid-s...
Molecular dynamics simulations
Tarmyshov, Konstantin B.
2007-01-01
Molecular simulations can provide a detailed picture of a desired chemical, physical, or biological process. It has been developed over last 50 years and is being used now to solve a large variety of problems in many different fields. In particular, quantum calculations are very helpful to study small systems at a high resolution where electronic structure of compounds is accounted for. Molecular dynamics simulations, in turn, are employed to study development of a certain molecular ensemble ...
Interactive Dynamic-System Simulation
Korn, Granino A
2010-01-01
Showing you how to use personal computers for modeling and simulation, Interactive Dynamic-System Simulation, Second Edition provides a practical tutorial on interactive dynamic-system modeling and simulation. It discusses how to effectively simulate dynamical systems, such as aerospace vehicles, power plants, chemical processes, control systems, and physiological systems. Written by a pioneer in simulation, the book introduces dynamic-system models and explains how software for solving differential equations works. After demonstrating real simulation programs with simple examples, the author
Chavanis, Pierre-Henri; Sire, Clement
2005-01-01
We derive the Virial theorem appropriate to the generalized Smoluchowski-Poisson system describing self-gravitating Brownian particles and bacterial populations (chemotaxis). We extend previous works by considering the case of an unbounded domain and an arbitrary equation of state. We use the Virial theorem to study the diffusion (evaporation) of an isothermal Brownian gas above the critical temperature T_c in dimension d=2 and show how the effective diffusion coefficient and the Einstein rel...
Interacting Brownian Swarms: Some Analytical Results
Guillaume Sartoretti
2016-01-01
Full Text Available We consider the dynamics of swarms of scalar Brownian agents subject to local imitation mechanisms implemented using mutual rank-based interactions. For appropriate values of the underlying control parameters, the swarm propagates tightly and the distances separating successive agents are iid exponential random variables. Implicitly, the implementation of rank-based mutual interactions, requires that agents have infinite interaction ranges. Using the probabilistic size of the swarm’s support, we analytically estimate the critical interaction range below that flocked swarms cannot survive. In the second part of the paper, we consider the interactions between two flocked swarms of Brownian agents with finite interaction ranges. Both swarms travel with different barycentric velocities, and agents from both swarms indifferently interact with each other. For appropriate initial configurations, both swarms eventually collide (i.e., all agents interact. Depending on the values of the control parameters, one of the following patterns emerges after collision: (i Both swarms remain essentially flocked, or (ii the swarms become ultimately quasi-free and recover their nominal barycentric speeds. We derive a set of analytical flocking conditions based on the generalized rank-based Brownian motion. An extensive set of numerical simulations corroborates our analytical findings.
Wu, Xiaobin; Qiu, Jian; Luo, Kaiqing; Han, Peng
2015-08-01
Dynamic Light Scattering is used for measuring particle size distribution of nano-particle under Brownian motion. Signal is detected through a photomultiplier and processed by correlation analysis, and results are inverted at last. Method by using CCD camera can record the procedure of motion. However, there are several weaknesses such as low refresh speed and noise from CCD camera, and this method depends on particle size and detecting angle. A simulation of nano-particle under Brownian motion is proposed to record dynamic images, studies contrast of dynamic images which can represent speed of diffusion, and its characteristic under different conditions. The results show that through contrast of dynamic images diffusion coefficient can be obtained, which is independent on density of scattering volume.
Molecular dynamics simulations
The molecular dynamics computer simulation discovery of the slow decay of the velocity autocorrelation function in fluids is briefly reviewed in order to contrast that long time tail with those observed for the stress autocorrelation function in fluids and the velocity autocorrelation function in the Lorentz gas. For a non-localized particle in the Lorentz gas it is made plausible that even if it behaved quantum mechanically its long time tail would be the same as the classical one. The generalization of Fick's law for diffusion for the Lorentz gas, necessary to avoid divergences due to the slow decay of correlations, is presented. For fluids, that generalization has not yet been established, but the region of validity of generalized hydrodynamics is discussed. 20 refs., 5 figs
Blending Brownian motion and heat equation
Cristiani, Emiliano
2015-01-01
In this short communication we present an original way to couple the Brownian motion and the heat equation. More in general, we suggest a way for coupling the Langevin equation for a particle, which describes a single realization of its trajectory, with the associated Fokker-Planck equation, which instead describes the evolution of the particle's probability density function. Numerical results show that it is indeed possible to obtain a regularized Brownian motion and a Brownianized heat equation still preserving the global statistical properties of the solutions. The results also suggest that the more macroscale leads the dynamics the more one can reduce the microscopic degrees of freedom.
Canonical active Brownian motion
Gluck, Alexander; Huffel, Helmuth; Ilijic, Sasa
2008-01-01
Active Brownian motion is the complex motion of active Brownian particles. They are active in the sense that they can transform their internal energy into energy of motion and thus create complex motion patterns. Theories of active Brownian motion so far imposed couplings between the internal energy and the kinetic energy of the system. We investigate how this idea can be naturally taken further to include also couplings to the potential energy, which finally leads to a general theory of cano...
A multiscale guide to Brownian motion
Grebenkov, Denis S.; Belyaev, Dmitry; Jones, Peter W.
2016-01-01
We revise the Lévy construction of Brownian motion as a simple though rigorous approach to operate with various Gaussian processes. A Brownian path is explicitly constructed as a linear combination of wavelet-based ‘geometrical features’ at multiple length scales with random weights. Such a wavelet representation gives a closed formula mapping of the unit interval onto the functional space of Brownian paths. This formula elucidates many classical results about Brownian motion (e.g., non-differentiability of its path), providing an intuitive feeling for non-mathematicians. The illustrative character of the wavelet representation, along with the simple structure of the underlying probability space, is different from the usual presentation of most classical textbooks. Similar concepts are discussed for the Brownian bridge, fractional Brownian motion, the Ornstein-Uhlenbeck process, Gaussian free fields, and fractional Gaussian fields. Wavelet representations and dyadic decompositions form the basis of many highly efficient numerical methods to simulate Gaussian processes and fields, including Brownian motion and other diffusive processes in confining domains.
A dynamic mesh refinement technique for Lattice Boltzmann simulations on octree-like grids
Neumann, Philipp
2012-04-27
In this contribution, we present our new adaptive Lattice Boltzmann implementation within the Peano framework, with special focus on nanoscale particle transport problems. With the continuum hypothesis not holding anymore on these small scales, new physical effects - such as Brownian fluctuations - need to be incorporated. We explain the overall layout of the application, including memory layout and access, and shortly review the adaptive algorithm. The scheme is validated by different benchmark computations in two and three dimensions. An extension to dynamically changing grids and a spatially adaptive approach to fluctuating hydrodynamics, allowing for the thermalisation of the fluid in particular regions of interest, is proposed. Both dynamic adaptivity and adaptive fluctuating hydrodynamics are validated separately in simulations of particle transport problems. The application of this scheme to an oscillating particle in a nanopore illustrates the importance of Brownian fluctuations in such setups. © 2012 Springer-Verlag.
Cooperative Transport of Brownian Particles
Derenyi, Imre; Vicsek, Tamas
1998-01-01
We consider the collective motion of finite-sized, overdamped Brownian particles (e.g., motor proteins) in a periodic potential. Simulations of our model have revealed a number of novel cooperative transport phenomena, including (i) the reversal of direction of the net current as the particle density is increased and (ii) a very strong and complex dependence of the average velocity on both the size and the average distance of the particles.
Popov, Ivan; Vitkin, Alex
2016-01-01
The study of flowing Brownian particles finds numerous biomedical applications, ranging from blood flow analysis to diffusion research. A mathematical model for the correlation function of laser radiation scattered by flowing Brownian particles measured with fiber-based optical coherence tomography (OCT), which accounts for the effects of sample arm optics, is presented. It is shown that the parameters of an OCT optical system of any complexity can be taken into account by using the ABCD ray tracing matrix approach. Specifically, the impact of any optical system can be characterized by the changes in the effective beam radius, which replaces the Gaussian beam radius in the existing mathematical models of scattered radiation. It is shown that the validity of the developed ABCD matrix formalism is governed by the condition that the source coherence length is much smaller than the Rayleigh range in the sample. The predictions of the developed model are compared with previously published theories and with experimental data and agree well with the latter.
Lawler, Gregory F.; Werner, Wendelin
2003-01-01
We define a natural conformally invariant measure on unrooted Brownian loops in the plane and study some of its properties. We relate this measure to a measure on loops rooted at a boundary point of a domain and show how this relation gives a way to ``chronologically add Brownian loops'' to simple curves in the plane.
Radiation Reaction on Brownian Motions
Seto, Keita
2016-01-01
Tracking the real trajectory of a quantum particle is one of the interpretation problem and it is expressed by the Brownian (stochastic) motion suggested by E. Nelson. Especially the dynamics of a radiating electron, namely, radiation reaction which requires us to track its trajectory becomes important in the high-intensity physics by PW-class lasers at present. It has been normally treated by the Furry picture in non-linear QED, but it is difficult to draw the real trajectory of a quantum particle. For the improvement of this, I propose the representation of a stochastic particle interacting with fields and show the way to describe radiation reaction on its Brownian motion.
Shen, Meng; Lueptow, Richard M
2016-01-01
The Angstrom-scale transport characteristics of water and six different solutes, methanol, ethanol, 2-propanol, urea, Na+, and Cl-, were studied for a polyamide reverse osmosis (RO) membrane, FT-30, using non-equilibrium molecular dynamics (NEMD) simulations. Results indicate that water transport increases with an increasing fraction of connected percolated free volume, or water-accessible open space, in the membrane polymer structure. This free volume is enhanced by the dynamic structure of the membrane at the molecular level as it swells when hydrated and vibrates due to molecular collisions allowing a continuous path connecting the opposite membrane surfaces. The tortuous paths available for transport of solutes result in Brownian motion of solute molecules and hopping from pore to pore as they pass through the polymer network structure of the membrane. The transport of alcohol solutes decreases for solutes with larger Van der Waals volume, which corresponds to less available percolated free volume, or sol...
Classical molecular dynamics simulation on the dynamical properties of H2 on silicene layer
Casuyac Miqueas
2016-01-01
Full Text Available This study investigates the diffusion of hydrogen molecule physisorbed on the surface of silicene nanoribbon (SiNRusing the classical molecular dynamic (MD simulation in LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator. The interactions between silicon atoms are modeled using the modified Tersoff potential, the Adaptive Intermolecular Reactive Empirical Bond Order (AIREBO potential for hydrogen – hydrogen interaction and the Lennard – Jones potential for the physisorbed H2 on SiNR. By varying the temperatures (60 K Δ 130 K, we observed that the Δxdisplacement of H2 on the surface SiNR shows a Brownian motion on a Lennard-Jones potential and a Gaussian probability distribution can be plotted describing the diffusion of H2. The calculated mean square displacement (MSD was approximately increasing in time and the activation energy barrier for diffusion has been found to be 43.23meV.
Molecular dynamics simulation of diffusivity
Juanfang LIU; Danling ZENG; Qin LI; Hong GAO
2008-01-01
Equilibrium molecular dynamics simulation was performed on water to calculate its diffusivity by adopting different potential models. The results show that the potential models have great influence on the simulated results. In addition, the diffusivities obtained by the SPCE model conform well to the experimental values.
Quantum trajectories for Brownian motion
Strunz, W T; Gisin, Nicolas; Yu, T; Strunz, Walter T.; Diosi, Lajos; Gisin, Nicolas
1999-01-01
We present the stochastic Schroedinger equation for the dynamics of a quantum particle coupled to a high temperature environment and apply it the dynamics of a driven, damped, nonlinear quantum oscillator. Apart from an initial slip on the environmental memory time scale, in the mean, our result recovers the solution of the known non-Lindblad quantum Brownian motion master equation. A remarkable feature of our approach is its localization property: individual quantum trajectories remain localized wave packets for all times, even for the classically chaotic system considered here, the localization being stronger the smaller $\\hbar$.
Brownian Motion in Minkowski Space
Paul O'Hara
2015-06-01
Full Text Available We construct a model of Brownian motion in Minkowski space. There are two aspects of the problem. The first is to define a sequence of stopping times associated with the Brownian “kicks” or impulses. The second is to define the dynamics of the particle along geodesics in between the Brownian kicks. When these two aspects are taken together, the Central Limit Theorem (CLT leads to temperature dependent four dimensional distributions defined on Minkowski space, for distances and 4-velocities. In particular, our processes are characterized by two independent time variables defined with respect to the laboratory frame: a discrete one corresponding to the stopping times when the impulses take place and a continuous one corresponding to the geodesic motion in-between impulses. The subsequent distributions are solutions of a (covariant pseudo-diffusion equation which involves derivatives with respect to both time variables, rather than solutions of the telegraph equation which has a single time variable. This approach simplifies some of the known problems in this context.
Role of Brownian Motion Hydrodynamics on Nanofluid Thermal Conductivity
W Evans, J Fish, P Keblinski
2005-11-14
We use a simple kinetic theory based analysis of heat flow in fluid suspensions of solid nanoparticles (nanofluids) to demonstrate that the hydrodynamics effects associated with Brownian motion have a minor effect on the thermal conductivity of the nanofluid. Our conjecture is supported by the results of molecular dynamics simulations of heat flow in a model nanofluid with well-dispersed particles. Our findings are consistent with the predictions of the effective medium theory as well as with recent experimental results on well dispersed metal nanoparticle suspensions.
Coupling of lever arm swing and biased Brownian motion in actomyosin.
Qing-Miao Nie
2014-04-01
Full Text Available An important unresolved problem associated with actomyosin motors is the role of Brownian motion in the process of force generation. On the basis of structural observations of myosins and actins, the widely held lever-arm hypothesis has been proposed, in which proteins are assumed to show sequential structural changes among observed and hypothesized structures to exert mechanical force. An alternative hypothesis, the Brownian motion hypothesis, has been supported by single-molecule experiments and emphasizes more on the roles of fluctuating protein movement. In this study, we address the long-standing controversy between the lever-arm hypothesis and the Brownian motion hypothesis through in silico observations of an actomyosin system. We study a system composed of myosin II and actin filament by calculating free-energy landscapes of actin-myosin interactions using the molecular dynamics method and by simulating transitions among dynamically changing free-energy landscapes using the Monte Carlo method. The results obtained by this combined multi-scale calculation show that myosin with inorganic phosphate (Pi and ADP weakly binds to actin and that after releasing Pi and ADP, myosin moves along the actin filament toward the strong-binding site by exhibiting the biased Brownian motion, a behavior consistent with the observed single-molecular behavior of myosin. Conformational flexibility of loops at the actin-interface of myosin and the N-terminus of actin subunit is necessary for the distinct bias in the Brownian motion. Both the 5.5-11 nm displacement due to the biased Brownian motion and the 3-5 nm displacement due to lever-arm swing contribute to the net displacement of myosin. The calculated results further suggest that the recovery stroke of the lever arm plays an important role in enhancing the displacement of myosin through multiple cycles of ATP hydrolysis, suggesting a unified movement mechanism for various members of the myosin family.
Simulations with dynamical HISQ quarks
Bazavov, A; DeTar, C; Freeman, W; Gottlieb, Steven; Heller, U M; Hetrick, J E; Laiho, J; Levkova, L; Oktay, M; Osborn, J; Sugar, R L; Toussaint, D; Van de Water, R S
2010-01-01
We report on the status of a program of generating and using configurations with four flavors of dynamical quarks, using the HISQ action. We study the lattice spacing dependence of physical quantities in these simulations, using runs at several lattice spacings, but with the light quark mass held fixed at two tenths of the strange quark mass. We find that the lattice artifacts in the HISQ simulations are much smaller than those in the asqtad simulations at the same lattice spacings and quark masses. We also discuss methods for setting the scale, or assigning a lattice spacing to ensembles run at unphysical parameters.
The key to obtaining the model-free description of the dynamics of a macromolecule is the optimisation of the model-free and Brownian rotational diffusion parameters using the collected R1, R2 and steady-state NOE relaxation data. The problem of optimising the chi-squared value is often assumed to be trivial, however, the long chain of dependencies required for its calculation complicates the model-free chi-squared space. Convolutions are induced by the Lorentzian form of the spectral density functions, the linear recombinations of certain spectral density values to obtain the relaxation rates, the calculation of the NOE using the ratio of two of these rates, and finally the quadratic form of the chi-squared equation itself. Two major topological features of the model-free space complicate optimisation. The first is a long, shallow valley which commences at infinite correlation times and gradually approaches the minimum. The most severe convolution occurs for motions on two timescales in which the minimum is often located at the end of a long, deep, curved tunnel or multidimensional valley through the space. A large number of optimisation algorithms will be investigated and their performance compared to determine which techniques are suitable for use in model-free analysis. Local optimisation algorithms will be shown to be sufficient for minimisation not only within the model-free space but also for the minimisation of the Brownian rotational diffusion tensor. In addition the performance of the programs Modelfree and Dasha are investigated. A number of model-free optimisation failures were identified: the inability to slide along the limits, the singular matrix failure of the Levenberg-Marquardt minimisation algorithm, the low precision of both programs, and a bug in Modelfree. Significantly, the singular matrix failure of the Levenberg-Marquardt algorithm occurs when internal correlation times are undefined and is greatly amplified in model-free analysis by both
Brownian particles in supramolecular polymer solutions
Gucht, van der J.; Besseling, N.A.M.; Knoben, W.; Bouteiller, L.; Cohen Stuart, M.A.
2003-01-01
The Brownian motion of colloidal particles embedded in solutions of hydrogen-bonded supramolecular polymers has been studied using dynamic light scattering. At short times, the motion of the probe particles is diffusive with a diffusion coefficient equal to that in pure solvent. At intermediate time
An efficient dissipative particle dynamics-based algorithm for simulating electrolyte solutions
Medina, Stefan; Wang, Zhen-Gang; Schmid, Friederike
2014-01-01
We propose an efficient simulation algorithm based on the dissipative particle dynamics (DPD) method for studying electrohydrodynamic phenomena in electrolyte fluids. The fluid flow is mimicked with DPD particles while the evolution of the concentration of the ionic species is described using Brownian pseudo particles. The method is designed especially for systems with high salt concentrations, as explicit treatment of the salt ions becomes computationally expensive. For illustration, we apply the method to electro-osmotic flow over patterned, superhydrophobic surfaces. The results are in good agreement with recent theoretical predictions.
Vehicle dynamics modeling and simulation
Schramm, Dieter; Bardini, Roberto
2014-01-01
The authors examine in detail the fundamentals and mathematical descriptions of the dynamics of automobiles. In this context different levels of complexity will be presented, starting with basic single-track models up to complex three-dimensional multi-body models. A particular focus is on the process of establishing mathematical models on the basis of real cars and the validation of simulation results. The methods presented are explained in detail by means of selected application scenarios.
Molecular dynamics simulation of benzene
Trumpakaj, Zygmunt; Linde, Bogumił B. J.
2016-03-01
Intermolecular potentials and a few models of intermolecular interaction in liquid benzene are tested by Molecular Dynamics (MD) simulations. The repulsive part of the Lennard-Jones 12-6 (LJ 12-6) potential is too hard, which yields incorrect results. The exp-6 potential with a too hard repulsive term is also often used. Therefore, we took an expa-6 potential with a small Gaussian correction plus electrostatic interactions. This allows to modify the curvature of the potential. The MD simulations are carried out in the temperature range 280-352 K under normal pressure and at experimental density. The Rayleigh scattering of depolarized light is used for comparison. The results of MD simulations are comparable with the experimental values.
Dynamic simulations of tissue welding
Maitland, D.J.; Eder, D.C.; London, R.A.; Glinsky, M.E. [and others
1996-02-01
The exposure of human skin to near-infrared radiation is numerically simulated using coupled laser, thermal transport and mass transport numerical models. The computer model LATIS is applied in both one-dimensional and two-dimensional geometries. Zones within the skin model are comprised of a topical solder, epidermis, dermis, and fatty tissue. Each skin zone is assigned initial optical, thermal and water density properties consistent with values listed in the literature. The optical properties of each zone (i.e. scattering, absorption and anisotropy coefficients) are modeled as a kinetic function of the temperature. Finally, the water content in each zone is computed from water diffusion where water losses are accounted for by evaporative losses at the air-solder interface. The simulation results show that the inclusion of water transport and evaporative losses in the model are necessary to match experimental observations. Dynamic temperature and damage distributions are presented for the skin simulations.
Atzberger, P. J.
2007-01-01
In this paper a direct correspondence is made between the effective stochastic dynamics of elastic structures of an Immersed Boundary Method incorporating thermal fluctuations and Stokesian-Browman Dynamics. The correspondence is made in the limit of small Reynolds number, in which the fluid relaxes rapidly on the time scale of the motion of the immersed structures, by performing an averaging procedure directly on the stochastic equations of the Immersed Boundary Method. It is found that the...
Noncommutative Brownian motion
Santos, Willien O; Souza, Andre M C
2016-01-01
We investigate the Brownian motion of a particle in a two-dimensional noncommutative (NC) space. Using the standard NC algebra embodied by the sympletic Weyl-Moyal formalism we find that noncommutativity induces a non-vanishing correlation between both coordinates at different times. The effect itself stands as a signature of spatial noncommutativity and offers further alternatives to experimentally detect the phenomena.
Meurs, P.; Broeck, C. Van Den
2005-01-01
Recently, a thermal Brownian motor was introduced [Van den Broeck, Kawai and Meurs, Phys. Rev. Lett. (2004)], for which an exact microscopic analysis is possible. The purpose of this paper is to review some further properties of this construction, and to discuss in particular specific issues including the relation with macroscopic response and the efficiency at maximum power.
Finding the dynamics of an entire macromolecule is a complex problem as the model-free parameter values are intricately linked to the Brownian rotational diffusion of the molecule, mathematically through the autocorrelation function of the motion and statistically through model selection. The solution to this problem was formulated using set theory as an element of the universal set U-the union of all model-free spaces (d'Auvergne EJ and Gooley PR (2007) Mol BioSyst 3(7), 483-494). The current procedure commonly used to find the universal solution is to initially estimate the diffusion tensor parameters, to optimise the model-free parameters of numerous models, and then to choose the best model via model selection. The global model is then optimised and the procedure repeated until convergence. In this paper a new methodology is presented which takes a different approach to this diffusion seeded model-free paradigm. Rather than starting with the diffusion tensor this iterative protocol begins by optimising the model-free parameters in the absence of any global model parameters, selecting between all the model-free models, and finally optimising the diffusion tensor. The new model-free optimisation protocol will be validated using synthetic data from Schurr JM et al. (1994) J Magn Reson B 105(3), 211-224 and the relaxation data of the bacteriorhodopsin (1-36)BR fragment from Orekhov VY (1999) J Biomol NMR 14(4), 345-356. To demonstrate the importance of this new procedure the NMR relaxation data of the Olfactory Marker Protein (OMP) of Gitti R et al. (2005) Biochem 44(28), 9673-9679 is reanalysed. The result is that the dynamics for certain secondary structural elements is very different from those originally reported
Wenzheng Cui
2015-09-01
Full Text Available Nanofluids are a new generation of high-efficiency refrigerant with abnormal increased thermal conductivity and convective heat transfer properties. In view of the paucity of research work on the contribution of nanoparticle Brownian motion for the thermal conductivity augmentation, the present paper carries out a series of MD simulations to explorer the order of magnitude of nanoparticle Brownian motion and discusses the effect of nanoparticle Brownian motion for thermal conductivity enhancement of nanofluids. Various influence factors including nanoparticle shapes, sizes, and materials are considered. The Brownian motion of nanoparticles is decomposed into rotation and migration and calculated by MD simulation. By means of Peclet number, the effect of nanoparticle Brownian motion for thermal conductivity enhancement of nanofluids is discussed.
Humanoid Robot Simulator: A Realistic Dynamics Approach
Lima, José; Gonçalves, José; Costa, Paulo; Moreira, António
2008-01-01
This paper describes a humanoid robot simulator with realistic dynamics. As simulation is a powerful tool for speeding up the control software development, the suggested accurate simulator allows to accomplish this goal. The simulator, based on the Open Dynamics Engine and GLScene graphics library, provides instant visual feedback and allows the user to test any control strategy without damaging the real robot in the early stages of the development. The proposed simulator also captures some c...
Brownian Motion, "Diverse and Undulating"
Duplantier, Bertrand
2016-01-01
We describe in detail the history of Brownian motion, as well as the contributions of Einstein, Sutherland, Smoluchowski, Bachelier, Perrin and Langevin to its theory. The always topical importance in physics of the theory of Brownian motion is illustrated by recent biophysical experiments, where it serves, for instance, for the measurement of the pulling force on a single DNA molecule. In a second part, we stress the mathematical importance of the theory of Brownian motion, illustrated by two chosen examples. The by-now classic representation of the Newtonian potential by Brownian motion is explained in an elementary way. We conclude with the description of recent progress seen in the geometry of the planar Brownian curve. At its heart lie the concepts of conformal invariance and multifractality, associated with the potential theory of the Brownian curve itself.
Trefan, Gyorgy
1993-01-01
The goal of this thesis is to contribute to the ambitious program of the foundation of developing statistical physics using chaos. We build a deterministic model of Brownian motion and provide a microscopic derivation of the Fokker-Planck equation. Since the Brownian motion of a particle is the result of the competing processes of diffusion and dissipation, we create a model where both diffusion and dissipation originate from the same deterministic mechanism--the deterministic interaction of that particle with its environment. We show that standard diffusion which is the basis of the Fokker-Planck equation rests on the Central Limit Theorem, and, consequently, on the possibility of deriving it from a deterministic process with a quickly decaying correlation function. The sensitive dependence on initial conditions, one of the defining properties of chaos insures this rapid decay. We carefully address the problem of deriving dissipation from the interaction of a particle with a fully deterministic nonlinear bath, that we term the booster. We show that the solution of this problem essentially rests on the linear response of a booster to an external perturbation. This raises a long-standing problem concerned with Kubo's Linear Response Theory and the strong criticism against it by van Kampen. Kubo's theory is based on a perturbation treatment of the Liouville equation, which, in turn, is expected to be totally equivalent to a first-order perturbation treatment of single trajectories. Since the boosters are chaotic, and chaos is essential to generate diffusion, the single trajectories are highly unstable and do not respond linearly to weak external perturbation. We adopt chaotic maps as boosters of a Brownian particle, and therefore address the problem of the response of a chaotic booster to an external perturbation. We notice that a fully chaotic map is characterized by an invariant measure which is a continuous function of the control parameters of the map
Combinatorial fractal Brownian motion model
朱炬波; 梁甸农
2000-01-01
To solve the problem of how to determine the non-scaled interval when processing radar clutter using fractal Brownian motion (FBM) model, a concept of combinatorial FBM model is presented. Since the earth (or sea) surface varies diversely with space, a radar clutter contains several fractal structures, which coexist on all scales. Taking the combination of two FBMs into account, via theoretical derivation we establish a combinatorial FBM model and present a method to estimate its fractal parameters. The correctness of the model and the method is proved by simulation experiments and computation of practial data. Furthermore, we obtain the relationship between fractal parameters when processing combinatorial model with a single FBM model. Meanwhile, by theoretical analysis it is concluded that when combinatorial model is observed on different scales, one of the fractal structures is more obvious.
Simulation of Fault Zone Dynamics
Mora, P.; Abe, S.; Place, D.
2002-12-01
Particle models such as the discrete element model for granular assemblies and the lattice solid model provide a means to study the dynamics of fault zones. The lattice solid model was developed with the aim of progressively building up the capacity to simulate all relevent physical processes in fault zones. The present implementation of the model is able to simulate the dynamics of a granular lattice consisting of bonded or unbonded circular (2D) or spherical (3D) particles. Thermal effects (frictional hear generation, thermal expansion, heat flow) and pore fluid effects (heat induced pore pressure gradients and the consequent Darcian flow and impact on effective friction) can be modelled. Past work involving both circular particles and non-circular grains constructed as groups of bonded particles have demonstrated that grain shape has a fundamental impact on zero-th order behaviour. When circular particles are used, rolling is the most efficient means to accomodate slip of a simulated fault gouge layer leading to unrealistically low friction, typically around 0.2. This is consistent with laboratory results by Mair and Marone which have demonstrated that gouge consisting entirely of spherical beads shows a lower coefficient of friction than gouge containing irregular shaped particles. Recent work comparing quasi-2D laboratory results using pasta (Marone) with 2D numerical results (Morgan) have confirmed that numerical and laboratory results with circular ``particles'' are in agreement. When irregular grains are modelled at the lowest scale, the friction of simulated gouge layers matches with laboratory observations of rock friction (μ ~ 0.6) and is insentitive to the value used for interparticle friction (Mora et al, 2000). This indicates a self-regulation mechanism is occurring in which the group behaviour of the gouge layer remains constant at around 0.6 by balancing the amount of slip and rolling of grains within the gouge layer. A limitation of these studies
Brownian semi-stationary processes, turbulence and smooth processes
Urbina, José Ulises Márquez
This thesis analysis the use of Brownian semi-stationary (BSS) processes to model the main statistical features present in turbulent time series, and some asymptotic properties of certain classes of smooth processes. Turbulence is a complex phenomena governed by the Navier-Stokes equations. These...... equations do not represent a fully functional model and, consequently, it has been necessary to develop phenomenological models capturing main aspects of turbulent dynamics. The BSS processes were proposed as an option to model turbulent time series. In this thesis we proved, through a simulation....... We also studied the distributional properties of the increments of BSS processes with the intent to better understand why the BSS processes seem to accurately reproduce the temporal turbulent dynamics. BSS processes in general are not semimartingales. However, there are conditions which make a BSS...
QUANTUM STOCHASTIC PROCESSES: BOSON AND FERMION BROWNIAN MOTION
A.E.Kobryn
2003-01-01
Full Text Available Dynamics of quantum systems which are stochastically perturbed by linear coupling to the reservoir can be studied in terms of quantum stochastic differential equations (for example, quantum stochastic Liouville equation and quantum Langevin equation. In order to work it out one needs to define the quantum Brownian motion. As far as only its boson version has been known until recently, in the present paper we present the definition which makes it possible to consider the fermion Brownian motion as well.
Holographic Brownian Motion at Finite Density
Banerjee, Pinaki
2015-01-01
We study holographic Brownian motion of a heavy charged particle at zero and small (but finite) temperature in presence of finite density. We are primarily interested in the dynamics at (near) zero temperature which is holographically described by motion of a fundamental string in an (near-) extremal RN black hole. We compute analytically the functional form of retarded Green's function and also compare that numerically at leading order in small frequency.
Brownian particles in supramolecular polymer solutions
Gucht, van der, J.; Besseling, N.A.M.; Knoben, W.; Bouteiller, L; Cohen Stuart, M. A.
2003-01-01
The Brownian motion of colloidal particles embedded in solutions of hydrogen-bonded supramolecular polymers has been studied using dynamic light scattering. At short times, the motion of the probe particles is diffusive with a diffusion coefficient equal to that in pure solvent. At intermediate time scales the particles are slowed down as a result of trapping in elastic cages formed by the polymer chains, while at longer times the motion is diffusive again, but with a much smaller diffusion c...
Enabling Strain Hardening Simulations with Dislocation Dynamics
Arsenlis, A; Cai, W
2006-12-20
Numerical algorithms for discrete dislocation dynamics simulations are investigated for the purpose of enabling strain hardening simulations of single crystals on massively parallel computers. The algorithms investigated include the /(N) calculation of forces, the equations of motion, time integration, adaptive mesh refinement, the treatment of dislocation core reactions, and the dynamic distribution of work on parallel computers. A simulation integrating all of these algorithmic elements using the Parallel Dislocation Simulator (ParaDiS) code is performed to understand their behavior in concert, and evaluate the overall numerical performance of dislocation dynamics simulations and their ability to accumulate percents of plastic strain.
Frustrated Brownian Motion of Nonlocal Solitary Waves
We investigate the evolution of solitary waves in a nonlocal medium in the presence of disorder. By using a perturbational approach, we show that an increasing degree of nonlocality may largely hamper the Brownian motion of self-trapped wave packets. The result is valid for any kind of nonlocality and in the presence of nonparaxial effects. Analytical predictions are compared with numerical simulations based on stochastic partial differential equations.
The open quantum Brownian motions
Using quantum parallelism on random walks as the original seed, we introduce new quantum stochastic processes, the open quantum Brownian motions. They describe the behaviors of quantum walkers—with internal degrees of freedom which serve as random gyroscopes—interacting with a series of probes which serve as quantum coins. These processes may also be viewed as the scaling limit of open quantum random walks and we develop this approach along three different lines: the quantum trajectory, the quantum dynamical map and the quantum stochastic differential equation. We also present a study of the simplest case, with a two level system as an internal gyroscope, illustrating the interplay between the ballistic and diffusive behaviors at work in these processes. Notation Hz: orbital (walker) Hilbert space, CZ in the discrete, L2(R) in the continuum Hc: internal spin (or gyroscope) Hilbert space Hsys=Hz⊗Hc: system Hilbert space Hp: probe (or quantum coin) Hilbert space, Hp=C2 ρttot: density matrix for the total system (walker + internal spin + quantum coins) ρ-bar t: reduced density matrix on Hsys: ρ-bar t=∫dxdy ρ-bar t(x,y)⊗|x〉z〈y| ρ-hat t: system density matrix in a quantum trajectory: ρ-hat t=∫dxdy ρ-hat t(x,y)⊗|x〉z〈y|. If diagonal and localized in position: ρ-hat t=ρt⊗|Xt〉z〈Xt| ρt: internal density matrix in a simple quantum trajectory Xt: walker position in a simple quantum trajectory Bt: normalized Brownian motion ξt, ξt†: quantum noises (paper)
Molecular dynamics simulation by atomic mass weighting
Mao, Boryeu; Friedman, Alan R.
1990-01-01
A molecular dynamics-based simulation method in which atomic masses are weighted is described. Results from this method showed that the capability for conformation search in molecular dynamics simulation of a short peptide (FMRF-amide) is significantly increased by mass weighting.
Satoshi Ota
2016-09-01
Full Text Available The dependence of magnetic relaxation on particle parameters, such as the size and anisotropy, has been conventionally discussed. In addition, the influences of external conditions, such as the intensity and frequency of the applied field, the surrounding viscosity, and the temperature on the magnetic relaxation have been researched. According to one of the basic theories regarding magnetic relaxation, the faster type of relaxation dominates the process. However, in this study, we reveal that Brownian and Néel relaxations coexist and that Brownian relaxation can occur after Néel relaxation despite having a longer relaxation time. To understand the mechanisms of Brownian rotation, alternating current (AC hysteresis loops were measured in magnetic fluids of different viscosities. These loops conveyed the amplitude and phase delay of the magnetization. In addition, the intrinsic loss power (ILP was calculated using the area of the AC hysteresis loops. The ILP also showed the magnetization response regarding the magnetic relaxation over a wide frequency range. To develop biomedical applications of magnetic nanoparticles, such as hyperthermia and magnetic particle imaging, it is necessary to understand the mechanisms of magnetic relaxation.
Accelerated molecular dynamics simulations of protein folding
Miao, Y.; Feixas, F; Eun, C; McCammon, JA
2015-01-01
© 2015 Wiley Periodicals, Inc. Folding of four fast-folding proteins, including chignolin, Trp-cage, villin headpiece and WW domain, was simulated via accelerated molecular dynamics (aMD). In comparison with hundred-of-microsecond timescale conventional molecular dynamics (cMD) simulations performed on the Anton supercomputer, aMD captured complete folding of the four proteins in significantly shorter simulation time. The folded protein conformations were found within 0.2-2.1 Å of the native ...
Speckle Patterns and 2-Dimensional Brownian Motion
We present the results of a Monte Carlo simulation of Brownian Motion on a 2-dimensional lattice with nearest-neighbor interactions described by a linear model. These nearest-neighbor interactions lead to a spatial variance structure on the lattice. The resulting Brownian pattern fluctuates in value from point to point in a manner characteristic of a stationary stochastic process. The value at a lattice point is interpreted as an intensity level. The difference in values in neighboring cells produces a fluctuating intensity pattern on the lattice. Changing the size of the mesh changes the relative size of the speckles. Increasing the mesh size tends to average out the intensity in the direction of the mean of the stationary process. (Author)
Modeling collective emotions: a stochastic approach based on Brownian agents
We develop a agent-based framework to model the emergence of collective emotions, which is applied to online communities. Agents individual emotions are described by their valence and arousal. Using the concept of Brownian agents, these variables change according to a stochastic dynamics, which also considers the feedback from online communication. Agents generate emotional information, which is stored and distributed in a field modeling the online medium. This field affects the emotional states of agents in a non-linear manner. We derive conditions for the emergence of collective emotions, observable in a bimodal valence distribution. Dependent on a saturated or a super linear feedback between the information field and the agent's arousal, we further identify scenarios where collective emotions only appear once or in a repeated manner. The analytical results are illustrated by agent-based computer simulations. Our framework provides testable hypotheses about the emergence of collective emotions, which can be verified by data from online communities. (author)
Atomic dynamics of alumina melt: A molecular dynamics simulation study
Jahn, S.; P. A. Madden
2008-01-01
The atomic dynamics of Al2O3 melt are studied by molecular dynamics simulation. The particle interactions are described by an advanced ionic interaction model that includes polarization effects and ionic shape deformations. The model has been shown to reproduce accurately the static structure factors S(Q) from neutron and x-ray diffraction and the dynamic structure factor S(Q,ω) from inelastic x-ray scattering. Analysis of the partial dynamic structure factors shows inelastic features in the ...
Dynamics in inhomogeneous liquids and glasses via the test particle limit
Archer, Andrew J.; Hopkins, Paul; Schmidt, Matthias
2006-01-01
We show that one may view the self and the distinct part of the van Hove dynamic correlation function of a simple fluid as the one-body density distributions of a binary mixture that evolve in time according to dynamical density functional theory. For a test case of soft core Brownian particles the theory yields results for the van Hove function that agree quantitatively with those of our Brownian dynamics computer simulations. At sufficiently high densities the free energy landscape underlyi...
Particle dynamics in two-dimensional random-energy landscapes: Experiments and simulations
Evers, Florian; Zunke, Christoph; Hanes, Richard D. L.; Bewerunge, Jörg; Ladadwa, Imad; Heuer, Andreas; Egelhaaf, Stefan U.
2013-08-01
The dynamics of individual colloidal particles in random potential energy landscapes was investigated experimentally and by Monte Carlo simulations. The value of the potential at each point in the two-dimensional energy landscape follows a Gaussian distribution. The width of the distribution, and hence the degree of roughness of the energy landscape, was varied and its effect on the particle dynamics studied. This situation represents an example of Brownian dynamics in the presence of disorder. In the experiments, the energy landscapes were generated optically using a holographic setup with a spatial light modulator, and the particle trajectories were followed by video microscopy. The dynamics is characterized using, e.g., the time-dependent diffusion coefficient, the mean squared displacement, the van Hove function, and the non-Gaussian parameter. In both experiments and simulations the dynamics is initially diffusive, showing an extended subdiffusive regime at intermediate times before diffusive motion is recovered at very long times. The dependence of the long-time diffusion coefficient on the width of the Gaussian distribution agrees with theoretical predictions. Compared to the dynamics in a one-dimensional potential energy landscape, the localization at intermediate times is weaker and the diffusive regime at long times reached earlier, which is due to the possibility to avoid local maxima in two-dimensional energy landscapes.
Intrinsic and extrinsic measurement for Brownian motion
Based upon the Smoluchowski equation on curved manifolds, three physical observables are considered for Brownian displacement, namely geodesic displacement s, Euclidean displacement δR, and projected displacement δR⊥. The Weingarten–Gauss equations are used to calculate the mean-square Euclidean displacements in the short-time regime. Our findings show that from an extrinsic point of view the geometry of the space affects the Brownian motion in such a way that the particle’s diffusion is decelerated, contrasting with the intrinsic point of view where dynamics is controlled by the sign of the Gaussian curvature (Castro-Villarreal, 2010 J. Stat. Mech. P08006). Furthermore, it is possible to give exact formulas for 〈δR〉 and 〈δR2〉 on spheres and minimal surfaces, which are valid for all values of time. In the latter case, surprisingly, Brownian motion corresponds to the usual diffusion in flat geometries, albeit minimal surfaces have non-zero Gaussian curvature. Finally, the two-dimensional case is emphasized due to its close relation to surface self-diffusion in fluid membranes. (paper)
Brownian Motion Theory and Experiment
Basu, K; Basu, Kasturi; Baishya, Kopinjol
2003-01-01
Brownian motion is the perpetual irregular motion exhibited by small particles immersed in a fluid. Such random motion of the particles is produced by statistical fluctuations in the collisions they suffer with the molecules of the surrounding fluid. Brownian motion of particles in a fluid (like milk particles in water) can be observed under a microscope. Here we describe a simple experimental set-up to observe Brownian motion and a method of determining the diffusion coefficient of the Brownian particles, based on a theory due to Smoluchowski. While looking through the microscope we focus attention on a fixed small volume, and record the number of particles that are trapped in that volume, at regular intervals of time. This gives us a time-series data, which is enough to determine the diffusion coefficient of the particles to a good degree of accuracy.
Agusdinata, Datu Buyung, E-mail: bagusdinata@niu.edu; Amouie, Mahbod [Northern Illinois University, Department of Industrial & Systems Engineering and Environment, Sustainability, & Energy Institute (United States); Xu, Tao [Northern Illinois University, Department of Chemistry and Biochemistry (United States)
2015-01-15
Due to their favorable electrical and optical properties, quantum dots (QDs) nanostructures have found numerous applications including nanomedicine and photovoltaic cells. However, increased future production, use, and disposal of engineered QD products also raise concerns about their potential environmental impacts. The objective of this work is to establish a modeling framework for predicting the diffusion dynamics and concentration of toxic materials released from Trioctylphosphine oxide-capped CdSe. To this end, an agent-based model simulation with reaction kinetics and Brownian motion dynamics was developed. Reaction kinetics is used to model the stability of surface capping agent particularly due to oxidation process. The diffusion of toxic Cd{sup 2+} ions in aquatic environment was simulated using an adapted Brownian motion algorithm. A calibrated parameter to reflect sensitivity to reaction rate is proposed. The model output demonstrates the stochastic spatial distribution of toxic Cd{sup 2+} ions under different values of proxy environmental factor parameters. With the only chemistry considered was oxidation, the simulation was able to replicate Cd{sup 2+} ion release from Thiol-capped QDs in aerated water. The agent-based method is the first to be developed in the QDs application domain. It adds both simplicity of the solubility and rate of release of Cd{sup 2+} ions and complexity of tracking of individual atoms of Cd at the same time.
Due to their favorable electrical and optical properties, quantum dots (QDs) nanostructures have found numerous applications including nanomedicine and photovoltaic cells. However, increased future production, use, and disposal of engineered QD products also raise concerns about their potential environmental impacts. The objective of this work is to establish a modeling framework for predicting the diffusion dynamics and concentration of toxic materials released from Trioctylphosphine oxide-capped CdSe. To this end, an agent-based model simulation with reaction kinetics and Brownian motion dynamics was developed. Reaction kinetics is used to model the stability of surface capping agent particularly due to oxidation process. The diffusion of toxic Cd2+ ions in aquatic environment was simulated using an adapted Brownian motion algorithm. A calibrated parameter to reflect sensitivity to reaction rate is proposed. The model output demonstrates the stochastic spatial distribution of toxic Cd2+ ions under different values of proxy environmental factor parameters. With the only chemistry considered was oxidation, the simulation was able to replicate Cd2+ ion release from Thiol-capped QDs in aerated water. The agent-based method is the first to be developed in the QDs application domain. It adds both simplicity of the solubility and rate of release of Cd2+ ions and complexity of tracking of individual atoms of Cd at the same time
On collisions of Brownian particles
Ichiba, Tomoyuki; Karatzas, Ioannis
2010-01-01
We examine the behavior of $n$ Brownian particles diffusing on the real line with bounded, measurable drift and bounded, piecewise continuous diffusion coefficients that depend on the current configuration of particles. Sufficient conditions are established for the absence and for the presence of triple collisions among the particles. As an application to the Atlas model for equity markets, we study a special construction of such systems of diffusing particles using Brownian motions with refl...
Archimedes’ principle for Brownian liquid
Burdzy, Krzysztof; Chen, Zhen-Qing; Pal, Soumik
2011-01-01
We consider a family of hard core objects moving as independent Brownian motions confined to a vessel by reflection. These are subject to gravitational forces modeled by drifts. The stationary distribution for the process has many interesting implications, including an illustration of the Archimedes' principle. The analysis rests on constructing reflecting Brownian motion with drift in a general open connected domain and studying its stationary distribution. In dimension two we utilize known ...
Archimedes' principle for Brownian liquid
Burdzy, Krzysztof; Pal, Soumik
2009-01-01
We consider a family of hard core objects moving as independent Brownian motions confined to a vessel by reflection. These are subject to gravitational forces modeled by drifts. The stationary distribution for the process has many interesting implications, including an illustration of the Archimedes' principle. The analysis rests on constructing reflecting Brownian motion with drift in a general open connected domain and studying its stationary distribution. In dimension two we utilize known results about sphere packing.
A qualitative dynamical model for cardiotocography simulation
Illanes, Alfredo; Haritopoulos, Michel; Robles, Felipe; Guerra, Francisco
2015-01-01
The purpose of this work is to present a new mathematical model for fetal monitoring simulation. It involves the simultaneous generation of fetal heart rate and maternal uterine contraction signals through a parametrical model. This model allows the generation of the main fetal monitoring dynamics including fetal movements, acceleration and deceleration of the heart rate and the dynami-cal adjustment of fetal heart rate following an uterine contraction. Simulated tracings were analyzed by spe...
Molecular Dynamic Simulation on High Performance Infrastrucutres
Bergant, Anže
2016-01-01
This thesis covers comparison between different computer platforms of high performance computing while performing molecular dynamics simulations, which falls under very complex problems and needs lots of processing power. Our goal was to critically evaluate different platforms while solving molecular dynamics, so we used 1 to 16 processor cores on a computer cluster and one and two graphics processing units (GPU) for simulations. The results will be used while planning on buying new computer ...
Visualizing Structure and Dynamics of Disaccharide Simulations
Matthews, J. F.; Beckham, G. T.; Himmel, M. E.; Crowley, M. F.
2012-01-01
We examine the effect of several solvent models on the conformational properties and dynamics of disaccharides such as cellobiose and lactose. Significant variation in timescale for large scale conformational transformations are observed. Molecular dynamics simulation provides enough detail to enable insight through visualization of multidimensional data sets. We present a new way to visualize conformational space for disaccharides with Ramachandran plots.
Molecular dynamics simulation of impact test
Akahoshi, Y. [Kyushu Inst. of Tech., Kitakyushu, Fukuoka (Japan); Schmauder, S.; Ludwig, M. [Stuttgart Univ. (Germany). Staatliche Materialpruefungsanstalt
1998-11-01
This paper describes an impact test by molecular dynamics (MD) simulation to evaluate embrittlement of bcc Fe at different temperatures. A new impact test model is developed for MD simulation. The typical fracture behaviors show transition from brittle to ductile fracture, and a history of the impact loads also demonstrates its transition. We conclude that the impact test by MD could be feasible. (orig.)
Molecular Dynamics Simulations of Simple Liquids
Speer, Owner F.; Wengerter, Brian C.; Taylor, Ramona S.
2004-01-01
An experiment, in which students were given the opportunity to perform molecular dynamics simulations on a series of molecular liquids using the Amber suite of programs, is presented. They were introduced to both physical theories underlying classical mechanics simulations and to the atom-atom pair distribution function.
Melting behaviour of gold-platinum nanoalloy clusters by molecular dynamics simulations
Ong, Yee Pin; Yoon, Tiem Leong [School of Physics, Universiti Sains Malaysia, 11800 USM, Penang (Malaysia); Lim, Thong Leng [Faculty of Engineering and Technology, Multimedia University, Melaka Campus, 75450 Melaka (Malaysia)
2015-04-24
The melting behavior of bimetallic gold-platinum nanoclusters is studied by applying Brownian-type isothermal molecular dynamics (MD) simulation, a program modified from the cubic coupling scheme (CCS). The process begins with the ground-state structures obtained from global minimum search algorithm and proceeds with the investigation of the effect of temperature on the thermal properties of gold-platinum nanoalloy clusters. N-body Gupta potential has been employed in order to account for the interactions between gold and platinum atoms. The ground states of the nanoalloy clusters, which are core-shell segregated, are heated until they become thermally segregated. The detailed melting mechanism of the nanoalloy clusters is studied via this approach to provide insight into the thermal stability of the nanoalloy clusters.
Melting behaviour of gold-platinum nanoalloy clusters by molecular dynamics simulations
The melting behavior of bimetallic gold-platinum nanoclusters is studied by applying Brownian-type isothermal molecular dynamics (MD) simulation, a program modified from the cubic coupling scheme (CCS). The process begins with the ground-state structures obtained from global minimum search algorithm and proceeds with the investigation of the effect of temperature on the thermal properties of gold-platinum nanoalloy clusters. N-body Gupta potential has been employed in order to account for the interactions between gold and platinum atoms. The ground states of the nanoalloy clusters, which are core-shell segregated, are heated until they become thermally segregated. The detailed melting mechanism of the nanoalloy clusters is studied via this approach to provide insight into the thermal stability of the nanoalloy clusters
Noise-to-signal transition of a Brownian particle in the cubic potential: I. general theory
Filip, Radim; Zemánek, Pavel
2016-06-01
The noise-to-signal transitions are very interesting processes in physics, as they might transform environmental noise to useful mechanical effects. We theoretically analyze stochastic noise-to-signal transition of overdamped Brownian motion of a particle in the cubic potential. The particle reaches thermal equilibrium with its environment in the quadratic potential which is suddenly swapped to the cubic potential. We predict a simultaneous increase of both the displacement and signal-to-noise ratio in the cubic potential for the position linearly powered by the temperature of the particle environment. The short-time analysis and numerical simulations fully confirm different dynamical regimes of this noise-to-signal transition.
Anomalous Brownian refrigerator
Rana, Shubhashis; Pal, P. S.; Saha, Arnab; Jayannavar, A. M.
2016-02-01
We present a detailed study of a Brownian particle driven by Carnot-type refrigerating protocol operating between two thermal baths. Both the underdamped as well as the overdamped limits are investigated. The particle is in a harmonic potential with time-periodic strength that drives the system cyclically between the baths. Each cycle consists of two isothermal steps at different temperatures and two adiabatic steps connecting them. Besides working as a stochastic refrigerator, it is shown analytically that in the quasistatic regime the system can also act as stochastic heater, depending on the bath temperatures. Interestingly, in non-quasistatic regime, our system can even work as a stochastic heat engine for certain range of cycle time and bath temperatures. We show that the operation of this engine is not reliable. The fluctuations of stochastic efficiency/coefficient of performance (COP) dominate their mean values. Their distributions show power law tails, however the exponents are not universal. Our study reveals that microscopic machines are not the microscopic equivalent of the macroscopic machines that we come across in our daily life. We find that there is no one to one correspondence between the performance of our system under engine protocol and its reverse.
Martínez, I. A.; Roldán, É.; Dinis, L.; Petrov, D.; Parrondo, J. M. R.; Rica, R. A.
2016-01-01
The Carnot cycle imposes a fundamental upper limit to the efficiency of a macroscopic motor operating between two thermal baths. However, this bound needs to be reinterpreted at microscopic scales, where molecular bio-motors and some artificial micro-engines operate. As described by stochastic thermodynamics, energy transfers in microscopic systems are random and thermal fluctuations induce transient decreases of entropy, allowing for possible violations of the Carnot limit. Here we report an experimental realization of a Carnot engine with a single optically trapped Brownian particle as the working substance. We present an exhaustive study of the energetics of the engine and analyse the fluctuations of the finite-time efficiency, showing that the Carnot bound can be surpassed for a small number of non-equilibrium cycles. As its macroscopic counterpart, the energetics of our Carnot device exhibits basic properties that one would expect to observe in any microscopic energy transducer operating with baths at different temperatures. Our results characterize the sources of irreversibility in the engine and the statistical properties of the efficiency--an insight that could inspire new strategies in the design of efficient nano-motors.
Molecular dynamics simulations using graphics processing units
Baker, J.A.; Hirst, J.D.
2011-01-01
It is increasingly easy to develop software that exploits Graphics Processing Units (GPUs). The molecular dynamics simulation community has embraced this recent opportunity. Herein, we outline the current approaches that exploit this technology. In the context of biomolecular simulations, we discuss some of the algorithms that have been implemented and some of the aspects that distinguish the GPU from previous parallel environments. The ubiquity of GPUs and the ingenuity of the simulation com...
Generalization of Brownian Motion with Autoregressive Increments
Fendick, Kerry
2011-01-01
This paper introduces a generalization of Brownian motion with continuous sample paths and stationary, autoregressive increments. This process, which we call a Brownian ray with drift, is characterized by three parameters quantifying distinct effects of drift, volatility, and autoregressiveness. A Brownian ray with drift, conditioned on its state at the beginning of an interval, is another Brownian ray with drift over the interval, and its expected path over the interval is a ray with a slope that depends on the conditioned state. This paper shows how Brownian rays can be applied in finance for the analysis of queues or inventories and the valuation of options. We model a queue's net input process as a superposition of Brownian rays with drift and derive the transient distribution of the queue length conditional on past queue lengths and on past states of the individual Brownian rays comprising the superposition. The transient distributions of Regulated Brownian Motion and of the Regulated Brownian Bridge are...
Nanofluidic Brownian Ratchet via atomically-stepped surfaces
Rahmani, Amir; Colosqui, Carlos
2015-11-01
Theoretical analysis and fully atomistic molecular dynamics simulations reveal a Brownian ratchet mechanism by which thermal motion can drive the directional displacement of liquids confined in micro- or nanoscale channels and pores. The particular systems discussed in this talk consist of two immiscible liquids confined in a slit-like nanochannel with atomically-stepped surfaces. Mean displacement rates reported in molecular dynamics simulations are in close agreement with theoretical predictions via analytical solution of a Smoluchowski equation for the probability density of the position of the liquid-liquid interface. The direction of the thermally-driven displacement of liquid is determined by the nanostructure surface geometry and thus imbibition or drainage can occur against the direction of action of capillary forces. The studied surface nanostructure with directional asymmetry can control the dynamics of wetting processes such as capillary filling, wicking, and imbibition in porous materials. The proposed physical mechanisms and derived analytical expressions can be applied to design nanofluidic and microfluidic devices for passive handling and separation.
Discrete dislocation dynamics simulations in a cylinder
Li, Maosheng; Gao, Chan; Xu, Jianing
2015-02-01
Mechanical properties of material are closely related to the motion of dislocations, and predicting the interactions and resulting collective motion of dislocations is a major task in understanding and modelling plastically deforming materials. A discrete dislocation dynamics model is used to describe the orientation substructure within the microstructure. Discrete dislocation dynamics simulations in three dimensions have been used to examine the role of dislocation multiplication and mobility on the plasticity in small samples under uniaxial compression. In this paper we describe the application of the dislocation dynamics simulations in a cylindrical geometry. The boundary conditions for the simulation were estimated from the distribution of the geometrically necessary dislocation density which was obtained from the orientation map. Numerical studies benchmark could validate the accuracy of the algorithms and the importance of handling the singularity correctly. The results of the simulation explain the formation of the experimentally observed substructure.
Spin dynamics simulations at AGS
Huang, H.; MacKay, W.W.; Meot, F.; Roser, T.
2010-05-23
To preserve proton polarization through acceleration, it is important to have a correct model of the process. It has been known that with the insertion of the two helical partial Siberian snakes in the Alternating Gradient Synchrotron (AGS), the MAD model of AGS can not deal with a field map with offset orbit. The stepwise ray-tracing code Zgoubi provides a tool to represent the real electromagnetic fields in the modeling of the optics and spin dynamics for the AGS. Numerical experiments of resonance crossing, including spin dynamics in presence of the snakes and Q-jump, have been performed in AGS lattice models, using Zgoubi. This contribution reports on various results so obtained.
Harmonic functions on Walsh's Brownian motion
Jehring, Kristin Elizabeth
2009-01-01
In this dissertation we examine a variation of two- dimensional Brownian motion introduced in 1978 by Walsh. Walsh's Brownian motion can be described as a Brownian motion on the spokes of a (rimless) bicycle wheel. We will construct such a process by randomly assigning an angle to the excursions of a reflecting Brownian motion from 0. With this construction we see that Walsh's Brownian motion in R² behaves like one-dimensional Brownian motion away from the origin, but at the origin behaves di...
Dynamics of individual colloidal particles in one-dimensional random potentials: a simulation study.
Hanes, Richard D L; Egelhaaf, Stefan U
2012-11-21
Using Monte Carlo simulations, individual Brownian particles have been investigated in a one-dimensional random energy landscape whose energy levels are selected from a Gaussian distribution. The standard deviation of the distribution determines the roughness of the noise-like potential and was varied in the simulations. After initialization, which was done by an instantaneous or infinitely slow (annealed) quench, the particle dynamics were followed. They were characterized by a number of parameters, such as the mean squared displacement, the time dependent diffusion coefficient, the non-Gaussian parameter, and the van Hove function. The dynamics exhibit different regimes: at very short times superdiffusion, followed by normal diffusion, and subsequently an extended period of subdiffusive dynamics due to localization within the minima of the potential, and finally, after a very slow approach towards the long-time limit, again diffusion with a significantly reduced diffusion coefficient. The long-time diffusion coefficient is consistent with theoretical predictions while no predictions exist for the intermediate times. Nevertheless, over the whole time range, the simulation results are in agreement with recent experimental findings on colloidal particles in a random potential created by a holographic optical setup. PMID:23114221
Híjar, Humberto
2015-02-01
We study the Brownian motion of a particle bound by a harmonic potential and immersed in a fluid with a uniform shear flow. We describe this problem first in terms of a linear Fokker-Planck equation which is solved to obtain the probability distribution function for finding the particle in a volume element of its associated phase space. We find the explicit form of this distribution in the stationary limit and use this result to show that both the equipartition law and the equation of state of the trapped particle are modified from their equilibrium form by terms increasing as the square of the imposed shear rate. Subsequently, we propose an alternative description of this problem in terms of a generalized Langevin equation that takes into account the effects of hydrodynamic correlations and sound propagation on the dynamics of the trapped particle. We show that these effects produce significant changes, manifested as long-time tails and resonant peaks, in the equilibrium and nonequilibrium correlation functions for the velocity of the Brownian particle. We implement numerical simulations based on molecular dynamics and multiparticle collision dynamics, and observe a very good quantitative agreement between the predictions of the model and the numerical results, thus suggesting that this kind of numerical simulations could be used as complement of current experimental techniques. PMID:25768490
Ilg, Patrick; Karlin, Iliya V.
2004-01-01
A method for the combination of microscopic and macroscopic simulations is developed which is based on the invariance of the macroscopic relative to the microscopic dynamics. The method recognizes the onset and breakdown of the macroscopic description during the integration. We apply this method to the case of ferrofluid dynamics, where it switches between direct Brownian dynamics simulations and integration of the constitutive equation.
Sensitivity Analysis of Fire Dynamics Simulation
Brohus, Henrik; Nielsen, Peter V.; Petersen, Arnkell J.; Sommerlund-Larsen, Kim
2007-01-01
In case of fire dynamics simulation requirements to reliable results are most often very high due to the severe consequences of erroneous results. At the same time it is a well known fact that fire dynamics simulation constitutes rather complex physical phenomena which apart from flow and energy...... equations require solution of the issues of combustion and gas radiation to mention a few. This paper performs a sensitivity analysis of a fire dynamics simulation on a benchmark case where measurement results are available for comparison. The analysis is performed using the method of Elementary Effects...... (Morris method). The parameters considered are selected among physical parameters and program specific parameters. The influence on the calculation result as well as the CPU time is considered. It is found that the result is highly sensitive to many parameters even though the sensitivity varies...
On the Generalized Brownian Motion and its Applications in Finance
Høg, Esben; Frederiksen, Per; Schiemert, Daniel
This paper deals with dynamic term structure models (DTSMs) and proposes a new way to handle the limitation of the classical affine models. In particular, the paper expands the exibility of the DTSMs by applying generalized Brownian motions with dependent increments as the governing force of the ...
Dynamical structure of fluid mercury: Molecular-dynamics simulations
Hoshino, Kozo; Tanaka, Shunichiro; Shimojo, Fuyuki
2007-01-01
We have carried out molecular-dynamics simulations for nonmetallic fluid mercury in liquid and vapor phases using a Lennard-Jones type effective potential and shown that the structure factors S(Q) and the dynamic structure factors S(Q, omega) of nonmetallic fluid mercury obtained by our MD simulations are in good agreement with recent X-ray diffraction and inelastic X-ray scattering experiments. We conclude from these results that, though the fluid mercury which shows a metal-nonmetal transit...
Entropic forces in Brownian motion
Roos, Nico
2013-01-01
The interest in the concept of entropic forces has risen considerably since E. Verlinde proposed to interpret the force in Newton s second law and Gravity as entropic forces. Brownian motion, the motion of a small particle (pollen) driven by random impulses from the surrounding molecules, may be the first example of a stochastic process in which such forces are expected to emerge. In this note it is shown that at least two types of entropic motion can be identified in the case of 3D Brownian motion (or random walk). This yields simple derivations of known results of Brownian motion, Hook s law and, applying an external (nonradial) force, Curie s law and the Langevin-Debye equation.
Entropic forces in Brownian motion
Roos, Nico
2014-12-01
Interest in the concept of entropic forces has risen considerably since Verlinde proposed in 2011 to interpret the force in Newton's second law and gravity as entropic forces. Brownian motion—the motion of a small particle (pollen) driven by random impulses from the surrounding molecules—may be the first example of a stochastic process in which such forces are expected to emerge. In this article, it is shown that at least two types of entropic force can be identified in three-dimensional Brownian motion. This analysis yields simple derivations of known results of Brownian motion, Hooke's law, and—applying an external (non-radial) force—Curie's law and the Langevin-Debye equation.
Methane in carbon nanotube - molecular dynamics simulation
Bartuś, Katarzyna; Bródka, Aleksander
2011-01-01
Abstract The behaviour of methane molecules inside carbon nanotube at room temperature is studied using classical molecular dynamics simulations. A methane molecule is represented either by a shapeless super-atom or by rigid set of 5 interaction centres localised on atoms. Different loadings of methane molecules ranging from the dense gas density to the liquid density, and the influence of flexibility of the CNT on structural and dynamics properties of confined molecules are consid...
Dynamic Fracture Simulations of Explosively Loaded Cylinders
Arthur, Carly W. [Univ. of California, Davis, CA (United States). Dept. of Civil and Environmental Engineering; Goto, D. M. [Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
2015-11-30
This report documents the modeling results of high explosive experiments investigating dynamic fracture of steel (AerMet® 100 alloy) cylinders. The experiments were conducted at Lawrence Livermore National Laboratory (LLNL) during 2007 to 2008 [10]. A principal objective of this study was to gain an understanding of dynamic material failure through the analysis of hydrodynamic computer code simulations. Two-dimensional and three-dimensional computational cylinder models were analyzed using the ALE3D multi-physics computer code.
Cyclotron beam dynamic simulations in MATLAB
MATLAB is useful for beam dynamic simulations in cyclotrons. Programming in an easy-to-use environment permits creation of models in a short space of time. Advanced graphical tools of MATLAB give good visualization features to created models. The beam dynamic modeling results with an example of two different cyclotron designs are presented. Programming with MATLAB opens wide possibilities of the development of the complex program, able to perform complete block of calculations for the design of the accelerators
Quantum Brownian motion model for the stock market
Meng, Xiangyi; Zhang, Jian-Wei; Guo, Hong
2016-06-01
It is believed by the majority today that the efficient market hypothesis is imperfect because of market irrationality. Using the physical concepts and mathematical structures of quantum mechanics, we construct an econophysical framework for the stock market, based on which we analogously map massive numbers of single stocks into a reservoir consisting of many quantum harmonic oscillators and their stock index into a typical quantum open system-a quantum Brownian particle. In particular, the irrationality of stock transactions is quantitatively considered as the Planck constant within Heisenberg's uncertainty relationship of quantum mechanics in an analogous manner. We analyze real stock data of Shanghai Stock Exchange of China and investigate fat-tail phenomena and non-Markovian behaviors of the stock index with the assistance of the quantum Brownian motion model, thereby interpreting and studying the limitations of the classical Brownian motion model for the efficient market hypothesis from a new perspective of quantum open system dynamics.
Simulating coronal condensation dynamics in 3D
Moschou, S P; Xia, C; Fang, X
2015-01-01
We present numerical simulations in 3D settings where coronal rain phenomena take place in a magnetic configuration of a quadrupolar arcade system. Our simulation is a magnetohydrodynamic simulation including anisotropic thermal conduction, optically thin radiative losses, and parametrised heating as main thermodynamical features to construct a realistic arcade configuration from chromospheric to coronal heights. The plasma evaporation from chromospheric and transition region heights eventually causes localised runaway condensation events and we witness the formation of plasma blobs due to thermal instability, that evolve dynamically in the heated arcade part and move gradually downwards due to interchange type dynamics. Unlike earlier 2.5D simulations, in this case there is no large scale prominence formation observed, but a continuous coronal rain develops which shows clear indications of Rayleigh-Taylor or interchange instability, that causes the denser plasma located above the transition region to fall do...
Dynamic Simulations of Advanced Fuel Cycles
Steven J. Piet; Brent W. Dixon; Jacob J. Jacobson; Gretchen E. Matthern; David E. Shropshire
2011-03-01
Years of performing dynamic simulations of advanced nuclear fuel cycle options provide insights into how they could work and how one might transition from the current once-through fuel cycle. This paper summarizes those insights from the context of the 2005 objectives and goals of the U.S. Advanced Fuel Cycle Initiative (AFCI). Our intent is not to compare options, assess options versus those objectives and goals, nor recommend changes to those objectives and goals. Rather, we organize what we have learned from dynamic simulations in the context of the AFCI objectives for waste management, proliferation resistance, uranium utilization, and economics. Thus, we do not merely describe “lessons learned” from dynamic simulations but attempt to answer the “so what” question by using this context. The analyses have been performed using the Verifiable Fuel Cycle Simulation of Nuclear Fuel Cycle Dynamics (VISION). We observe that the 2005 objectives and goals do not address many of the inherently dynamic discriminators among advanced fuel cycle options and transitions thereof.
Massively-Parallel Dislocation Dynamics Simulations
Cai, W; Bulatov, V V; Pierce, T G; Hiratani, M; Rhee, M; Bartelt, M; Tang, M
2003-06-18
Prediction of the plastic strength of single crystals based on the collective dynamics of dislocations has been a challenge for computational materials science for a number of years. The difficulty lies in the inability of the existing dislocation dynamics (DD) codes to handle a sufficiently large number of dislocation lines, in order to be statistically representative and to reproduce experimentally observed microstructures. A new massively-parallel DD code is developed that is capable of modeling million-dislocation systems by employing thousands of processors. We discuss the general aspects of this code that make such large scale simulations possible, as well as a few initial simulation results.
Winding statistics of a Brownian particle on a ring
We consider a Brownian particle moving on a ring. We study the probability distributions of the total number of turns and the net number of counter-clockwise turns the particle makes until time t. Using a method based on the renewal properties of a Brownian walker, we find exact analytical expressions of these distributions. This method serves as an alternative to the standard path integral techniques which are not always easily adaptable for certain observables. For large t, we show that these distributions have Gaussian scaling forms. We also compute large deviation functions associated to these distributions characterizing atypically large fluctuations. We provide numerical simulations in support of our analytical results. (paper)
Dynamic simulation of a high efficiency building
Dot Bardolet, Núria
2009-01-01
The dynamic simulation of the heat transfer between buildings and their environment is situdied through the software TRNSYS. First, a very simple building is considered and the energy demand for heating evaluated through TRNSYS is compared with that obtained by apllying the Italian Standard UNI TS 11300-1:2008. This excersise has the purpose of verifying the correct use of the simulation program. Then, the annual energy demand for heating, cooling and humidity control, with reference to a wel...
Feynman Rules For Brownian Motion
Hatamian, S T
2003-01-01
We present a perturbation theory extending a prescription due to Feynman for computing the probability density function in Brownian-motion. The method used, can be applied to a wide variety of otherwise difficult circumstances in Brownian-motion. The exact moments and kurtosis, if not the distribution itself, for many important cases can be summed for arbitrary times. As expected, the behavior at early time regime, for the sample processes considered, deviate significantly from the usual diffusion theory; a fact with important consequences in various applications such as financial physics. A new class of functions dubbed "Damped-exponential-integrals" are also identified.
Brownian movement and molecular reality
Perrin, Jean
2005-01-01
How do we know that molecules really exist? An important clue came from Brownian movement, a concept developed in 1827 by botanist Robert Brown, who noticed that tiny objects like pollen grains shook and moved erratically when viewed under a microscope. Nearly 80 years later, in 1905, Albert Einstein explained this ""Brownian motion"" as the result of bombardment by molecules. Einstein offered a quantitative explanation by mathematically estimating the average distance covered by the particles over time as a result of molecular bombardment. Four years later, Jean Baptiste Perrin wrote Brownia
Brownian molecular rotors: Theoretical design principles and predicted realizations
Schönborn, Jan Boyke; Herges, Rainer; Hartke, Bernd
2009-01-01
We propose simple design concepts for molecular rotors driven by Brownian motion and external photochemical switching. Unidirectionality and efﬁciency of the motion is measured by explicit simulations. Two different molecular scaffolds are shown to yield viable molecular rotors when decorated with suitable substituents.
Brownian colloidal particles: Ito, Stratonovich, or a different stochastic interpretation
Sancho, J. M.
2011-12-01
Recent experiments on Brownian colloidal particles have been studied theoretically in terms of overdamped Langevin equations with multiplicative white noise using an unconventional stochastic interpretation. Complementary numerical simulations of the same system are well described using the conventional Stratonovich interpretation. Here we address this dichotomy from a more generic starting point: the underdamped Langevin equation and its corresponding Fokker-Planck equation.
The valuation of currency options by fractional Brownian motion.
Shokrollahi, Foad; Kılıçman, Adem
2016-01-01
This research aims to investigate a model for pricing of currency options in which value governed by the fractional Brownian motion model (FBM). The fractional partial differential equation and some Greeks are also obtained. In addition, some properties of our pricing formula and simulation studies are presented, which demonstrate that the FBM model is easy to use. PMID:27504243
Induction generator models in dynamic simulation tools
Knudsen, Hans; Akhmatov, Vladislav
1999-01-01
found to be possible to include a transient model in dynamic stability tools and, then, obtain correct results also in dynamic tools. The representation of the rotating system influences on the voltage recovery shape which is an important observation in case of windmills, where a heavy mill is connected......For AC network with large amount of induction generators (windmills) the paper demonstrates a significant discrepancy in the simulated voltage recovery after fault in weak networks when comparing dynamic and transient stability descriptions and the reasons of discrepancies are explained. It is...
Some Brownian functionals and their laws
Donati-Martin, C.; Yor, M.
1997-01-01
We develop some topics about Brownian motion with a particular emphasis on the study of principal values of Brownian local times. We show some links between principal values and Doob’s $h$-transforms of Brownian motion, for nonpositive harmonic functions $h$. We also give a survey and complement some martingale approaches to Ray–Knight theorems for local times.
Fluid Dynamics Theory, Computation, and Numerical Simulation
Pozrikidis, Constantine
2009-01-01
Fluid Dynamics: Theory, Computation, and Numerical Simulation is the only available book that extends the classical field of fluid dynamics into the realm of scientific computing in a way that is both comprehensive and accessible to the beginner. The theory of fluid dynamics, and the implementation of solution procedures into numerical algorithms, are discussed hand-in-hand and with reference to computer programming. This book is an accessible introduction to theoretical and computational fluid dynamics (CFD), written from a modern perspective that unifies theory and numerical practice. There are several additions and subject expansions in the Second Edition of Fluid Dynamics, including new Matlab and FORTRAN codes. Two distinguishing features of the discourse are: solution procedures and algorithms are developed immediately after problem formulations are presented, and numerical methods are introduced on a need-to-know basis and in increasing order of difficulty. Matlab codes are presented and discussed for ...
Fluid dynamics theory, computation, and numerical simulation
Pozrikidis, C
2001-01-01
Fluid Dynamics Theory, Computation, and Numerical Simulation is the only available book that extends the classical field of fluid dynamics into the realm of scientific computing in a way that is both comprehensive and accessible to the beginner The theory of fluid dynamics, and the implementation of solution procedures into numerical algorithms, are discussed hand-in-hand and with reference to computer programming This book is an accessible introduction to theoretical and computational fluid dynamics (CFD), written from a modern perspective that unifies theory and numerical practice There are several additions and subject expansions in the Second Edition of Fluid Dynamics, including new Matlab and FORTRAN codes Two distinguishing features of the discourse are solution procedures and algorithms are developed immediately after problem formulations are presented, and numerical methods are introduced on a need-to-know basis and in increasing order of difficulty Matlab codes are presented and discussed for a broad...
Zimmermann, Urs; Smallenburg, Frank; Löwen, Hartmut
2016-06-22
Using both dynamical density functional theory and particle-resolved Brownian dynamics simulations, we explore the flow of two-dimensional colloidal solids and fluids driven through a linear channel with a constriction. The flow is generated by a constant external force acting on all colloids. The initial configuration is equilibrated in the absence of flow and then the external force is switched on instantaneously. Upon starting the flow, we observe four different scenarios: a complete blockade, a monotonic decay to a constant particle flux (typical for a fluid), a damped oscillatory behaviour in the particle flux, and a long-lived stop-and-go behaviour in the flow (typical for a solid). The dynamical density functional theory describes all four situations but predicts infinitely long undamped oscillations in the flow which are always damped in the simulations. We attribute the mechanisms of the underlying stop-and-go flow to symmetry conditions on the flowing solid. Our predictions are verifiable in real-space experiments on magnetic colloidal monolayers which are driven through structured microchannels and can be exploited to steer the flow throughput in microfluidics. PMID:27116706
Monte Carlo simulation for soot dynamics
Zhou Kun
2012-01-01
Full Text Available A new Monte Carlo method termed Comb-like frame Monte Carlo is developed to simulate the soot dynamics. Detailed stochastic error analysis is provided. Comb-like frame Monte Carlo is coupled with the gas phase solver Chemkin II to simulate soot formation in a 1-D premixed burner stabilized flame. The simulated soot number density, volume fraction, and particle size distribution all agree well with the measurement available in literature. The origin of the bimodal distribution of particle size distribution is revealed with quantitative proof.
Monte carlo simulation for soot dynamics
Zhou, Kun
2012-01-01
A new Monte Carlo method termed Comb-like frame Monte Carlo is developed to simulate the soot dynamics. Detailed stochastic error analysis is provided. Comb-like frame Monte Carlo is coupled with the gas phase solver Chemkin II to simulate soot formation in a 1-D premixed burner stabilized flame. The simulated soot number density, volume fraction, and particle size distribution all agree well with the measurement available in literature. The origin of the bimodal distribution of particle size distribution is revealed with quantitative proof.
Quantum Brownian motion in a Landau level
Cobanera, E.; Kristel, P.; Morais Smith, C.
2016-06-01
Motivated by questions about the open-system dynamics of topological quantum matter, we investigated the quantum Brownian motion of an electron in a homogeneous magnetic field. When the Fermi length lF=ℏ /(vFmeff) becomes much longer than the magnetic length lB=(ℏc /e B ) 1 /2 , then the spatial coordinates X ,Y of the electron cease to commute, [X ,Y ] =i lB2 . As a consequence, localization of the electron becomes limited by Heisenberg uncertainty, and the linear bath-electron coupling becomes unconventional. Moreover, because the kinetic energy of the electron is quenched by the strong magnetic field, the electron has no energy to give to or take from the bath, and so the usual connection between frictional forces and dissipation no longer holds. These two features make quantum Brownian motion topological, in the regime lF≫lB , which is at the verge of current experimental capabilities. We model topological quantum Brownian motion in terms of an unconventional operator Langevin equation derived from first principles, and solve this equation with the aim of characterizing diffusion. While diffusion in the noncommutative plane turns out to be conventional, with the mean displacement squared being proportional to tα and α =1 , there is an exotic regime for the proportionality constant in which it is directly proportional to the friction coefficient and inversely proportional to the square of the magnetic field: in this regime, friction helps diffusion and the magnetic field suppresses all fluctuations. We also show that quantum tunneling can be completely suppressed in the noncommutative plane for suitably designed metastable potential wells, a feature that might be worth exploiting for storage and protection of quantum information.
Molecular dynamic simulation of directional crystal growth
Costa, B. V.; Coura, P. Z.; Mesquita, O. N.
1999-01-01
We use molecular dynamic to simulate the directional growth of binary mixtures. our results compare very well with analitical and experimental results. This opens up the possibility to probe growth situations which are difficult to reach experimentally, being an important tool for further experimental and theoretical developments in the area of crystal growth.
Molecular dynamics - NMR experiments and simulations
Dračínský, Martin; Hodgkinson, P.; Kessler, Jiří; Bouř, Petr
Brno : Masaryk University Press, 2015 - (Sklenář, V.). s. 277-278 ISBN 978-80-210-7890-1. [EUROMAR 2015. 05.07.2015-10.07.2015, Praha] Institutional support: RVO:61388963 Keywords : molecular dynamics * NMR spectroscopy * MD simulations Subject RIV: CF - Physical ; Theoretical Chemistry
Object Oriented Modelling and Dynamical Simulation
Wagner, Falko Jens; Poulsen, Mikael Zebbelin
1998-01-01
This report with appendix describes the work done in master project at DTU.The goal of the project was to develop a concept for simulation of dynamical systems based on object oriented methods.The result was a library of C++-classes, for use when both building componentbased models and when...
Dynamic simulation of a reverse Brayton refrigerator
Peng, N.; Xiong, L. Y.; Dong, B.; Liu, L. Q. [State Key Laboratory of Technologies in Space Cryogenic Propellants, Technical Institute of Physics and Chemistry, CAS, Beijing, 100190 (China); Lei, L. L.; Tang, J. C. [State Key Laboratory of Technologies in Space Cryogenic Propellants, Technical Institute of Physics and Chemistry, CAS, Beijing, 100190 China and Graduate University of Chinese Academy of Sciences, Beijing, 100190 (China)
2014-01-29
A test refrigerator based on the modified Reverse Brayton cycle has been developed in the Chinese Academy of Sciences recently. To study the behaviors of this test refrigerator, a dynamic simulation has been carried out. The numerical model comprises the typical components of the test refrigerator: compressor, valves, heat exchangers, expander and heater. This simulator is based on the oriented-object approach and each component is represented by a set of differential and algebraic equations. The control system of the test refrigerator is also simulated, which can be used to optimize the control strategies. This paper describes all the models and shows the simulation results. Comparisons between simulation results and experimental data are also presented. Experimental validation on the test refrigerator gives satisfactory results.
Dynamic modeling and simulation of wind turbines
Using wind energy for generating electricity in wind turbines is a good way for using renewable energies. It can also help to protect the environment. The main objective of this paper is dynamic modeling by energy method and simulation of a wind turbine aided by computer. In this paper, the equations of motion are extracted for simulating the system of wind turbine and then the behavior of the system become obvious by solving the equations. The turbine is considered with three blade rotor in wind direction, induced generator that is connected to the network and constant revolution for simulation of wind turbine. Every part of the wind turbine should be simulated for simulation of wind turbine. The main parts are blades, gearbox, shafts and generator
Dynamic simulation of a reverse Brayton refrigerator
A test refrigerator based on the modified Reverse Brayton cycle has been developed in the Chinese Academy of Sciences recently. To study the behaviors of this test refrigerator, a dynamic simulation has been carried out. The numerical model comprises the typical components of the test refrigerator: compressor, valves, heat exchangers, expander and heater. This simulator is based on the oriented-object approach and each component is represented by a set of differential and algebraic equations. The control system of the test refrigerator is also simulated, which can be used to optimize the control strategies. This paper describes all the models and shows the simulation results. Comparisons between simulation results and experimental data are also presented. Experimental validation on the test refrigerator gives satisfactory results
Molecular Dynamics Simulations of Janus Particle Dynamics in Uniform Flow
Archereau, Aurelien Y M; Willmott, Geoff R
2016-01-01
We use molecular dynamics simulations to study the dynamics of Janus particles, micro- or nanoparticles which are not spherically symmetric, in the uniform flow of a simple liquid. In particular we consider spheres with an asymmetry in the solid-liquid interaction over their surfaces and calculate the forces and torques experienced by the particles as a function of their orientation with respect to the flow. We also examine particles that are deformed slightly from a spherical shape. We compare the simulation results to the predictions of a previously introduced theoretical approach, which computes the forces and torques on particles with variable slip lengths or aspherical deformations that are much smaller than the particle radius. We find that there is good agreement between the forces and torques computed from our simulations and the theoretical predictions, when the slip condition is applied to the first layer of liquid molecules adjacent to the surface.
Computer simulation of confined liquid crystal dynamics
Webster, R E
2001-01-01
are performed of the formation of structures in confined smectic systems where layer tilt is induced by an imposed surface pretilt. Results show that bookshelf, chevron and tilled layer structures are observable in a confined Gay-Berne system. The formation and stability of the chevron structure are shown to be influenced by surface slip. Results are presented from a series of simulations undertaken to determine whether dynamic processes observed in device-scale liquid crystal cells confined between aligning substrates can be simulated in a molecular system using parallel molecular dynamics of the Gay-Berne model. In a nematic cell, on removal of an aligning field, initial near-surface director relaxation can induce flow, termed 'backflow' in the liquid. This, in turn, can cause director rotation, termed 'orientational kickback', in the centre of the cell. Simulations are performed of the relaxation in nematic systems confined between substrates with a common alignment on removal of an aligning field. Results...
Probing Cellular Dynamics with Mesoscopic Simulations
Shillcock, Julian C.
2010-01-01
Cellular processes span a huge range of length and time scales from the molecular to the near-macroscopic. Understanding how effects on one scale influence, and are themselves influenced by, those on lower and higher scales is a critical issue for the construction of models in Systems Biology....... Advances in computing hardware and software now allow explicit simulation of some aspects of cellular dynamics close to the molecular scale. Vesicle fusion is one example of such a process. Experiments, however, typically probe cellular behavior from the molecular scale up to microns. Standard particle...... soon be coupled to Mass Action models allowing the parameters in such models to be continuously tuned according to the finer resolution simulation. This will help realize the goal of a computational cellular simulation that is able to capture the dynamics of membrane-associated processes such as...
Perturbative theory for Brownian vortexes.
Moyses, Henrique W; Bauer, Ross O; Grosberg, Alexander Y; Grier, David G
2015-06-01
Brownian vortexes are stochastic machines that use static nonconservative force fields to bias random thermal fluctuations into steadily circulating currents. The archetype for this class of systems is a colloidal sphere in an optical tweezer. Trapped near the focus of a strongly converging beam of light, the particle is displaced by random thermal kicks into the nonconservative part of the optical force field arising from radiation pressure, which then biases its diffusion. Assuming the particle remains localized within the trap, its time-averaged trajectory traces out a toroidal vortex. Unlike trivial Brownian vortexes, such as the biased Brownian pendulum, which circulate preferentially in the direction of the bias, the general Brownian vortex can change direction and even topology in response to temperature changes. Here we introduce a theory based on a perturbative expansion of the Fokker-Planck equation for weak nonconservative driving. The first-order solution takes the form of a modified Boltzmann relation and accounts for the rich phenomenology observed in experiments on micrometer-scale colloidal spheres in optical tweezers. PMID:26172698
Atomic dynamics of alumina melt: A molecular dynamics simulation study
S.Jahn
2008-03-01
Full Text Available The atomic dynamics of Al2O3 melt are studied by molecular dynamics simulation. The particle interactions are described by an advanced ionic interaction model that includes polarization effects and ionic shape deformations. The model has been shown to reproduce accurately the static structure factors S(Q from neutron and x-ray diffraction and the dynamic structure factor S(Q,ω from inelastic x-ray scattering. Analysis of the partial dynamic structure factors shows inelastic features in the spectra up to momentum transfers, Q, close to the principal peaks of partial static structure factors. The broadening of the Brillouin line widths is discussed in terms of a frequency dependent viscosity η(ω.
Engineering dynamics from the Lagrangian to simulation
Gans, Roger F
2013-01-01
This engineering dynamics textbook is aimed at beginning graduate students in mechanical engineering and other related engineering disciplines who need training in dynamics as applied to engineering mechanisms. It introduces the formal mathematical development of Lagrangian mechanics (and its corollaries), while solving numerous engineering applications. The author’s goal is to instill an understanding of the basic physics required for engineering dynamics, while providing a recipe (algorithm) for the simulation of engineering mechanisms such as robots. The book is reasonably self-contained so that the practicing engineer interested in this area can also make use of it. This book is made accessible to the widest possible audience by numerous, solved examples and diagrams that apply the principles to real engineering applications. • Provides an applied textbook for intermediate/advanced engineering dynamics courses; • Discusses Lagrangian mechanics in the context of numerous engineering applications...
秦天奇; 王飞; 杨博; 罗懋康
2015-01-01
Based on the theory of fractional integration, direct transport behaviors of coupled Brownian motors with feedback control in viscoelastic media are investigated. The mathematical model of fractional overdamped coupled Brownian motors is established by adopting the power function as damping kernel function of general Langevin equation due to the power-law memory characteristics of cytosol in biological cells. Numerical solution is observed by fractional difference method and the influence of model parameters on cooperative direct transport of the coupled Brownian motors is discussed in detail by numerical simulation. The research shows that the memory of the fractional dynamical system can affect the direct transport phenomenon of the coupled Brownian motors through changing the on-off switching frequency of the ratchet potential with feedback control. To be more specific, in a proper range of the fractional order, the memory of the dynamical system can increase the on-off switching frequency of the ratchet potential, which can lead to the velocity increase of the direct transport. Furthermore, in the case of small fractional order, since the coupled Brownian motors move under the competition between the damping force with memory and the potential force with feedback control, the resultant force exerted on the coupled particles is always positive when the ratchet potential with feedback control is on although the fractional damping force is large, which leads to the result that the coupled Brownian motors move in the positive direction in the mass. On the contrary, in the case of large fractional order, the on-off switching frequency of potential with feedback control becomes small, as a result of which the main influential factor of the direct transport becomes the potential depth. Therefore the coupled Brownian motors are more likely to stay in the potential wells for a long time because the probability that describes the possibility that the coupled Brownian
Simulation studies using multibody dynamics code DART
Keat, James E.
1989-01-01
DART is a multibody dynamics code developed by Photon Research Associates for the Air Force Astronautics Laboratory (AFAL). The code is intended primarily to simulate the dynamics of large space structures, particularly during the deployment phase of their missions. DART integrates nonlinear equations of motion numerically. The number of bodies in the system being simulated is arbitrary. The bodies' interconnection joints can have an arbitrary number of degrees of freedom between 0 and 6. Motions across the joints can be large. Provision for simulating on-board control systems is provided. Conservation of energy and momentum, when applicable, are used to evaluate DART's performance. After a brief description of DART, studies made to test the program prior to its delivery to AFAL are described. The first is a large angle reorientating of a flexible spacecraft consisting of a rigid central hub and four flexible booms. Reorientation was accomplished by a single-cycle sine wave shape torque input. In the second study, an appendage, mounted on a spacecraft, was slewed through a large angle. Four closed-loop control systems provided control of this appendage and of the spacecraft's attitude. The third study simulated the deployment of the rim of a bicycle wheel configuration large space structure. This system contained 18 bodies. An interesting and unexpected feature of the dynamics was a pulsing phenomena experienced by the stays whole playout was used to control the deployment. A short description of the current status of DART is given.
Müller, Kei W.; Meier, Christoph; Wall, Wolfgang A.
2015-12-01
Networks of crosslinked biopolymer filaments such as the cytoskeleton are the subject of intense research. Oftentimes, mechanics on the scale of single monomers (∼ 5 nm) govern the mechanics of the entire network (∼ 10 μm). Until now, one either resolved the small scales and lost the big (network) picture or focused on mechanics above the single-filament scale and neglected the molecular architecture. Therefore, the study of network mechanics influenced by the entire spectrum of relevant length scales has been infeasible so far. We propose a method that reconciles both small and large length scales without the otherwise inevitable loss in either numerical efficiency or geometrical (molecular) detail. Both explicitly modeled species, filaments and their crosslinkers, are discretized with geometrically exact beam finite elements of Simo-Reissner type. Through specific coupling conditions between the elements of the two species, mechanical joints can be established anywhere along a beam's centerline, enabling arbitrary densities of chemical binding sites. These binding sites can be oriented to model the monomeric architecture of polymers. First, we carefully discuss the method and then demonstrate its capabilities by means of a series of numerical examples.
Cleuren, Bart; Van den Broeck, Christian
2007-01-01
An asymmetric object, undergoing dissipative collisions with surrounding particles, acquires a nonzero average velocity. The latter is calculated analytically by an expansion of the Boltzmann equation and the result is compared with Monte Carlo simulations. Copyright (c) EPLA, 2007
Hydrogen Bond in Liquid Water as a Brownian Oscillator
Woutersen, Sander; Bakker, Huib J.
1999-09-01
We present the first experimental observation of a vibrational dynamic Stokes shift. This dynamic Stokes shift is observed in a femtosecond pump-probe study on the OH-stretch vibration of HDO dissolved in D2O. We find that the Stokes shift has a value of approximately 70 cm-1 and occurs with a time constant of approximately 500 femtoseconds. The measurements can be accurately described by modeling the hydrogen bond in liquid water as a Brownian oscillator.
Simulation of Boiling Water Reactor dynamics
This master thesis describes a mathematical model of a boiling water reactor and address the dynamic behaviour of the neutron kinetics, boilding dynamics and pressur stability. The simulation have been done using the SIMNON-program. The meaning were that the result from this work possibly would be adjust to supervision methods suitable for application in computer systems. This master thesis in automatic control has been done at the Department of Automatic Control, Lund Institute of Technology. The initiative to the work came from Sydkraft AB. (author)
Dynamic simulation of the mastication muscles
Weingaertner, Tim; Albrecht, Jochen
1998-05-01
The purpose of a simulated operation system in craniofacial surgery is to evaluate and visualize the results of operations on the overall facial shape of the patient and on the functionality of his jaw. This paper presents the analyzation of muscle movements in the mastication system by applying real jaw movements to the simulation. With this method an accurate modeling of the mastication muscles can be performed which is a prerequisite for a realistic simulation and precise intra- operative registration. According to this results a large- scale musculoskeletal model of the mastication system is generated including kinematic and dynamic parameters. By integrating distance sensors in the simulation of a segmented CT (computer tomograph) image of the maxilla and mandible the motions of the masticatory muscles during different kinds of jaw movements have been analyzed. The data for this motions have been recorded by a commercial system (CONDYLOCOMP LR3) on a test person and transformed to the graphical simulation system. This method for the first time allows to observe the dynamics of the mastication muscles and their different parameters like muscle length ratio and velocity. The integration of a kinematic model for the jaw movement makes it possible to analyze non traced movements.
Dynamic Simulation of a Helium Liquefier
Dynamic behavior of a helium liquefier has been studied in detail with a Cryogenic Process REal-time SimulaTor (C-PREST) at the National Institute for Fusion Science (NIFS). The C-PREST is being developed to integrate large-scale helium cryogenic plant design, operation and maintenance for optimum process establishment. As a first step of simulations of cooldown to 4.5 K with the helium liquefier model is conducted, which provides a plant-process validation platform. The helium liquefier consists of seven heat exchangers, a liquid-nitrogen (LN2) precooler, two expansion turbines and a liquid-helium (LHe) reservoir. Process simulations are fulfilled with sequence programs, which were implemented with C-PREST based on an existing liquefier operation. The interactions of a JT valve, a JT-bypass valve and a reservoir-return valve have been dynamically simulated. The paper discusses various aspects of refrigeration process simulation, including its difficulties such as a balance between complexity of the adopted models and CPU time
Parallelization of quantum molecular dynamics simulation code
A quantum molecular dynamics simulation code has been developed for the analysis of the thermalization of photon energies in the molecule or materials in Kansai Research Establishment. The simulation code is parallelized for both Scalar massively parallel computer (Intel Paragon XP/S75) and Vector parallel computer (Fujitsu VPP300/12). Scalable speed-up has been obtained with a distribution to processor units by division of particle group in both parallel computers. As a result of distribution to processor units not only by particle group but also by the particles calculation that is constructed with fine calculations, highly parallelization performance is achieved in Intel Paragon XP/S75. (author)
Description of the grout system dynamic simulation
The grout system dynamic computer simulation was created to allow investigation of the ability of the grouting system to meet established milestones, for various assumed system configurations and parameters. The simulation simulates the movement of tank waste through the system versus time, from initial storage tanks, through feed tanks and the grout plant, then finally to a grout vault. The simulation properly accounts for the following (1) time required to perform various actions or processes, (2) delays involved in gaining regulatory approval, (3) random system component failures, (4) limitations on equipment capacities, (5) available parallel components, and (6) different possible strategies for vault filling. The user is allowed to set a variety of system parameters for each simulation run. Currently, the output of a run primarily consists of a plot of projected grouting campaigns completed versus time, for comparison with milestones. Other outputs involving any model component can also be quickly created or deleted as desired. In particular, sensitivity runs where the effect of varying a model parameter (flow rates, delay times, number of feed tanks available, etc.) on the ability of the system to meet milestones can be made easily. The grout system simulation was implemented using the ITHINK* simulation language for Macintosh** computers
Molecular dynamics simulation for modelling plasma spectroscopy
Talin, B; Calisti, A; Gigosos, M A; González, M A; Gaztelurrutia, T R; Dufty, J W
2003-01-01
The ion-electron coupling properties for an ion impurity in an electron gas and for a two-component plasma are carried out on the basis of a regularized electron-ion potential removing the short-range Coulomb divergence. This work is largely motivated by the study of radiator dipole relaxation in plasmas which makes a real link between models and experiments. Current radiative property models for plasmas include single electron collisions neglecting charge-charge correlations within the classical quasi-particle approach commonly used in this field. The dipole relaxation simulation based on electron-ion molecular dynamics proposed here will provide a means to benchmark and improve model developments. Benefiting from a detailed study of a single ion embedded in an electron plasma, the challenging two-component ion-electron molecular dynamics simulations are proved accurate. They open new possibilities of obtaining reference lineshape data.
Simulation of Abrasive Machining Using Molecular Dynamics
Oluwajobi, Akinjide O.; Chen, Xun
2009-01-01
The development of ultra–precision processes which can achieve excellent surface finish and tolerance at nanometre level is now a critical requirement for many applications in medical, electronics and energy industry. Presently, it is very difficult to observe the diverse microscopic physical phenomena occurring in nanometric machining through experiments. The use of Molecular Dynamics (MD) simulation has proved to be an effective tool for the prediction and the analysis ...
Parallel beam dynamics simulation of linear accelerators
Qiang, Ji; Ryne, Robert D.
2002-01-01
In this paper we describe parallel particle-in-cell methods for the large scale simulation of beam dynamics in linear accelerators. These techniques have been implemented in the IMPACT (Integrated Map and Particle Accelerator Tracking) code. IMPACT is being used to study the behavior of intense charged particle beams and as a tool for the design of next-generation linear accelerators. As examples, we present applications of the code to the study of emittance exchange in high intensity b...
Molecular dynamics simulation of expanding infinite matter
Multi-fragmentation occurred in an expanding infinite system is studied by using molecular dynamics simulation. To evaluate the secondary decay effect, the time evolution of expanding system is proceeded till all fragments are stabilized completely. The fragment mass distribution from the expansion is compared with a percolation model and the cause of the exponential shape is clarified. The cause of small critical temperature is also discussed. (author)
Schwinger model simulations with dynamical overlap fermions
Bietenholz, Wolfgang; Shcheredin, Stanislav; Volkholz, Jan
2007-01-01
We present simulation results for the 2-flavour Schwinger model with dynamical overlap fermions. In particular we apply the overlap hypercube operator at seven light fermion masses. In each case we collect sizable statistics in the topological sectors 0 and 1. Since the chiral condensate Sigma vanishes in the chiral limit, we observe densities for the microscopic Dirac spectrum, which have not been addressed yet by Random Matrix Theory (RMT). Nevertheless, by confronting the averages of the l...
Temperature Dependent Molecular Dynamic Simulation of Friction
Dias, R A; Coura, P Z; Costa, B V
2006-01-01
In this work we present a molecular dynamics simulation of a FFM experiment. The tip-sample interaction is studied by varying the normal force in the tip and the temperature of the surface. The friction force, cA, at zero load and the friction coefficient, $\\mu$, were obtained. Our results strongly support the idea that the effective contact area, A, decreases with increasing temperature and the friction coefficient presents a clear signature of the premelting process of the surface.
Dislocation dynamics simulations in a cylinder
Weinberger, Christopher R; Aubry, Sylvie; Cai, Wei [Department of Mechanical Engineering, Stanford University, Stanford, CA 94305-4040 (United States); Lee, Seok-Woo, E-mail: cweinber@stanford.edu, E-mail: caiwei@stanford.edu [Department of Materials Science and Engineering, Stanford University, Stanford, CA 94305-4034 (United States)
2009-07-15
In this work we describe how to perform dislocation dynamics simulations in a cylindrical geometry. An algorithm for computing the image stress is given in detail including methods for handling the singularity. Additional remesh rules address the problems of the cylindrical geometry and the required self consistency with mobility laws. Numerical studies benchmark the accuracy of the algorithms and the importance of handling the singularity correctly.
Simulating granular flow with molecular dynamics
Ristow, Gerald
1992-01-01
We investigate by means of Molecular Dynamics simulations an assembly of spheres to model a granular medium flowing from an upper rectangular chamber through a hole into a lower chamber. Two different two dimensional models are discussed one of them including rotations of the individual spheres. The outflow properties are investigated and compared to experimental data. The qualitative agreement suggests that our models contain the necessary ingredients to describe the outflow properties of gr...
Simulation of counterflow pedestrian dynamics using spheropolygons
Alonso-Marroquín, Fernando; Busch, Jonathan; Chiew, Coraline; Lozano, Celia; Ramírez-Gómez, Álvaro
2014-12-01
Pedestrian dynamic models are typically designed for comfortable walking or slightly congested conditions and typically use a single disk or combination of three disks for the shape of a pedestrian. Under crowd conditions, a more accurate pedestrian shape has advantages over the traditional single or three-disks model. We developed a method for simulating pedestrian dynamics in a large dense crowd of spheropolygons adapted to the cross section of the chest and arms of a pedestrian. Our numerical model calculates pedestrian motion from Newton's second law, taking into account viscoelastic contact forces, contact friction, and ground-reaction forces. Ground-reaction torque was taken to arise solely from the pedestrians' orientation toward their preferred destination. Simulations of counterflow pedestrians dynamics in corridors were used to gain insight into a tragic incident at the Madrid Arena pavilion in Spain, where five girls were crushed to death. The incident took place at a Halloween Celebration in 2012, in a long, densely crowded hallway used as entrance and exit at the same time. Our simulations reconstruct the mechanism of clogging in the hallway. The hypothetical case of a total evacuation order was also investigated. The results highlights the importance of the pedestrians' density and the effect of counterflow in the onset of avalanches and clogging and provides an estimation of the number of injuries based on a calculation of the contact-force network between the pedestrians.
Electronic continuum model for molecular dynamics simulations.
Leontyev, I V; Stuchebrukhov, A A
2009-02-28
A simple model for accounting for electronic polarization in molecular dynamics (MD) simulations is discussed. In this model, called molecular dynamics electronic continuum (MDEC), the electronic polarization is treated explicitly in terms of the electronic continuum (EC) approximation, while the nuclear dynamics is described with a fixed-charge force field. In such a force-field all atomic charges are scaled to reflect the screening effect by the electronic continuum. The MDEC model is rather similar but not equivalent to the standard nonpolarizable force-fields; the differences are discussed. Of our particular interest is the calculation of the electrostatic part of solvation energy using standard nonpolarizable MD simulations. In a low-dielectric environment, such as protein, the standard MD approach produces qualitatively wrong results. The difficulty is in mistreatment of the electronic polarizability. We show how the results can be much improved using the MDEC approach. We also show how the dielectric constant of the medium obtained in a MD simulation with nonpolarizable force-field is related to the static (total) dielectric constant, which includes both the nuclear and electronic relaxation effects. Using the MDEC model, we discuss recent calculations of dielectric constants of alcohols and alkanes, and show that the MDEC results are comparable with those obtained with the polarizable Drude oscillator model. The applicability of the method to calculations of dielectric properties of proteins is discussed. PMID:19256627
Generator dynamics in aeroelastic analysis and simulations
Larsen, Torben J.; Hansen, Morten Hartvig; Iov, F.
2003-01-01
of the first order terms in the model as well as the influence on drive train eigenfrequencies and damping has been investigated. Load response during timesimulation of wind turbine response have been compared to simulations with a traditional static generator model based entirely on the slip angle....... A 2 MW turbine has been modelled in the aeroelastic code HAWC. When using the new dynamic generator modelthere is an interesting coupling between the generator dynamics and a global turbine vibration mode at 4.5 Hz, which only occurs when a dynamic formulation of the generator equations is applied....... This frequency can especially be seen in the electricalpower of the generator and the rotational speed of the generator, but also as torque variations in the drive train....
Adaptive wavelet simulation of global ocean dynamics
N. K.-R. Kevlahan
2015-07-01
Full Text Available In order to easily enforce solid-wall boundary conditions in the presence of complex coastlines, we propose a new mass and energy conserving Brinkman penalization for the rotating shallow water equations. This penalization does not lead to higher wave speeds in the solid region. The error estimates for the penalization are derived analytically and verified numerically for linearized one dimensional equations. The penalization is implemented in a conservative dynamically adaptive wavelet method for the rotating shallow water equations on the sphere with bathymetry and coastline data from NOAA's ETOPO1 database. This code could form the dynamical core for a future global ocean model. The potential of the dynamically adaptive ocean model is illustrated by using it to simulate the 2004 Indonesian tsunami and wind-driven gyres.
Dynamic simulation analyzes expanded refinery steam system
Gandhi, S.L.; Graham, J.; Duffield, M.A.; Cortes, R.M. [M.W. Kellogg Co., Houston, TX (United States)
1995-11-01
There is a rising demand for more-stable process operations with increasing facility modernization and automation. This requires steam system controls that provide effective and stable responses to operational disturbances. Today it is possible to use dynamic simulation analysis for designing steam systems. The main objectives were to: develop a high-fidelity dynamic simulation model of the steam system and its associated controls; develop dynamic system responses to major plant disturbances; and assist in designing new control strategies and modifying existing controls to minimize the impact of plant upsets on system performance. The steam system for this project included two power plant supplying steam and power and a steam pipe network comprising three levels of steam with stream flows to and from process units spread over several square kilometers. Whether the steam system and its controls would respond effectively to major operational upsets, such as boiler, steam turbogenerator or process gas compressor trips, or a steep increase in steam demand, can only be determined from evaluating dynamic responses to these incidents. The model consisted of about 5,000 equations, 1,000 parameters and data items, and was used for evaluating steam system responses to various operational upset scenarios. The responses were used to identify limitations of the as designed control system. The model also proved to be an effective tool for verifying adequacy of the proposed recommendations. These included modifications to control valve letdown stations, steam turbogenerator governor controls and the boiler master-pressure controller.
SIMULATION OF INTERLINE DYNAMIC VOLTAGE RESTORER
J.Singaravelan
2011-08-01
Full Text Available This paper presents a new approach for the dynamic control of a current source inverter (CSI using Super Conductive Magnetic energy storage (SMES based Interline DVR. The dynamic voltage restorer (DVR provides a technically advanced and economical solution to voltage-sag problem. As the voltage-restoration process involves the real-power injection into the distribution system, the capability ofa DVR, especially for compensating long-duration voltage sags, it depends on the energy storage capacity of the DVR. The interline DVR proposed in this paper provides a way to replenish Dc-link energy storage dynamically. The IDVR consists of several DVRs connected to different distribution feeders in the power system. The DVRs in the IDVR system shares the common energy storage. When one of the DVRcompensates for voltage sag appearing in that feeder, the other DVRs replenish the energy in the common dc-link dynamically. Thus, one DVR in the IDVR system works in voltage-sag compensation mode whilethe other DVRs in the IDVR system operate in power-flow control mode. The proposed topology is simulated using Matlab/Simulink and total IDVR system is simulated using Matlab/Simulink.
Dynamic Shear Modulus of Polymers from Molecular Dynamics Simulations
Byutner, Oleksiy; Smith, Grant
2001-03-01
In this work we describe the methodology for using equilibrium molecular dynamics simulations (MD) simulations to obtain the viscoelastic properties of polymers in the glassy regime. Specifically we show how the time dependent shear stress modulus and frequency dependent complex shear modulus in the high-frequency regime can be determined from the off-diagonal terms of the stress-tensor autocorrelation function obtained from MD trajectories using the Green-Kubo method and appropriate Fourier transforms. In order to test the methodology we have performed MD simulations of a low-molecular-weight polybutadiene system using quantum chemistry based potential functions. Values of the glassy modulus and the maximum loss frequency were found to be in good agreement with experimental data for polybutadiene at 298 K.
A molecular dynamics simulation code ISIS
Computer simulation based on the molecular dynamics (MD) method has become an important tool complementary to experiments and theoretical calculations in a wide range of scientific fields such as physics, chemistry, biology, and so on. In the MD method, the Newtonian equations-of-motion of classical particles are integrated numerically to reproduce a phase-space trajectory of the system. In the 1980's, several new techniques have been developed for simulation at constant-temperature and/or constant-pressure in convenient to compare result of computer simulation with experimental results. We first summarize the MD method for both microcanonical and canonical simulations. Then, we present and overview of a newly developed ISIS (Isokinetic Simulation of Soft-spheres) code and its performance on various computers including vector processors. The ISIS code has a capability to make a MD simulation under constant-temperature condition by using the isokinetic constraint method. The equations-of-motion is integrated by a very accurate fifth-order finite differential algorithm. The bookkeeping method is also utilized to reduce the computational time. Furthermore, the ISIS code is well adopted for vector processing: Speedup ratio ranged from 16 to 24 times is obtained on a VP2600/10 vector processor. (author)
Dynamic simulator for PEFC propulsion plant
Hiraide, Masataka; Kaneda, Eiichi; Sato, Takao [Mitsui Engineering & Shipbuilding Co., Ltd., Tokyo (Japan)] [and others
1996-12-31
This report covers part of a joint study on a PEFC propulsion system for surface ships, summarized in a presentation to this Seminar, entitled {open_quote}Study on a Polymer Electrolyte Fuel Cell (PEFC) Propulsion System for Surface Ships{close_quotes}, and which envisages application to a 1,500 DWT cargo vessel. The work presented here focuses on a simulation study on PEFC propulsion plant performance, and particularly on the system response to changes in load. Using a dynamic simulator composed of system components including fuel cell, various simulations were executed, to examine the performance of the system as a whole and of the individual system components under quick and large load changes such as occasioned by maneuvering operations and by racing when the propeller emerges above water in heavy sea.
Efficiency of Brownian heat engines.
Derényi, I; Astumian, R D
1999-06-01
We study the efficiency of one-dimensional thermally driven Brownian ratchets or heat engines. We identify and compare the three basic setups characterized by the type of the connection between the Brownian particle and the two heat reservoirs: (i) simultaneous, (ii) alternating in time, and (iii) position dependent. We make a clear distinction between the heat flow via the kinetic and the potential energy of the particle, and show that the former is always irreversible and it is only the third setup where the latter is reversible when the engine works quasistatically. We also show that in the third setup the heat flow via the kinetic energy can be reduced arbitrarily, proving that even for microscopic heat engines there is no fundamental limit of the efficiency lower than that of a Carnot cycle. PMID:11969723
Brownian motion and stochastic calculus
Karatzas, Ioannis
1998-01-01
This book is designed as a text for graduate courses in stochastic processes. It is written for readers familiar with measure-theoretic probability and discrete-time processes who wish to explore stochastic processes in continuous time. The vehicle chosen for this exposition is Brownian motion, which is presented as the canonical example of both a martingale and a Markov process with continuous paths. In this context, the theory of stochastic integration and stochastic calculus is developed. The power of this calculus is illustrated by results concerning representations of martingales and change of measure on Wiener space, and these in turn permit a presentation of recent advances in financial economics (option pricing and consumption/investment optimization). This book contains a detailed discussion of weak and strong solutions of stochastic differential equations and a study of local time for semimartingales, with special emphasis on the theory of Brownian local time. The text is complemented by a large num...
Hänggi, Peter; Marchesoni, Fabio
2005-01-01
In the year 1905 Albert Einstein published four papers that raised him to a giant in the history of science of all times. These works encompass the photon hypothesis (for which he obtained the Nobel prize in 1921), his first two papers on (special) relativity theory and, of course, his first paper on Brownian motion, entitled "\\"Uber die von der molekularkinetischen Theorie der W\\"arme geforderte Bewegung von in ruhenden Fl\\"ussigkeiten suspendierten Teilchen'' (submitted on May 11, 1905). Th...
Kingman's coalescent and Brownian motion
Berestycki, J.; Berestycki, N
2009-01-01
We describe a simple construction of Kingman's coalescent in terms of a Brownian excursion. This construction is closely related to, and sheds some new light on, earlier work by Aldous and Warren. Our approach also yields some new results: for instance, we obtain the full multifractal spectrum of Kingman's coalescent. This complements earlier work on Beta-coalescents by the authors and Schweinsberg. Surprisingly, the thick part of the spectrum is not obtained by taking the limit as $\\alpha \\t...
Brownian ratchets and Parrondo's games
Harmer, Gregory P.; Abbott, Derek; Taylor, Peter G.; Parrondo, Juan M. R.
2001-09-01
Parrondo's games present an apparently paradoxical situation where individually losing games can be combined to win. In this article we analyze the case of two coin tossing games. Game B is played with two biased coins and has state-dependent rules based on the player's current capital. Game B can exhibit detailed balance or even negative drift (i.e., loss), depending on the chosen parameters. Game A is played with a single biased coin that produces a loss or negative drift in capital. However, a winning expectation is achieved by randomly mixing A and B. One possible interpretation pictures game A as a source of "noise" that is rectified by game B to produce overall positive drift—as in a Brownian ratchet. Game B has a state-dependent rule that favors a losing coin, but when this state dependence is broken up by the noise introduced by game A, a winning coin is favored. In this article we find the parameter space in which the paradoxical effect occurs and carry out a winning rate analysis. The significance of Parrondo's games is that they are physically motivated and were originally derived by considering a Brownian ratchet—the combination of the games can be therefore considered as a discrete-time Brownian ratchet. We postulate the use of games of this type as a toy model for a number of physical and biological processes and raise a number of open questions for future research.