Source author record

Søren Toxvaerd

Søren Toxvaerd appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

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

13 published item(s)

preprint2023arXiv

The formation and stability of homochiral peptides in aqueous prebiological environment in the Earth's crust

The oldest forms of living organisms on Earth are about 3.5 billion years old, and they are found in hydrothermal deposits, and it is often hypothesized that life originated there. But the hydrothermal systems with a fairly strong flow of chemical components are not the optimal place for the prebiological self-assembly of biomolecules and for the emergence of homochirality. This article examines the possibility for that the self-assembly of homochiral molecules took place in an aqueous environment in the Earth's crust. Based on the latest literature regarding the conditions in the lithosphere there are several factors that point to that the crust could be the location for the prebiological self-assembly of biomolecules, and there is nothing against it. The crust and the mantle contain a substantial amount of water, and at the time prior to the emergence of life the crust contained most likely the necessary chemical substances for the synthesis of the biomolecules and an aqueous environments where the homochirality could be established.

preprint2022arXiv

An algorithm for coalescence of classical objects and formation of planetary systems

Isaac Newton formulated the central difference algorithm (Eur. Phys. J. Plus (2020) 135:267) when he derived his second law. The algorithm is under various names ("Verlet, leap-frog,...") the most used algorithm in simulations of complex systems in Physics and Chemistry, and it is also applied in Astrophysics. His discrete dynamics has the same qualities as his exact analytic dynamics for contineus space and time with time reversibility, symplecticity and conservation of momentum, angular momentum and energy. Here the algorithm is extended to include the fusion of objects at collisions. The extended algorithm is used to obtain the self-assembly of celestial objects at the emergence of planetary systems. The emergence of twelve planetary systems is obtained. The systems are stable over very long times, even when two "planets" collide or if a planet is engulfed by its sun.

preprint2022arXiv

Planetary systems with forces other than gravitational forces

A discrete and exact algorithm for obtaining planetary systems is derived in a recent article (Eur. Phys. J. Plus 2022, 137:99). Here the algorithm is used to obtain planetary systems with forces different from the Newtonian inverse square gravitational forces. A Newtonian planetary system exhibits regular elliptical orbits, and here it is demonstrated that a planetary system with pure inverse forces also is stable and with regular orbits, whereas a planetary system with inverse cubic forces is unstable and without regular orbits. The regular orbits in a planetary system with inverse forces deviate, however, from the usual elliptical orbits by having revolving orbits with tendency to orbits with three or eight loops. Newton's Proposition 45 in $\textit{Principia}$ for the Moon's revolving orbits caused by an additional attraction to the gravitational attraction is confirmed, but whereas the additional inverse forces stabilize the planetary system, the additional inverse cubic forces can destabilize the planetary system at a sufficient strength.

preprint2021arXiv

The role of local bond-order at crystallization in a simple supercooled liquid

Large scale Molecular Dynamics simulations of sixty-five systems with N=80000 Lennard-Jones particles at two different supercooled liquid state points reveals, that the supercooled states contain spatially heterogeneous distributed subdomains of particles with significant higher bond-order than the mean bond-order in the supercooled liquid. The onset of crystallization starts in such an area with relative high six-ford bond-order for a supercooled state, but low bond-order for a fcc crystal. The critical nucleus of N \approx 70 particles is surrounded by many hundred of particles with relative high supercooled liquid bond-order and many of these particles are aligned with the crystal ordered particles in the critical nucleus. The crystallizations are very fast and supported by a fast growth of the supercooled areas with relative high liquid bond-order. The crystallization are to fcc crystals, but with significant part of the crystallizations exhibit five fold arrangements of polycrystalline subdomains mainly with fcc crystal order and sign of hcp crystallites.

preprint2020arXiv

Newton's discrete dynamics

