BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Niels Grønbech-Jensen

Niels Grønbech-Jensen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.softcond-mat.mtrl-sciphysics.comp-phBiological Physicscond-mat.supr-connlin.CDquant-ph

Catalog footprint

What is connected

27works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

27 published item(s)

preprint2020arXiv

Defining velocities for accurate kinetic statistics in the GJF thermostat

We expand on two previous developments in the modeling of discrete-time Langevin systems. One is the well-documented Grønbech-Jensen Farago (GJF) thermostat, which has been demonstrated to give robust and accurate configurational sampling of the phase space. Another is the recent discovery that also kinetics can be accurately sampled for the GJF method. Through a complete investigation of all possible finite difference approximations to the velocity, we arrive at two main conclusions:~1) It is not possible to define a so-called on-site velocity such that kinetic temperature will be correct and independent of the time step, and~2) there exists a set of infinitely many possibilities for defining a two-point (leap-frog) velocity that measures kinetic energy correctly for linear systems in addition to the correct configurational statistics obtained from the GJF algorithm. We give explicit expressions for the possible definitions, and we incorporate these into convenient and practical algorithmic forms of the normal Verlet-type algorithms along with a set of suggested criteria for selecting a useful definition of velocity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Modeling Escape from a One-Dimensional Potential Well at Zero or Very Low Temperatures

The process of activation out a one-dimensional potential is investigated systematically in zero and nonzero temperature conditions. The features of the potential are traced through statistical escape out of its wells whose depths are tuned in time by a forcing term. The process is carried out on the damped pendulum system imposing specific initial conditions on the potential variable. While for relatively high values of the dissipation the statistical properties follow a behavior that can be derived from the standard Kramers model, decreasing the dissipation we observe responses/deviations which have regular dependencies on initial conditions, temperature, and loss parameter itself. It is shown that failures of the thermal activation model are originated at low temperatures, and very low dissipation, by the initial conditions and intrinsic, namely T=0, characteristic oscillations of the potential-generated dynamical equation.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Effects of Intrinsic Dynamical Ghost Modes in Discrete-Time Langevin Simulations

Using the recently published GJF-2GJ Langevin thermostat, which can produce time-step-independent statistical measures even for large time steps, we analyze and discuss the causes for abrupt deviations in statistical data as the time step is increased for some simulations of nonlinear oscillators. Exemplified by the pendulum, we identify a couple of discrete-time dynamical modes in the purely damped pendulum equation as the cause of the observed discrepancies in statistics. The existence, stability and kinetics of the modes are consistent with the acquired velocity distribution functions from Langevin simulations, and we conclude that the simulation deviations from physical expectations are not due to normal, systematic algorithmic time-step errors, but instead due to the inherent properties of discrete time in nonlinear dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A survey of classical and quantum interpretations of experiments on Josephson junctions at very low temperatures

For decades following its introduction in 1968, the resistively and capacitively shunted junction (RCSJ) model, sometimes referred to as the Stewart-McCumber model, was successfully applied to study the dynamics of Josephson junctions embedded in a variety of superconducting circuits. In 1980 a theoretical conjecture by A.J. Leggett suggested a possible new and quite different behavior for Josephson junctions at very low temperatures. A number of experiments seemed to confirm this prediction and soon it was taken as a given that junctions at tens of millikelvins should be regarded as macroscopic quantum entities. As such, they would possess discrete levels in their effective potential wells, and would escape from those wells (with the appearance of a finite junction voltage) via a macroscopic quantum tunneling process. A zeal to pursue this new physics led to a virtual abandonment of the RCSJ model in this low temperature regime. In this paper we consider a selection of essentially prototypical experiments that were carried out with the intention of confirming aspects of anticipated macroscopic quantum behavior in Josephson junctions. We address two questions: (1) How successful is the non-quantum theory (RCSJ model) in replicating those experiments? (2) How strong is the evidence that data from these same experiments does indeed reflect macroscopic quantum behavior?

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Isothermal Langevin dynamics in systems with power-law spatially-dependent friction

