Non-equilibrium statistical mechanics of classical nuclei interacting with the quantum electron gas

Kinetic equations governing time evolution of positions and momenta of atoms in extended systems are derived using quantum-classical ensembles within the Non-Equilibrium Statistical Operator Method (NESOM). Ions are treated classically, while their electrons quantum mechanically; however, the statistical operator is not factorised in any way and no simplifying assumptions are made concerning the electronic subsystem. Using this method, we derive kinetic equations of motion for the classical degrees of freedom (atoms) which account fully for the interaction and energy exchange with the quantum variables (electrons). Our equations, alongside the usual Newtonian-like terms normally associated with the Ehrenfest dynamics, contain additional terms, proportional to the atoms velocities, which can be associated with the electronic friction. Possible ways of calculating the friction forces which are shown to be given via complicated non-equilibrium correlation functions, are discussed. In particular, we demonstrate that the correlation functions are directly related to the thermodynamic Matsubara Green's functions, and this relationship allows for the diagrammatic methods to be used in treating electron-electron interaction perturbatively when calculating the correlation functions. This work also generalises previous attempts, mostly based on model systems, of introducing the electronic friction into Molecular Dynamics equations of atoms.Comment: 18 page

Calculation of electron density of periodic systems using non-orthogonal localised orbitals

Methods for calculating an electron density of a periodic crystal constructed using non-orthogonal localised orbitals are discussed. We demonstrate that an existing method based on the matrix expansion of the inverse of the overlap matrix into a power series can only be used when the orbitals are highly localised (e.g. ionic systems). In other cases including covalent crystals or those with an intermediate type of chemical bonding this method may be either numerically inefficient or fail altogether. Instead, we suggest an exact and numerically efficient method which can be used for orbitals of practically arbitrary localisation. Theory is illustrated by numerical calculations on a model system.Comment: 12 pages, 4 figure

Mobility coefficients in the systems of magnetic dipolar particles

In this paper, we present our first results on the mobility coefficients in the systems of magnetic dipolar particles. In our study, we investigate the influence of chain formation and polydispersity of particles on self-diffusion. The work is purely theoretical and combines direct calculations with the density functional approach to calculate equilibrium densities of chains. We study mainly bulk systems. It is shown that the formation of chains leads to the average decrease of mobility in monodisperse systems, but in the case of bidisperse particle size distribution, the particle mobility becomes a function of the fractional composition. The mobility coefficients obtained here are important for calculating the diffusion coefficients in case of gradient-induced diffusion (be that of the field or density gradient) in magnetic fluids with chain aggregates

Nonequilibrium processes from Generalised Langevin Equations: realistic nanoscale systems connected to two thermal baths

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.Comment: present version accepted for publication in Phys. Rev. B (Apr 2016

Studying synthesis confinement effects on the internal structure of nanogels in computer simulations

We study the effects of droplet finite size on the structure of nanogel particles synthesized by random crosslinking of molecular polymers diluted in nanoemulsions. For this, we use a bead-spring computer model of polymer-like structures that mimics the confined random crosslinking process corresponding to irradiation- or electrochemically-induced crosslinking methods. Our results indicate that random crosslinking under strong confinement can lead to unusual nanogel internal structures, with a central region less dense than the external one, whereas under moderate confinement the resulting structure has a denser central region. We analyze the topology of the polymer networks forming nanogel particles with both types of architectures, their overall structural parameters, their response to the quality of the solvent and compare the cases of non-ionic and ionic systems

Nonequilibrium Generalised Langevin Equation for the calculation of heat transport properties in model 1D atomic chains coupled to two 3D thermal baths

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 $\Delta 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 ($\Delta 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 ($\Delta 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.Comment: Accepted for publication in J. Chem. Phy

Curved Noncommutative Tori as Leibniz Quantum Compact Metric Spaces

We prove that curved noncommutative tori, introduced by Dabrowski and Sitarz, are Leibniz quantum compact metric spaces and that they form a continuous family over the group of invertible matrices with entries in the commutant of the quantum tori in the regular representation, when this group is endowed with a natural length function.Comment: 16 Pages, v3: accepted in Journal of Math. Physic

Kinetic and ion pairing contributions in the dielectric spectra of electrolyte aqueous solutions

Understanding dielectric spectra can reveal important information about the dynamics of solvents and solutes from the dipolar relaxation times down to electronic ones. In the late 1970s, Hubbard and Onsager predicted that adding salt ions to a polar solution would result in a reduced dielectric permittivity that arises from the unexpected tendency of solvent dipoles to align opposite to the applied field. So far, this effect has escaped an experimental verification, mainly because of the concomitant appearance of dielectric saturation from which the Hubbard-Onsager decrement cannot be easily separated. Here we develop a novel non-equilibrium molecular dynamics simulation approach to determine this decrement accurately for the first time. Using a thermodynamic consistent all-atom force field we show that for an aqueous solution containing sodium chloride around 4.8 Mol/l, this effect accounts for 12\% of the total dielectric permittivity. The dielectric decrement can be strikingly different if a less accurate force field for the ions is used. Using the widespread GROMOS parameters, we observe in fact an {\it increment} of the dielectric permittivity rather than a decrement. We can show that this increment is caused by ion pairing, introduced by a too low dispersion force, and clarify the microscopic connection between long-living ion pairs and the appearance of specific features in the dielectric spectrum of the solution

On the Regularity of Optimal Transportation Potentials on Round Spheres

In this paper the regularity of optimal transportation potentials defined on round spheres is investigated. Specifically, this research generalises the calculations done by Loeper, where he showed that the strong (A3) condition of Trudinger and Wang is satisfied on the round sphere, when the cost-function is the geodesic distance squared. In order to generalise Loeper's calculation to a broader class of cost-functions, the (A3) condition is reformulated via a stereographic projection that maps charts of the sphere into Euclidean space. This reformulation subsequently allows one to verify the (A3) condition for any case where the cost-fuction of the associated optimal transportation problem can be expressed as a function of the geodesic distance between points on a round sphere. With this, several examples of such cost-functions are then analysed to see whether or not they satisfy this (A3) condition.Comment: 24 pages, 4 figure