Free energy surface reconstruction from umbrella samples using Gaussian process regression II: Multiple collective variables

We demonstrate how a prior assumption of smoothness can be used to enhance the reconstruction of free energy profiles from multiple umbrella sampling simulations using the Bayesian Gaussian process regression approach. The method we derive allows the concurrent use of histograms and free energy gradients and can easily be extended to include further data. In Part I we review the necessary theory and test the method for one collective variable. We demonstrate improved performance with respect to the weighted histogram analysis method and obtain meaningful error bars without any significant additional computation. In Part II we consider the case of multiple collective variables and compare to a reconstruction using least squares fitting of radial basis functions. We find substantial improvements in the regimes of spatially sparse data or short sampling trajectories. A software implementation is made available on www.libatoms.org.

Reactive Many-Body Expansion for a Protonated Water Cluster.

We generalize the standard many-body expansion technique that is used to approximate the total energy of a molecular system to enable the treatment of chemical reactions by quantum chemical techniques. By considering all possible assignments of atoms to monomer units of the many-body expansion and associating suitable weights with each, we construct a potential energy surface that is a smooth function of the nuclear positions. We derive expressions for this reactive many-body expansion energy and describe an algorithm for its evaluation, which scales polynomially with system size, and therefore will make the method feasible for future condensed phase simulations. We demonstrate the accuracy and smoothness of the resulting potential energy surface on a molecular dynamics trajectory of the protonated water hexamer, using the Hartree-Fock method for the many-body term and Møller-Plesset theory for the low order terms of the many-body expansion.

First principles interatomic potential for tungsten based on Gaussian process regression

An accurate description of atomic interactions, such as that provided by first principles quantum mechanics, is fundamental to realistic prediction of the properties that govern plasticity, fracture or crack propagation in metals. However, the computational complexity associated with modern schemes explicitly based on quantum mechanics limits their applications to systems of a few hundreds of atoms at most. This thesis investigates the application of the Gaussian Approximation Potential (GAP) scheme to atomistic modelling of tungsten - a bcc transition metal which exhibits a brittle-to-ductile transition and whose plasticity behaviour is controlled by the properties of $\frac{1}{2} \langle 111 \rangle$ screw dislocations. We apply Gaussian process regression to interpolate the quantum-mechanical (QM) potential energy surface from a set of points in atomic configuration space. Our training data is based on QM information that is computed directly using density functional theory (DFT). To perform the fitting, we represent atomic environments using a set of rotationally, permutationally and reflection invariant parameters which act as the independent variables in our equations of non-parametric, non-linear regression. We develop a protocol for generating GAP models capable of describing lattice defects in metals by building a series of interatomic potentials for tungsten. We then demonstrate that a GAP potential based on a Smooth Overlap of Atomic Positions (SOAP) covariance function provides a description of the $\frac{1}{2} \langle 111 \rangle$ screw dislocation that is in agreement with the DFT model. We use this potential to simulate the mobility of $\frac{1}{2} \langle 111 \rangle$ screw dislocations by computing the Peierls barrier and model dislocation-vacancy interactions to QM accuracy in a system containing more than 100,000 atoms.

A lattice dynamical treatment for the total potential energy of single-walled carbon nanotubes and its applications: relaxed equilibrium structure, elastic properties, and vibrational modes of ultra-narrow tubes

In this paper, we propose a lattice dynamic treatment for the total potential energy of single-walled carbon nanotubes (SWCNTs) which is, apart from a parameter for the nonlinear effects, extracted from the vibrational energy of the planar graphene sheet. The energetics, elasticity and lattice dynamics are treated in terms of the same set of force constants, independently of the tube structures. Based upon this proposal, we have investigated systematically the relaxed lattice configuration for narrow SWCNTs, the strain energy, the Young's modulus and Poisson ratio, and the lattice vibrational properties with respect to the relaxed equilibrium tubule structure. Our calculated results for various physical quantities are nicely in consistency with existing experimental measurements. In particular, we verified that the relaxation effect makes the bond length longer and the frequencies of various optical vibrational modes softer. Our calculation provides evidence that the Young's modulus of an armchair tube exceeds that of the planar graphene sheet, and that the large diameter limits of the Young's modulus and Poisson ratio are in agreement with the experimental values of graphite; the calculated radial breathing modes for ultra-narrow tubes with diameters ranging between 2 and 5 angstrom coincide with the experimental results and the existing ab initio calculations with satisfaction. For narrow tubes with a diameter of 20 angstrom, the calculated frequencies of optical modes in the tubule's tangential plane, as well as those of radial breathing modes, are also in good agreement with the experimental measurements. In addition, our calculation shows that various physical quantities of relaxed SWCNTs can actually be expanded in terms of the chiral angle defined for the corresponding ideal SWCNTs.

Circumstantial Evidence for a Soft Nuclear Symmetry Energy at Suprasaturation Densities

Within an isospin- and momentum-dependent hadronic transport model, it is shown that the recent FOPI data on the pi(-)/pi(+) ratio in central heavy-ion collisions at SIS/GSI energies [Willy Reisdorf , Nucl. Phys. A 781, 459 (2007)] provide circumstantial evidence suggesting a rather soft nuclear symmetry energy E-sym(rho) at rho >= 2 rho(0) compared to the Akmal-Pandharipande-Ravenhall prediction. Some astrophysical implications and the need for further experimental confirmations are discussed.

Shell-model Hamiltonian from self-consistent mean-field model: N = Z nuclei

