-
Driving Spin-Boson Models From Equilibrium Using Exact Quantum Dynamics
Authors:
G. M. G. McCaul,
C. D. Lorenz,
L. Kantorovich
Abstract:
We present an application of the Extended Stochastic Liouville Equation (ESLE) Phys. Rev. B 95, 125124, which gives an exact solution for the reduced density matrix of an open system surrounded by a harmonic heat bath. This method considers the extended system (the open system and the bath) being thermally equilibrated prior to the action of a time dependent perturbation, as opposed to the usual a…
▽ More
We present an application of the Extended Stochastic Liouville Equation (ESLE) Phys. Rev. B 95, 125124, which gives an exact solution for the reduced density matrix of an open system surrounded by a harmonic heat bath. This method considers the extended system (the open system and the bath) being thermally equilibrated prior to the action of a time dependent perturbation, as opposed to the usual assumption that system and bath are initially partitioned. This is an exact technique capable of accounting for arbitrary parameter regimes of the model. Here we present our first numerical implementation of the method in the simplest case of a Caldeira-Leggett representation of the bath Hamiltonian, and apply it to a spin-Boson system driven from coupled equilibrium. We observe significant behaviours in both the transient dynamics and asymptotic states of the reduced density matrix not present in the usual approximation.
△ Less
Submitted 18 December, 2017;
originally announced December 2017.
-
Modelling a Bistable System Strongly Coupled to a Debye Bath: A Quasiclassical Approach Based on the Generalised Langevin Equation
Authors:
L. Stella,
H. Ness,
C. D. Lorenz,
L. Kantorovich
Abstract:
Bistable systems present two degenerate metastable configurations separated by an energy barrier. Thermal or quantum fluctuations can promote the transition between the configurations at a rate which depends on the dynamical properties of the local environment (i.e., a thermal bath). In the case of classical systems, strong system-bath interaction has been successfully modelled by the Generalised…
▽ More
Bistable systems present two degenerate metastable configurations separated by an energy barrier. Thermal or quantum fluctuations can promote the transition between the configurations at a rate which depends on the dynamical properties of the local environment (i.e., a thermal bath). In the case of classical systems, strong system-bath interaction has been successfully modelled by the Generalised Langevin Equation (GLE) formalism. Here we show that the efficient GLE algorithm introduced in Phys. Rev. B 89, 134303 (2014) can be extended to include some crucial aspects of the quantum fluctuations. In particular, the expected isotopic effect is observed along with the convergence of the quantum and classical transition rates in the strong coupling limit. Saturation of the transition rates at low temperature is also retrieved, in qualitative, yet not quantitative, agreement with the analytic predictions. The discrepancies in the tunnelling regime are due to an incorrect sampling close to the barrier top. The domain of applicability of the quasiclassical GLE is also discussed.
△ Less
Submitted 12 April, 2017;
originally announced April 2017.
-
Partition-free approach to open quantum systems in harmonic environments: an exact stochastic Liouville equation
Authors:
G. M. G. McCaul,
C. D. Lorenz,
L. Kantorovich
Abstract:
We present a partition-free approach to the evolution of density matrices for open quantum systems coupled to a harmonic environment. The influence functional formalism combined with a two-time Hubbard-Stratonovich transformation allows us to derive a set of exact differential equations for the reduced density matrix of an open system, termed the Extended Stochastic Liouville-von Neumann equation.…
▽ More
We present a partition-free approach to the evolution of density matrices for open quantum systems coupled to a harmonic environment. The influence functional formalism combined with a two-time Hubbard-Stratonovich transformation allows us to derive a set of exact differential equations for the reduced density matrix of an open system, termed the Extended Stochastic Liouville-von Neumann equation. Our approach generalises previous work based on Caldeira-Leggett models and a partitioned initial density matrix. This provides a simple, yet exact, closed-form description for the evolution of open systems from equilibriated initial conditions. The applicability of this model and the potential for numerical implementations are also discussed.
△ Less
Submitted 3 January, 2017; v1 submitted 16 December, 2016;
originally announced December 2016.
-
Nonequilibrium Generalised Langevin Equation for the calculation of heat transport properties in model 1D atomic chains coupled to two 3D thermal baths
Authors:
H. Ness,
L. Stella,
C. D. Lorenz,
L. Kantorovich
Abstract:
We use a Generalised Langevin Equation (GLE) scheme to study the thermal transport of low dimensional systems. In this approach, the central classical region is connected to two realistic thermal baths kept at two different temperatures [H. Ness et al., Phys. Rev. B {\bf 93}, 174303 (2016)]. We consider model Al systems, i.e. one-dimensional atomic chains connected to three-dimensional baths. The…
▽ More
We use a Generalised Langevin Equation (GLE) scheme to study the thermal transport of low dimensional systems. In this approach, the central classical region is connected to two realistic thermal baths kept at two different temperatures [H. Ness et al., Phys. Rev. B {\bf 93}, 174303 (2016)]. We consider model Al systems, i.e. one-dimensional atomic chains connected to three-dimensional baths. The thermal transport properties are studied as a function of the chain length $N$ and the temperature difference $ΔT$ between the baths. We calculate the transport properties both in the linear response regime and in the non-linear regime. Two different laws are obtained for the linear conductance versus the length of the chains. For large temperatures ($T \gtrsim 500$ K) and temperature differences ($ΔT \gtrsim 500$ K), the chains, with $N > 18$ atoms, present a diffusive transport regime with the presence of a temperature gradient across the system. For lower temperatures($T \lesssim 500$ K) and temperature differences ($ΔT \lesssim 400$ K), a regime similar to the ballistic regime is observed. Such a ballistic-like regime is also obtained for shorter chains ($N \le 15 $). Our detailed analysis suggests that the behaviour at higher temperatures and temperature differences is mainly due to anharmonic effects within the long chains.
△ Less
Submitted 8 April, 2017; v1 submitted 3 December, 2016;
originally announced December 2016.
-
Nonequilibrium processes from Generalised Langevin Equations: realistic nanoscale systems connected to two thermal baths
Authors:
H. Ness,
A. Genina,
L. Stella,
C. D. Lorenz,
L. Kantorovich
Abstract:
We extend the Generalised Langevin Equation (GLE) method [Phys. Rev. B 89, 134303 (2014)] to model a central classical region connected to two realistic thermal baths at two different temperatures. In such nonequilibrium conditions a heat flow is established, via the central system, in between the two baths. The GLE-2B (GLE two baths) scheme permits us to have a realistic description of both the d…
▽ More
We extend the Generalised Langevin Equation (GLE) method [Phys. Rev. B 89, 134303 (2014)] to model a central classical region connected to two realistic thermal baths at two different temperatures. In such nonequilibrium conditions a heat flow is established, via the central system, in between the two baths. The GLE-2B (GLE two baths) scheme permits us to have a realistic description of both the dissipative central system and its surrounding baths. Following the original GLE approach, the extended Langevin dynamics scheme is modified to take into account two sets of auxiliary degrees of freedom corresponding to the mapping of the vibrational properties of each bath. These auxiliary variables are then used to solve the non-Markovian dissipative dynamics of the central region. The resulting algorithm is used to study a model of a short Al nanowire connected to two baths. The results of the simulations using the GLE-2B approach are compared to the results of other simulations that were carried out using standard thermostatting approaches (based on Markovian Langevin and Nose-Hoover thermostats). We concentrate on the steady state regime and study the establishment of a local temperature profile within the system. The conditions for obtaining a flat profile or a temperature gradient are examined in detail, in agreement with earlier studies. The results show that the GLE-2B approach is able to treat, within a single scheme, two widely different thermal transport regimes, i.e. ballistic systems, with no temperature gradient, and diffusive systems with a temperature gradient.
△ Less
Submitted 30 April, 2016; v1 submitted 11 December, 2015;
originally announced December 2015.
-
Applications of the Generalised Langevin Equation: towards a realistic description of the baths
Authors:
H. Ness,
L. Stella,
C. D. Lorenz,
L. Kantorovich
Abstract:
The Generalised Langevin Equation (GLE) method, as developed in Ref. [Phys. Rev. B 89, 134303 (2014)], is used to calculate the dissipative dynamics of systems described at the atomic level. The GLE scheme goes beyond the commonly used bilinear coupling between the central system and the bath, and permits us to have a realistic description of both the dissipative central system and its surrounding…
▽ More
The Generalised Langevin Equation (GLE) method, as developed in Ref. [Phys. Rev. B 89, 134303 (2014)], is used to calculate the dissipative dynamics of systems described at the atomic level. The GLE scheme goes beyond the commonly used bilinear coupling between the central system and the bath, and permits us to have a realistic description of both the dissipative central system and its surrounding bath. We show how to obtain the vibrational properties of a realistic bath and how to convey such properties into an extended Langevin dynamics by the use of the mapping of the bath vibrational properties onto a set of auxiliary variables. Our calculations for a model of a Lennard-Jones solid show that our GLE scheme provides a stable dynamics, with the dissipative/relaxation processes properly described. The total kinetic energy of the central system always thermalises toward the expected bath temperature, with appropriate fluctuation around the mean value. More importantly, we obtain a velocity distribution for the individual atoms in the central system which follows the expected canonical distribution at the corresponding temperature. This confirms that both our GLE scheme and our mapping procedure onto an extended Langevin dynamics provide the correct thermostat. We also examined the velocity autocorrelation functions and compare our results with more conventional Langevin dynamics.
△ Less
Submitted 18 December, 2014;
originally announced December 2014.
-
The Generalized Langevin Equation: An efficient approach to non-equilibrium molecular dynamics of open systems
Authors:
L. Stella,
C. D. Lorenz,
L. Kantorovich
Abstract:
The Generalized Langevin Equation (GLE) has been recently suggested to simulate the time evolution of classical solid and molecular systems when considering general non-equilibrium processes. In this approach, a part of the whole system (an open system), which interacts and exchanges energy with its dissipative environment, is studied. Because the GLE is derived by projecting out exactly the harmo…
▽ More
The Generalized Langevin Equation (GLE) has been recently suggested to simulate the time evolution of classical solid and molecular systems when considering general non-equilibrium processes. In this approach, a part of the whole system (an open system), which interacts and exchanges energy with its dissipative environment, is studied. Because the GLE is derived by projecting out exactly the harmonic environment, the coupling to it is realistic, while the equations of motion are non-Markovian. Although the GLE formalism has already found promising applications, e.g., in nanotribology and as a powerful thermostat for equilibration in classical molecular dynamics simulations, efficient algorithms to solve the GLE for realistic memory kernels are highly non-trivial, especially if the memory kernels decay non-exponentially. This is due to the fact that one has to generate a colored noise and take account of the memory effects in a consistent manner. In this contribution, we present a simple, yet efficient, algorithm for solving the GLE for practical memory kernels and we demonstrate its capability for the exactly solvable case of a harmonic oscillator coupled to a Debye bath.
△ Less
Submitted 18 March, 2014; v1 submitted 17 December, 2013;
originally announced December 2013.
-
On the estimation of the curvatures and bending rigidity of membrane networks via a local maximum-entropy approach
Authors:
Fernando Fraternali,
Christan D. Lorenz,
Gianluca Marcelli
Abstract:
We present a meshfree method for the curvature estimation of membrane networks based on the Local Maximum Entropy approach recently presented in (Arroyo and Ortiz, 2006). A continuum regularization of the network is carried out by balancing the maximization of the information entropy corresponding to the nodal data, with the minimization of the total width of the shape functions. The accuracy and…
▽ More
We present a meshfree method for the curvature estimation of membrane networks based on the Local Maximum Entropy approach recently presented in (Arroyo and Ortiz, 2006). A continuum regularization of the network is carried out by balancing the maximization of the information entropy corresponding to the nodal data, with the minimization of the total width of the shape functions. The accuracy and convergence properties of the given curvature prediction procedure are assessed through numerical applications to benchmark problems, which include coarse grained molecular dynamics simulations of the fluctuations of red blood cell membranes (Marcelli et al., 2005; Hale et al., 2009). We also provide an energetic discrete-to-continuum approach to the prediction of the zero-temperature bending rigidity of membrane networks, which is based on the integration of the local curvature estimates. The Local Maximum Entropy approach is easily applicable to the continuum regularization of fluctuating membranes, and the prediction of membrane and bending elasticities of molecular dynamics models.
△ Less
Submitted 16 May, 2011;
originally announced May 2011.
-
Forces between functionalized silica nanoparticles in solution
Authors:
J. Matthew D. Lane,
Ahmed E. Ismail,
Michael Chandross,
Christian D. Lorenz,
Gary S. Grest
Abstract:
To prevent the flocculation and phase separation of nanoparticles in solution, nanoparticles are often functionalized with short chain surfactants. Here we present fully-atomistic molecular dynamics simulations which characterize how these functional coatings affect the interactions between nanoparticles and with the surrounding solvent. For 5 nm diameter silica nanoparticles coated with poly(et…
▽ More
To prevent the flocculation and phase separation of nanoparticles in solution, nanoparticles are often functionalized with short chain surfactants. Here we present fully-atomistic molecular dynamics simulations which characterize how these functional coatings affect the interactions between nanoparticles and with the surrounding solvent. For 5 nm diameter silica nanoparticles coated with poly(ethylene oxide) (PEO) oligomers in water, we determined the hydrodynamic drag on two approaching nanoparticles moving through solvent and on a single nanoparticle as it approaches a planar surface. In most circumstances, acroscale fluid theory accurately predicts the drag on these nano-scale particles. Good agreement is seen with Brenner's analytical solutions for wall separations larger than the soft nanoparticle radius. For two approaching coated nanoparticles, the solvent-mediated (velocity-independent) and lubrication (velocity-dependent) forces are purely repulsive and do not exhibit force oscillations that are typical of uncoated rigid spheres.
△ Less
Submitted 16 April, 2009;
originally announced April 2009.
-
Micellar Crystals in Solution from Molecular Dynamics Simulations
Authors:
J. A. Anderson,
C. D. Lorenz,
A. Travesset
Abstract:
Polymers with both soluble and insoluble blocks typically self-assemble into micelles, aggregates of a finite number of polymers where the soluble blocks shield the insoluble ones from contact with the solvent. Upon increasing concentration, these micelles often form gels that exhibit crystalline order in many systems. In this paper, we present a study of both the dynamics and the equilibrium pr…
▽ More
Polymers with both soluble and insoluble blocks typically self-assemble into micelles, aggregates of a finite number of polymers where the soluble blocks shield the insoluble ones from contact with the solvent. Upon increasing concentration, these micelles often form gels that exhibit crystalline order in many systems. In this paper, we present a study of both the dynamics and the equilibrium properties of micellar crystals of triblock polymers using molecular dynamics simulations. Our results show that equilibration of single micelle degrees of freedom and crystal formation occurs by polymer transfer between micelles, a process that is described by transition state theory. Near the disorder (or melting) transition, bcc lattices are favored for all triblocks studied. Lattices with fcc ordering are also found, but only at lower kinetic temperatures and for triblocks with short hydrophilic blocks. Our results lead to a number of theoretical considerations and suggest a range of implications to experimental systems with a particular emphasis on Pluronic polymers.
△ Less
Submitted 18 January, 2008;
originally announced January 2008.
-
Charge Inversion of Divalent Ionic Solutions in Silica Channels
Authors:
Christian D. Lorenz,
Alex Travesset
Abstract:
Recent experiments (F.H.J. Van Der Heyden et al., PRL 96, 224502 (2006)) of streaming currents in silica nanochannels with divalent ions report charge inversion, i.e. interfacial charges attracting counterions in excess of their own nominal charge, in conflict with existing theoretical and simulation results. We reveal the mechanism of charge inversion by using all-atomic molecular dynamics simu…
▽ More
Recent experiments (F.H.J. Van Der Heyden et al., PRL 96, 224502 (2006)) of streaming currents in silica nanochannels with divalent ions report charge inversion, i.e. interfacial charges attracting counterions in excess of their own nominal charge, in conflict with existing theoretical and simulation results. We reveal the mechanism of charge inversion by using all-atomic molecular dynamics simulations. Our results show excellent agreement with experiments, both qualitatively and quantitatively. We further discuss the implications of our study for the general problem of ionic correlations in solutions as well as in regards of the properties of silica-water interfaces.
△ Less
Submitted 24 June, 2007;
originally announced June 2007.
-
Excess number of percolation clusters on the surface of a sphere
Authors:
Christian D. Lorenz,
Robert M. Ziff
Abstract:
Monte Carlo simulations were performed in order to determine the excess number of clusters b and the average density of clusters n_c for the two-dimensional "Swiss cheese" continuum percolation model on a planar L x L system and on the surface of a sphere. The excess number of clusters for the L x L system was confirmed to be a universal quantity with a value b = 0.8841 as previously predicted a…
▽ More
Monte Carlo simulations were performed in order to determine the excess number of clusters b and the average density of clusters n_c for the two-dimensional "Swiss cheese" continuum percolation model on a planar L x L system and on the surface of a sphere. The excess number of clusters for the L x L system was confirmed to be a universal quantity with a value b = 0.8841 as previously predicted and verified only for lattice percolation. The excess number of clusters on the surface of a sphere was found to have the value b = 1.215(1) for discs with the same coverage as the flat critical system. Finally, the average critical density of clusters was calculated for continuum systems n_c = 0.0408(1).
△ Less
Submitted 1 February, 2001;
originally announced February 2001.
-
Similarity of percolation thresholds on the hcp and fcc lattices
Authors:
Christian D. Lorenz,
Raechelle May,
Robert M. Ziff
Abstract:
Extensive Monte-Carlo simulations were performed in order to determine the precise values of the critical thresholds for site ($p^{hcp}_{c,S} = 0.199 255 5 \pm 0.000 001 0$) and bond ($p^{hcp}_{c,B} = 0.120 164 0 \pm 0.000 001 0$) percolation on the hcp lattice to compare with previous precise measuremens on the fcc lattice. Also, exact enumeration of the hcp and fcc lattices was performed and y…
▽ More
Extensive Monte-Carlo simulations were performed in order to determine the precise values of the critical thresholds for site ($p^{hcp}_{c,S} = 0.199 255 5 \pm 0.000 001 0$) and bond ($p^{hcp}_{c,B} = 0.120 164 0 \pm 0.000 001 0$) percolation on the hcp lattice to compare with previous precise measuremens on the fcc lattice. Also, exact enumeration of the hcp and fcc lattices was performed and yielded generating functions and series for the zeroth, first, and second moments of both lattices. When these series and the values of $p_c$ are compared to those for the fcc lattice, it is apparent that the site percolation thresholds are different; however, the bond percolation thresholds are equal within error bars, and the series only differ slightly in the higher order terms, suggesting the actual values are very close to each other, if not identical.
△ Less
Submitted 2 August, 1999;
originally announced August 1999.
-
Shape-dependent universality in percolation
Authors:
Robert M. Ziff,
Christian D. Lorenz,
Peter Kleban
Abstract:
The shape-dependent universality of the excess percolation cluster number and cross-configuration probability on a torus is discussed. Besides the aspect ratio of the torus, the universality class depends upon the twist in the periodic boundary conditions, which for example are generally introduced when triangular lattices are used in simulations.
The shape-dependent universality of the excess percolation cluster number and cross-configuration probability on a torus is discussed. Besides the aspect ratio of the torus, the universality class depends upon the twist in the periodic boundary conditions, which for example are generally introduced when triangular lattices are used in simulations.
△ Less
Submitted 9 November, 1998;
originally announced November 1998.
-
Universality of the excess number of clusters and the crossing probability function in three-dimensional percolation
Authors:
Christian D. Lorenz,
Robert M. Ziff
Abstract:
Extensive Monte-Carlo simulations were performed to evaluate the excess number of clusters and the crossing probability function for three-dimensional percolation on the simple cubic (s.c.), face-centered cubic (f.c.c.), and body-centered cubic (b.c.c.) lattices. Systems L x L x L' with L' >> L were studied for both bond (s.c., f.c.c., b.c.c.) and site (f.c.c.) percolation. The excess number of…
▽ More
Extensive Monte-Carlo simulations were performed to evaluate the excess number of clusters and the crossing probability function for three-dimensional percolation on the simple cubic (s.c.), face-centered cubic (f.c.c.), and body-centered cubic (b.c.c.) lattices. Systems L x L x L' with L' >> L were studied for both bond (s.c., f.c.c., b.c.c.) and site (f.c.c.) percolation. The excess number of clusters $\tilde {b}$ per unit length was confirmed to be a universal quantity with a value $\tilde {b} \approx 0.412$. Likewise, the critical crossing probability in the L' direction, with periodic boundary conditions in the L x L plane, was found to follow a universal exponential decay as a function of r = L'/L for large r. Simulations were also carried out to find new precise values of the critical thresholds for site percolation on the f.c.c. and b.c.c. lattices, yielding $p_c(f.c.c.)= 0.199 236 5 \pm 0.000 001 0$, $p_c(b.c.c.)= 0.245 961 5\pm 0.000 001 0$.
△ Less
Submitted 29 June, 1998; v1 submitted 18 June, 1998;
originally announced June 1998.
-
Precise determination of the bond percolation thresholds and finite-size scaling corrections for the s.c., f.c.c., and b.c.c. lattices
Authors:
Christian D. Lorenz,
Robert M. Ziff
Abstract:
Extensive Monte-Carlo simulations were performed to study bond percolation on the simple cubic (s.c.), face-centered cubic (f.c.c.), and body-centered cubic (b.c.c.) lattices, using an epidemic kind of approach. These simulations provide very precise values of the critical thresholds for each of the lattices: pc(s.c.) = 0.248 812 6(5), pc(f.c.c.) = 0.120 163 5(10), and pc(b.c.c.) = 0.180 287 5(1…
▽ More
Extensive Monte-Carlo simulations were performed to study bond percolation on the simple cubic (s.c.), face-centered cubic (f.c.c.), and body-centered cubic (b.c.c.) lattices, using an epidemic kind of approach. These simulations provide very precise values of the critical thresholds for each of the lattices: pc(s.c.) = 0.248 812 6(5), pc(f.c.c.) = 0.120 163 5(10), and pc(b.c.c.) = 0.180 287 5(10). For p close to pc, the results follow the expected finite-size and scaling behavior, with values for the Fisher exponent $tau$ (2.189(2)), the finite-size correction exponent $omega$ (0.64(2)), and the scaling function exponent $sigma$ (0.445(1)) confirmed to be universal.
△ Less
Submitted 3 October, 1997;
originally announced October 1997.