We study the dynamics of Brownian particles in a heterogeneous one-dimensional medium with a spatially-dependent diffusion coefficient of the form $D(x)\sim |x|^c$, at constant temperature. The particle's probability distribution function (PDF) is calculated both analytically, by solving Fick's diffusion equation, and from numerical simulations of the underdamped Langevin equation. At large times, the PDFs calculated by both approaches yield identical results, corresponding to subdiffusion for $c<0$, and superdiffusion for $0<c<1$. For $c>1$, the diffusion equation predicts that the particles accelerate. Here, we show that this phenomenon, previously considered in several works as an illustration for the possible dramatic effects of spatially-dependent thermal noise, is unphysical. We argue that in an isothermal medium, the motion cannot exceed the ballistic limit ($\left\langle x^2\right\rangle \sim t^2$). The ballistic limit is reached when the friction coefficient drops sufficiently fast at large distances from the origin, and is correctly captured by Langevin's equation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the connection between dissipative particle dynamics and the Itô-Stratonovich dilemma

Dissipative Particle Dynamics (DPD) is a popular simulation model for investigating hydrodynamic behavior of systems with non-negligible equilibrium thermal fluctuations. DPD employs soft core repulsive interactions between the system particles, thus allowing them to overlap. This supposedly permits relatively large integration time steps, which is an important feature for simulations on large temporal scales. In practice, however, an increase in the integration time step leads to increasingly larger systematic errors in the sampling statistics. Here, we demonstrate that the prime origin of these systematic errors is the multiplicative nature of the thermal noise term in Langevin's equation, i.e., the fact that it depends on the instantaneous coordinates of the particles. This lead to an ambiguity in the interpretation of the stochastic differential Langevin equation, known as the Itô-Stratonovich dilemma. Based on insights from previous studies of the dilemma, we propose a novel algorithm for DPD simulations exhibiting almost an order of magnitude improvement in accuracy, and nearly twice the efficiency of commonly used DPD Langevin thermostats.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The G-JF Thermostat for Accurate Configurational Sampling in Soft-Matter Simulations

We implement the statistically sound G-JF thermostat for Langevin Dynamics simulations into the ESPREesSo molecular package for large-scale simulations of soft matter systems. The implemented integration method is tested against the integrator currently used by the molecular package in simulations of a fluid bilayer membrane. While the latter exhibits deviations in the sampling statistics that increase with the integration time step dt, the former reproduces near-correct configurational statistics for all dt within the stability range of the simulations. We conclude that, with very modest revisions to existing codes, one can significantly improve the performance of statistical sampling using Langevin thermostats.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Thermodynamics of Endosomal Escape and DNA Release from Lipoplexes

Complexes of cationic and neutral lipids and DNA (lipoplexes) are emerging as promising vectors for gene therapy applications. Their appeal stems from their non pathogenic nature and the fact that they self-assemble under conditions of thermal equilibrium. Lipoplex adhesion to the cell plasma membrane initiates a three-stage process termed transfection, consisting of (i) endocytosis, (ii) lipoplex breakdown, and (iii) DNA release followed by gene expression. As successful transfection requires lipoplex degradation, it tends to be hindered by the lipoplex thermodynamic stability; nevertheless, it is known that the transfection process may proceed spontaneously. Here, we use a simple model to study the thermodynamic driving forces governing transfection. We demonstrate that after endocytosis [stage (i)], the lipoplex becomes inherently unstable. This instability, which is triggered by interactions between the cationic lipids of the lipoplex and the anionic lipids of the enveloping plasma membrane, is entropically controlled involving both remixing of the lipids and counterions release. Our detailed calculation shows that the free energy gain during stage (ii) is approximately linear in $Φ_+$, the mole fraction of cationic lipids in the lipoplex. This free energy gain, $ΔF$, reduces the barrier for fusion between the enveloping and the lipoplex bilayers, which produces a hole allowing for DNA release [stage (iii)]. The linear relationship between $ΔF$ and the fraction of cationic lipids explains the experimentally observed exponential increase of transfection efficiency with $Φ_+$ in lamellar lipoplexes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Constant pressure and temperature discrete-time Langevin molecular dynamics