We propose a procedure to determine the effective nuclear shell-model Hamiltonian in a truncated space from a self-consistent mean-field model, e.g., the Skyrme model. The parameters of pairing plus quadrupole-quadrupole interaction with monopole force are obtained so that the potential energy surface of the Skyrme Hartree-Fock + BCS calculation is reproduced. We test our method for N = Z nuclei in the fpg- and sd-shell regions. It is shown that the calculated energy spectra with these parameters are in a good agreement with experimental data, in which the importance of the monopole interaction is discussed. This method may represent a practical way of defining the Hamiltonian for general shell-model calculations. (C) 2009 Elsevier B.V. All rights reserved.

Triton-He-3 relative and differential flows as probes of the nuclear symmetry energy at supra-saturation densities

Using a transport model coupled with a phase-space coalescence afterburner, we study the triton-He-3 (t-He-3) ratio with both relative and differential transverse flows in semicentral Sn-132 + Sn-124 reactions at a beam energy of 400 MeV/nucleon. The neutron-proton ratios with relative and differential flows are also discussed as a reference. We find that similar to the neutron-proton pairs, the t-He-3 pairs also carry interesting information regarding the density dependence of the nuclear symmetry energy. Moreover, the nuclear symmetry energy affects more strongly the t-He-3 relative and differential flows than the pi(-)/pi(+) ratio in the same reaction. The t-He-3 relative flow can be used as a particularly powerful probe of the high-density behavior of the nuclear symmetry energy.

DETERMINING THE DENSITY DEPENDENCE OF THE NUCLEAR SYMMETRY ENERGY USING HEAVY-ION REACTIONS

We review recent progress in the determination of the subsaturation density behavior of the nuclear symmetry energy from heavy-ion collisions as well as the theoretical progress in probing the high density behavior of the symmetry energy in heavy-ion reactions induced by high energy radioactive beams. We further discuss the implications of these results for the nuclear effective interactions and the neutron skin thickness of heavy nuclei.

The neutron/proton ratio of squeezed-out nucleons and the high density behavior of the nuclear symmetry energy

Within a transport model it is shown that the neutron/proton ratio of squeezed-out nucleons perpendicular to the reaction plane, especially at high transverse momenta, in heavy-ion reactions induced by high energy neutron-rich nuclei can be a useful tool for studying the high density behavior of the nuclear symmetry energy.

Double neutron-proton differential transverse flow as a probe for the high density behavior of the nuclear symmetry energy

The double neutron-proton differential transverse flow taken from two reaction systems using different isotopes of the same element is studied at incident beam energies of 400 and 800 MeV/nucleon within the framework of an isospin- and momentum-dependent hadronic transport model IBUU04. The double differential flow is found to retain about the same sensitivity to the density dependence of the nuclear symmetry energy as the single differential flow in the more neutron-rich reaction. Because the double differential flow reduces significantly both the systematic errors and the influence of the Coulomb force, it is thus more effective probe for the high-density behavior of the nuclear symmetry energy.

Single and double pi(-)/pi(+) ratios in heavy-ion reactions as probes of the high-density behavior of the nuclear symmetry energy

Based on the isospin- and momentum-dependent hadronic transport model IBUU04, effects of the nuclear symmetry energy on the single and double pi(-)/pi(+) ratios in central reactions of Sn-132+Sn-124 and Sn-112+Sn-112 at a beam energy of 400 MeV/nucleon are studied. It is found that around the Coulomb peak of the single pi(-)/pi(+) ratio the double pi(-)/pi(+) ratio taken from the two isotopic reactions retains about the same sensitivity to the density dependence of nuclear symmetry energy. Because the double pi(-)/pi(+) ratio can significantly reduce the systematic errors, it is thus a more effective probe for the high-density behavior of the nuclear symmetry energy.

Double neutron/proton ratio of nucleon emissions in isotopic reaction systems as a robust probe of nuclear symmetry energy

The double neutron/proton ratio of nucleon emissions taken from two reaction systems using four isotopes of the same element, namely, the neutron/proton ratio in the neutron-rich system over that in the more symmetric system, has the advantage of reducing systematically the influence of the Coulomb force and the normally poor efficiencies of detecting low energy neutrons. The double ratio thus suffers less systematic errors. Within the IBUU04 transport model the double neutron/proton ratio is shown to have about the same sensitivity to the density dependence of nuclear symmetry energy as the single neutron/proton ratio in the neutron-rich system involved. The double neutron/proton ratio is therefore more useful for further constraining the symmetry energy of neutron-rich matter.

System size and beam energy effects on probing the high-density behavior of nuclear symmetry energy with pion ratio

Based on the isospin-and momentum-dependent hadronic transport model IBUU04, we have investigated the pi(-)/pi(+) ratio in the following three reactions: Ca-48+Ca-48, Sn-124 +Sn-124 and Au-197+Au-197 with nearly the same isospin asymmetry but different masses, at the bombarding energies from 0.25 to 0.6 A GeV. It is shown that the sensitivity of probing the E-sym (rho) with pi(-)/pi(+) increases with increasing the system size or decreasing the beam energy, showing a correlation to the degree of isospin fractionation. Therefore, with a given isospin asymmetry, heavier system at energies near the pion threshold is preferential to study the behavior Of nuclear symmetry energy at supra-saturation densities.

TRITON-(HE)-H-3 RELATIVE AND DIFFERENTIAL FLOWS AND THE HIGH DENSITY BEHAVIOR OF NUCLEAR SYMMETRY ENERGY

Using a transport model coupled with a phase-space coalescence after-burner we study the triton-He-3 relative and differential transverse flows in semi-central Sn-132 + Sn-124 reactions at a beam energy of 400 MeV/nucleon. We find that the triton-He-3 pairs carry interesting information about the density dependence of the nuclear symmetry energy. The t-He-3 relative flow can be used as a particularly powerful probe of the high-density behavior of the nuclear symmetry energy.