In 1687 Isaac Newton published PHILOSOPHIÆ\ NATURALIS PRINCIPIA MATHEMATICA, where the classical analytic dynamics was formulated. But Newton also formulated a discrete dynamics, which is the central difference algorithm, known as the Verlet algorithm. In fact Newton used the central difference to derive his second law. The central difference algorithm is used in computer simulations,where almost all Molecular Dynamics simulations are performed with the Verlet algorithm or other reformulations of the central difference algorithm. Here we show, that the discrete dynamics obtained by Newtons algorithm for Kepler's equation has the same solutions as the analytic dynamics. The discrete positions of a celestial body are located on an ellipse, which is the exact solution for a shadow Hamiltonian nearby the Hamiltonian for the analytic solution.

preprint2016arXiv

Nucleation and droplet growth from supersaturated vapor at temperatures below the triple point temperature

In 1897 Ostwald formulated his step rule for formation of the most stable crystal state for a system with crystal polymorphism. The rule describes the irreversible way a system converts to the crystal with lowest free energy. But in fact the irreversible way a supercooled gas below the triple point temperature $T_{tr.p.}$ crystallizes via a liquid droplet is an example of Ostwald's step rule. The homogeneous nucleation in the supersaturated gas is not to a crystal, but to a liquid-like critical nucleus. We have for the first time performed constant energy (NVE) Molecular Dynamics (MD) of homogeneous nucleation without the use of a thermostat. The simulations of homogeneous nucleation in a Lennard-Jones system from supersaturated vapor at temperatures below $T_{tr.p.}$ reveals that the nucleation to a liquid-like critical nucleus is initiated by a small cold cluster [S. Toxvaerd, J. Chem. Phys. \textbf{143} 154705 (2015)]. The release of latent heat at the subsequent droplet growth increases the temperature in the liquid-like droplet, which for not deep supercooling and/or low supersaturation, can exceed $T_{tr.p.}$. The temperature of the liquid-like droplet increases less for a low supersaturation and remains below $T_{tr.p.}$, but without a crystallization of the droplet for long times. The dissipation of the latent heat into the surrounding gas is affected by a traditional MD thermostat, with the consequence that droplet growth is different for (NVE) MD and constant temperature (NVT) MD.

preprint2015arXiv

Dynamics of homogeneous nucleation

The classical nucleation theory for homogeneous nucleation is formulated as a theory for a density fluctuation in a supersaturated gas at a given temperature. But Molecular Dynamics simulations reveal that it is small cold clusters which initiates the nucleation. The temperature in the nucleating clusters fluctuate, but the mean temperature remains below the temperature in the supersaturated gas until they reach the critical nucleation size. The critical nuclei have, however, a temperature equal to the supersaturated gas. The kinetics of homogeneous nucleation is not only caused by a grow or shrink by accretion or evaporation of monomers only, but by an exponentially declining change in cluster size per time step equal to the cluster distribution in the supersaturated gas.

preprint2012arXiv

Do the repulsive and attractive pair forces play separate roles for the physics of liquids?

According to standard liquid-state theory repulsive and attractive pair forces play distinct roles for the physics of liquids. This paradigm is put into perspective here by demonstrating a continuous series of pair potentials that have virtually the same structure and dynamics, although only some of them have attractive forces of significance. Our findings reflect the fact that the motion of a given particle is determined by the total force on it, whereas the quantity usually discussed in liquid-state theory is the individual pair force.

preprint2011arXiv

NVU dynamics. I. Geodesic motion on the constant-potential-energy hypersurface

An algorithm is derived for computer simulation of geodesics on the constant potential-energy hypersurface of a system of N classical particles. First, a basic time-reversible geodesic algorithm is derived by discretizing the geodesic stationarity condition and implementing the constant potential energy constraint via standard Lagrangian multipliers. The basic NVU algorithm is tested by single-precision computer simulations of the Lennard-Jones liquid. Excellent numerical stability is obtained if the force cutoff is smoothed and the two initial configurations have identical potential energy within machine precision. Nevertheless, just as for NVE algorithms, stabilizers are needed for very long runs in order to compensate for the accumulation of numerical errors that eventually lead to "entropic drift" of the potential energy towards higher values. A modification of the basic NVU algorithm is introduced that ensures potential-energy and step-length conservation; center-of-mass drift is also eliminated. Analytical arguments confirmed by simulations demonstrate that the modified NVU algorithm is absolutely stable. Finally, simulations show that the NVU algorithm and the standard leap-frog NVE algorithm have identical radial distribution functions for the Lennard-Jones liquid.