We present a new and improved method for simultaneous control of temperature and pressure in molecular dynamics simulations with periodic boundary conditions. The thermostat-barostat equations are build on our previously developed stochastic thermostat, which has been shown to provide correct statistical configurational sampling for any time step that yields stable trajectories. Here, we extend the method and develop a set of discrete-time equations of motion for both particle dynamics and system volume in order to seek pressure control that is insensitive to the choice of the numerical time step. The resulting method is simple, practical, and efficient. The method is demonstrated through direct numerical simulations of two characteristic model systems - a one dimensional particle chain for which exact statistical results can be obtained and used as benchmarks, and a three dimensional system of Lennard-Jones interacting particles simulated in both solid and liquid phases. The results, which are compared against the method of Kolb & Dunweg, show that the new method behaves according to the objective, namely that acquired statistical averages and fluctuations of configurational measures are accurate and robust against the chosen time step applied to the simulation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Effects of applied strain on radiation damage generation in body-centered cubic iron

Radiation damage in body-centered cubic (BCC) Fe has been extensively studied by computer simulations to quantify effects of temperature, impinging particle energy, and the presence of extrinsic particles. However, limited investigation has been conducted into the effects of mechanical stresses and strain. In a reactor environment, structural materials are often mechanically strained, and an expanded understanding of how this strain affects the generation of defects may be important for predicting microstructural evolution and damage accumulation under such conditions. In this study, we have performed molecular dynamics simulations in which various types of homogeneous strains are applied to BCC Fe and the effect on defect generation is examined. It is found that volume-conserving shear strains yield no statistically significant variations in the stable number of defects created via cascades in BCC Fe. However, strains that result in volume changes are found to produce significant effects on defect generation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Elasticity and mechanical instability of charged lipid bilayers in ionic solutions

We use coarse-grained Monte Carlo simulations to study the elastic properties of charged membranes in solutions of monovalent and pentavalent counterions. The simulation results of the two cases reveal trends opposite to each other. The bending rigidity and projected area increase with the membrane charge density for monovalent counterions, while they decrease for the pentavalent ions. These observations can be related to the counterion screening of the lipid charges. While the monovalent counterions only weakly screen the Coulomb interactions, which implies a repulsive Coulomb system, the multivalent counterions condense on the membrane and, through spatial charge correlations, make the effective interactions due the charged lipids attractive. The differences in the elastic properties of the charged membranes in monovalent and multivalent counterion solutions are reflected in the mechanisms leading to their mechanical instability at high charge densities. In the former case, the membranes develop pores to relieve the electrostatic tensile stresses, while in the latter case, the membrane exhibit large wavelength bending instability.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fluctuation-dissipation relation for systems with spatially varying friction

When a particle diffuses in a medium with spatially dependent friction coefficient $α(r)$ at constant temperature $T$, it drifts toward the low friction end of the system even in the absence of any real physical force $f$. This phenomenon, which has been previously studied in the context of non-inertial Brownian dynamics, is termed "spurious drift", although the drift is real and stems from an inertial effect taking place at the short temporal scales. Here, we study the diffusion of particles in inhomogeneous media within the framework of the inertial Langevin equation. We demonstrate that the quantity which characterizes the dynamics with non-uniform $α(r)$ is not the displacement of the particle $Δr=r-r^0$ (where $r^0$ is the initial position), but rather $ΔA(r)=A(r)-A(r^0)$, where $A(r)$ is the primitive function of $α(r)$. We derive expressions relating the mean and variance of $ΔA$ to $f$, $T$, and the duration of the dynamics $Δt$. For a constant friction coefficient $α(r)=α$, these expressions reduce to the well known forms of the force-drift and fluctuation-dissipation relations. We introduce a very accurate method for Langevin dynamics simulations in systems with spatially varying $α(r)$, and use the method to validate the newly derived expressions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Switching Current Distributions in Josephson Junctions at Very Low Temperatures

Swept bias experiments carried out on Josephson junctions yield the distributions of the probabilities of early switching from the zero voltage state. Kramers' theory of thermally activated escape from a one-dimensional potential is well known to fall short of explaining such experiments when the junctions are at millikelvin temperatures. We propose a simple revision of the theory which is shown to yield extremely good agreement with experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A simple and effective Verlet-type algorithm for simulating Langevin dynamics

We present a revision to the well known Stormer-Verlet algorithm for simulating second order differential equations. The revision addresses the inclusion of linear friction with associated stochastic noise, and we analytically demonstrate that the new algorithm correctly reproduces diffusive behavior of a particle in a flat potential. For a harmonic oscillator, our algorithm provides the exact Boltzmann distribution for any value of damping, frequency, and time step for both underdamped and over damped behavior within the usual the stability limit of the Verlet algorithm. Given the structure and simplicity of the method we conclude this approach can trivially be adapted for contemporary applications, including molecular dynamics with extensions such as molecular constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Application of the G-JF Discrete-Time Thermostat for Fast and Accurate Molecular Simulations

A new Langevin-Verlet thermostat that preserves the fluctuation-dissipation relationship for discrete time steps, is applied to molecular modeling and tested against several popular suites (AMBER, GROMACS, LAMMPS) using a small molecule as an example that can be easily simulated by all three packages. Contrary to existing methods, the new thermostat exhibits no detectable changes in the sampling statistics as the time step is varied in the entire numerical stability range. The simple form of the method, which we express in the three common forms (Velocity-Explicit, Stormer-Verlet, and Leap-Frog), allows for easy implementation within existing molecular simulation packages to achieve faster and more accurate results with no cost in either computing time or programming complexity.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Langevin dynamics in inhomogeneous media: Re-examining the Itô-Stratonovich Dilemma

The diffusive dynamics of a particle in a medium with space-dependent friction coefficient is studied within the framework of the inertial Langevin equation. In this description, the ambiguous interpretation of the stochastic integral, known as the Itô-Stratonovich dilemma, is avoided since all interpretations converge to the same solution in the limit of small time steps. We use a newly developed method for Langevin simulations to measure the probability distribution of a particle diffusing in a flat potential. Our results reveal that both the Itô and Stratonovich interpretations converge very slowly to the uniform equilibrium distribution for vanishing time step sizes. Three other conventions exhibit significantly improved accuracy: (i) the "isothermal" (Hänggi) convention, (ii) the Stratonovich convention corrected by a drift term, and (iii) a new convention employing two different effective friction coefficients representing two different averages of the friction function during the time step. We argue that the most physically accurate dynamical description is provided by the third convention, in which the particle experiences a drift originating from the dissipation instead of the fluctuation term. This feature is directly related to the fact that the drift is a result of an inertial effect that cannot be well understood in the Brownian, overdamped limit of the Langevin equation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A classical statistical model for distributions of escape events in swept-bias Josephson junctions

We have developed a model for experiments in which the bias current applied to a Josephson junction is slowly increased from zero until the junction switches from its superconducting zero-voltage state, and the bias value at which this occurs is recorded. Repetition of such measurements yields experimentally determined probability distributions for the bias current at the moment of escape. Our model provides an explanation for available data on the temperature dependence of these escape peaks. When applied microwaves are included we observe an additional peak in the escape distributions and demonstrate that this peak matches experimental observations. The results suggest that experimentally observed switching distributions, with and without applied microwaves, can be understood within classical mechanics and may not exhibit phenomena that demand an exclusively quantum mechanical interpretation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Breakdown of the classical description of a local system

We provide a straightforward demonstration of a fundamental difference between classical and quantum mechanics for a single local system; namely the absence of a joint probability distribution of the position $x$ and momentum $p$. Elaborating on a recently reported criterion by Bednorz and Belzig [Phys. Rev. A {\bf 83}, 52113] we derive a simple criterion that must be fulfilled for any joint probability distribution in classical physics. We demonstrate the violation of this criterion using homodyne measurement of a single photon state, thus proving a straightforward signature of the breakdown of a classical description of the underlying state. Most importantly, the criterion used does not rely on quantum mechanics and can thus be used to demonstrate non-classicality of systems not immediately apparent to exhibit quantum behavior. The criterion is directly applicable any system described by the continuous canonical variables x and p, such as a mechanical or an electrical oscillator and a collective spin of a large ensemble.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Molecular simulation analysis of structural variations in lipoplexes