preprint2011arXiv

NVU dynamics. II. Comparing to four other dynamics

In the companion paper [Ingebrigtsen et al., arXiv:1012.3447] an algorithm was developed for tracing out a geodesic curve on the constant-potential-energy hypersurface. Here simulations of this NVU dynamics are compared to results for four other dynamics, both deterministic and stochastic. First, NVU dynamics is compared to the standard energy-conserving Newtonian NVE dynamics by simulations of the Kob-Andersen binary Lennard-Jones liquid, its WCA version (i.e., with cut-off's at the pair potential minima), and the Gaussian Lennard-Jones liquid. We find identical results for all quantities probed: radial distribution functions, incoherent intermediate scattering functions, and mean-square displacement as function of time. Arguments are then presented for the equivalence of NVU and NVE dynamics in the thermodynamic limit; in particular to leading order in 1/N these two dynamics give identical time-autocorrelation functions. In the final section NVU dynamics is compared to Monte Carlo dynamics, to a diffusive dynamics of small-step random walks on the constant-potential-energy hypersurface, and to Nose-Hoover NVT dynamics. If time is scaled for the two stochastic dynamics to make their single-particle diffusion constants identical to those of NVE dynamics, the simulations show that all five dynamics are equivalent at low temperatures except at short times.

preprint2011arXiv

Role of the first coordination shell in determining the equilibrium structure and dynamics of simple liquids

The traditional view that the physical properties of a simple liquid are determined primarily by its repulsive forces was recently challenged by Berthier and Tarjus, who showed that in some cases ignoring the attractions leads to large errors in the dynamics [L. Berthier and G. Tarjus, Phys. Rev. Lett. 103, 170601 (2009); J. Chem. Phys. 134, 214503 (2011)]. We present simulations of the standard Lennard-Jones liquid at several condensed-fluid state points, including a fairly low density state and a very high density state, as well as simulations of the Kob-Andersen binary Lennard-Jones mixture at several temperatures. By varying the range of the forces, results for the thermodynamics, dynamics, and structure show that the determining factor for getting the correct statics and dynamics is not whether or not the attractive forces {\it per se} are included in the simulations. What matters is whether or not interactions are included from all particles within the first coordination shell (FCS) - the attractive forces can thus be ignored, but only at extremely high densities. The recognition of the importance of a local shell in condensed fluids goes back to van der Waals; our results confirm this idea and thereby the basic picture of the old hole- and cell theories for simple condensed fluids.

preprint2011arXiv

Shifted forces in molecular dynamics

Simulations involving the Lennard-Jones potential usually employ a cut-off at $r=2.5σ$. This paper investigates the possibility of reducing the cut-off. Two different cut-off implementations are compared, the standard shifted potential cut-off and the less commonly used shifted forces cut-off. The first has correct forces below the cut-off, whereas the shifted forces cut-off modifies Newton's equations at all distances. The latter is nevertheless superior; we find that for most purposes realistic simulations may be obtained using a shifted-forces cut-off at $r=1.5σ$, even though the pair force is here 30 times larger than at $r=2.5σ$.

preprint2008arXiv

Crystallization of the Kob-Andersen binary Lennard-Jones liquid

At density 1.2 the Kob-Andersen binary Lennard-Jones liquid partly crystallizes in the temperature interval [0.39, 0.44] after typically 10-100 microseconds (Argon units). The crystallization is initiated by a phase separation where the large (A) particles cluster in a volume void of B particles. We investigate a modification of the Kob-Andersen system where the attraction between particles of the same type is removed, thus disfavoring phase separation. We have not been able to crystallize this new system.