We use a coarse-grained molecular model to study the self-assembly process of complexes of cationic and neutral lipids with DNA molecules ("lipoplexes") - a promising nonviral carrier of DNA for gene therapy. We identify the resulting structures through direct visualization of the molecular arrangements and through calculations of the corresponding scattering plots. The latter approach provides a means for comparison with published data from X-ray scattering experiments. Consistent with experimental results, we find that upon increasing the stiffness of the lipid material, the system tends to form lamellar structures. Two characteristic distances can be extracted from the scattering plots of lamellar complexes - the lamellar (interlayer) spacing and the DNA-spacing within each layer. We find a remarkable agreement between the computed values of these two quantities and the experimental data [J. O. Rädler, I. Koltover, T. Salditt and C. R. Safinya, Science Vol. 275, 810 (1997)] over the entire range of mole fractions of charged lipids (CLs) studied experimentally. A visual inspection of the simulated systems reveals that, for very high fractions of CLs, disordered structures consisting of DNA molecules bound to small membrane fragments are spontaneously formed. The diffraction plots of these non-lamellar disordered complexes appear very similar to that of the lamellar structure, which makes the interpretation of the X-ray data ambiguous. The loss of lamellar order may be the origin of the observed increase in the efficiency of lipoplexes as gene delivery vectors at high charge densities.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Critical velocities for deflagration and detonation triggered by voids in a REBO high explosive

The effects of circular voids on the shock sensitivity of a two-dimensional model high explosive crystal are considered. We simulate a piston impact using molecular dynamics simulations with a Reactive Empirical Bond Order (REBO) model potential for a sub-micron, sub-ns exothermic reaction in a diatomic molecular solid. The probability of initiating chemical reactions is found to rise more suddenly with increasing piston velocity for larger voids that collapse more deterministically. A void with radius as small as 10 nm reduces the minimum initiating velocity by a factor of 4.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Interatomic potentials for mixed oxide (MOX) nuclear fuels

We extend our recently developed interatomic potentials for UO_{2} to the mixed oxide fuel system (U,Pu,Np)O_{2}. We do so by fitting against an extensive database of ab initio results as well as to experimental measurements. The applicability of these interactions to a variety of mixed environments beyond the fitting domain is also assessed. The employed formalism makes these potentials applicable across all interatomic distances without the need for any ambiguous splining to the well-established short-range Ziegler-Biersack-Littmark universal pair potential. We therefore expect these to be reliable potentials for carrying out damage simulations (and Molecular Dynamics simulations in general) in nuclear fuels of varying compositions for all relevant atomic collision energies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Sensitivity effects of void density and arrangement in a REBO high explosive

The shock response of two-dimensional model high explosive crystals with various arrangements of circular voids is explored. We simulate a piston impact using molecular dynamics simulations with a Reactive Empirical Bond Order (REBO) model potential for a sub-micron, sub-ns exothermic reaction in a diatomic molecular solid. In square lattices of voids (of equal size), reducing the size of the voids or increasing the porosity while holding the other parameter fixed causes the hotspots to consume the material more quickly and detonation to occur sooner and at lower piston velocities. The early time behavior is seen to follow a very simple ignition and growth model. The hotspots are seen to collectively develop a broad pressure wave (a sonic, diffuse deflagration front) that, upon merging with the lead shock, transforms it into a detonation. The reaction yields produced by triangular lattices are not significantly different. With random void arrangements, the mean time to detonation is 15.5\% larger than with the square lattice; the standard deviation of detonation delays is just 5.1%.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Simulation of ion track ranges in uranium oxide

Direct comparisons between statistically sound simulations of ion-tracks and published experimental measurements of range densities of iodine implants in uranium dioxide have been made with implant energies in the range of 100-800 keV. Our simulations are conducted with REED-MD (Rare Event Enhanced Domain-following Molecular Dynamics) in order to account for the materials structure in both single crystalline and polycrystalline experimental samples. We find near-perfect agreement between REED-MD results and experiments for polycrystalline target materials.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Is state tomography an unambiguous test of quantum entanglement?

We provide an alternative interpretation of recently published experimental results that were represented as demonstrating entanglement between two macroscopic quantum Josephson oscillators. We model the experimental system using the well-established classical equivalent circuit of a resistively and capacitively shunted junction. Simulation results are used to generate the corresponding density matrix, which is strikingly similar to the previously published matrix that has been declared to be an unambiguous demonstration of quantum entanglement. Since our data are generated by a classical model, we therefore submit that state tomography cannot be used to determine absolutely whether or not quantum entanglement has taken place. Analytical arguments are given for why the classical analysis provides an adequate explanation of the experimental results.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Transitions between distinct compaction regimes in complexes of multivalent cationic lipids and DNA

We use X-ray scattering and molecular simulations to investigate the structural properties of complexes of multivalent cationic lipids and DNA molecules. At low mole fraction of neutral lipids (NLs), $Φ_{\rm NL}$, the complexes show dramatic DNA compaction down to essentially close packed DNA arrays with a DNA interaxial spacing $d_{\rm DNA}=25Å$. A gradual increase in $Φ_{\rm NL}$ does not lead to a continuous increase in $d_{\rm DNA}$ as observed for DNA complexes of monovalent cationic lipids (CLs). Instead, distinct spacing regimes exist, with sharp transitions between them. Three packing states have been identified: (i) close packed, (ii) condensed, but not close packed, with $d_{\rm DNA}=27-28Å$, and (iii) an expanded state, where $d_{\rm DNA}$ increases gradually with $Φ_{\rm NL}$. Based on our experimental and computational results, we conclude that the DNA condensation is mediated by the multivalent cationic lipids, which assemble between the negatively charged DNA rods. Quite remarkably, the computational results show that the less tightly packed structure in regime (ii) is thermodynamically more stable than the close packed structure in regime (i). Accordingly, the constant DNA spacing observed in regime (ii) is attributed to lateral phase coexistence between this stable CL-DNA complex and neutral membranes. This finding may explain the reduced transfection efficiency measured for such complexes: Transfection involves endosomal escape and disassembly of the complex, and these processes are inhibited by the high thermodynamic stability. Our results, which demonstrate the existence of an inverse correlation between the stability and transfection activity of lamellar CL-DNA complexes are, therefore, consistent with a recently proposed model of cellular entry.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Non-iterative and exact method for constraining particles in a linear geometry

We present a practical numerical method for evaluating the Lagrange multipliers necessary for maintaining a constrained linear geometry of particles in dynamical simulations. The method involves no iterations, and is limited in accuracy only by the numerical methods for solving small systems of linear equations. As a result of the non-iterative and exact (within numerical accuracy) nature of the procedure there is no drift in the constrained geometry, and the method is therefore readily applied to molecular dynamics simulations of, e.g., rigid linear molecules or materials of non-spherical grains. We illustrate the approach through implementation in the commonly used second-order velocity explicit Verlet method.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Vacancy assisted arsenic diffusion and time dependent clustering effects in silicon

We present results of kinetic lattice Monte Carlo (KLMC) simulations of substitutional arsenic diffusion in silicon mediated by lattice vacancies. Large systems are considered, with 1000 dopant atoms and long range \textit{ab initio} interactions, to the 18th nearest lattice neighbor, and the diffusivity of each defect species over time is calculated. The concentration of vacancies is greater than equilibrium concentrations in order to simulate conditions shortly after ion implantation. A previously unreported time dependence in the applicability of the pair diffusion model, even at low temperatures, is demonstrated. Additionally, long range interactions are shown to be of critical importance in KLMC simulations; when shorter interaction ranges are considered only clusters composed entirely of vacancies form. An increase in arsenic diffusivity for arsenic concentrations up to $10^{19} \text{cm}^{-3}$ is observed, along with a decrease in arsenic diffusivity for higher arsenic concentrations, due to the formation of arsenic dominated clusters. Finally, the effect of vacancy concentration on diffusivity and clustering is studied, and increasing vacancy concentration is found to lead to a greater number of clusters, more defects per cluster, and a greater vacancy fraction within the clusters.

0 citations0 reviewsNo code linkedOpen access