Molecular simulation study of the heat capacity of metastable water between 100K and 300K

Molecular simulation study of the heat capacity of metastable water between 100K and 300K Molecular simulations have been used to study the heat capacity of metastable liquid water at low temperature adsorbed on a smooth surface. These calculations aim at modelling water properties measured by experiments performed on water films adsorbed on Vycor nanoporous silica at low temperature. In particular, the study focuses on the non-monotonous variation of the heat capacity around between 100 and 300 K

ECHO: an Eulerian Conservative High Order scheme for general relativistic magnetohydrodynamics and magnetodynamics

We present a new numerical code, ECHO, based on an Eulerian Conservative High Order scheme for time dependent three-dimensional general relativistic magnetohydrodynamics (GRMHD) and magnetodynamics (GRMD). ECHO is aimed at providing a shock-capturing conservative method able to work at an arbitrary level of formal accuracy (for smooth flows), where the other existing GRMHD and GRMD schemes yield an overall second order at most. Moreover, our goal is to present a general framework, based on the 3+1 Eulerian formalism, allowing for different sets of equations, different algorithms, and working in a generic space-time metric, so that ECHO may be easily coupled to any solver for Einstein's equations. Various high order reconstruction methods are implemented and a two-wave approximate Riemann solver is used. The induction equation is treated by adopting the Upwind Constrained Transport (UCT) procedures, appropriate to preserve the divergence-free condition of the magnetic field in shock-capturing methods. The limiting case of magnetodynamics (also known as force-free degenerate electrodynamics) is implemented by simply replacing the fluid velocity with the electromagnetic drift velocity and by neglecting the matter contribution to the stress tensor. ECHO is particularly accurate, efficient, versatile, and robust. It has been tested against several astrophysical applications, including a novel test on the propagation of large amplitude circularly polarized Alfven waves. In particular, we show that reconstruction based on a Monotonicity Preserving filter applied to a fixed 5-point stencil gives highly accurate results for smooth solutions, both in flat and curved metric (up to the nominal fifth order), while at the same time providing sharp profiles in tests involving discontinuities.Comment: 20 pages, revised version submitted to A&

The electric dipole moment of the nucleon from simulations at imaginary vacuum angle theta

We compute the electric dipole moment of proton and neutron from lattice QCD simulations with N_f=2 flavors of dynamical quarks at imaginary vacuum angle theta. The calculation proceeds via the CP odd form factor F_3. A novel feature of our calculation is that we use partially twisted boundary conditions to extract F_3 at zero momentum transfer. As a byproduct, we test the QCD vacuum at nonvanishing theta.Comment: 22 pages, 10 figure

Improving the lattice axial vector current

For Wilson and clover fermions traditional formulations of the axial vector current do not respect the continuum Ward identity which relates the divergence of that current to the pseudoscalar density. Here we propose to use a point-split or one-link axial vector current whose divergence exactly satisfies a lattice Ward identity, involving the pseudoscalar density and a number of irrelevant operators. We check in one-loop lattice perturbation theory with SLiNC fermion and gauge plaquette action that this is indeed the case including order $O(a)$ effects. Including these operators the axial Ward identity remains renormalisation invariant. First preliminary results of a nonperturbative check of the Ward identity are also presented.Comment: 7 pages, 3 figures, Proceedings of the 33rd International Symposium on Lattice Field Theory, 14-18 July 2015, Kobe, Japa

Precision electromagnetic structure of decuplet baryons in the chiral regime

The electromagnetic properties of the baryon decuplet are calculated in quenched QCD on a 20^3 x 40 lattice with a lattice spacing of 0.128 fm using the fat-link irrelevant clover (FLIC) fermion action with quark masses providing a pion mass as low as 300 MeV. Magnetic moments and charge radii are extracted from the electric and magnetic form factors for each individual quark sector. From these, the corresponding baryon properties are constructed. We present results for the higher order moments of the spin-3/2 baryons, including the electric quadrupole moment E2 and the magnetic octupole moment M3. The world's first determination of a non-zero M3 form factor for the Delta baryon is presented. With these results we provide a conclusive analysis which shows that decuplet baryons are deformed. We compare the decuplet baryon results from a similar lattice calculation of the octet baryons. We establish that the environment sensitivity is far less pronounced in the case of the decuplet baryons compared to that in the octet baryons. A surprising result is that the charge radii of the decuplet baryons are generally smaller than that of the octet baryons. The magnetic moment of the Delta^+ reveals a turn over in the low quark mass region, making it smaller than the proton magnetic moment. These results are consistent with the expectations of quenched chiral perturbation theory. A similar turn over is also noticed in the magnetic moment of the Sigma^*0, but not for Xi^* where only kaon loops can appear in quenched QCD. The electric quadrupole moment of the Omega^- baryon is positive when the negative charge factor is included, and is equal to 0.86 +- 0.12 x 10^-2 fm^2, indicating an oblate shape.Comment: 30 pages, 32 figure

Some continuum physics results from the lattice V-A correlator

We present preliminary results on extractions of the chiral LECs L_10 and C_87 and constraints on the excited pseudoscalar state pi(1300) and pi(1800) decay constants obtained from an analysis of lattice data for the flavor ud light quark V-A correlator. A comparison of the results for the correlator to the corresponding mildly-model-dependent continuum results (based primarily on experimental hadronic tau decay data) is also givenComment: 7 pages, 3 figures. Prepared for the Proceedings of the 30th International Symposium on Lattice Field Theory, Cairns, Australia, June 24-29, 2012; expanded version of Reference 1

New results from the lattice on the theoretical inputs to the hadronic tau determination of V_us

Recent sum rule determinations of |V_us|, employing flavor-breaking combinations of hadronic tau decay data, are significantly lower than either expectations based on 3-family unitarity or determinations from K_ell3 and Gamma[K_mu2]/Gamma[pi_mu2]. We use lattice data to investigate the accuracy/reliability of the OPE representation of the flavor-breaking correlator combination entering the tau decay analyses. The behavior of an alternate correlator combination, constructed to reduce problems associated with the slow convergence of the D = 2 OPE series, and entering an alternate sum rule requiring both electroproduction cross-section and hadronic tau decay data, is also investigated. Preliminary updates of both analyses, with the lessons learned from the lattice data in mind, are also presented.Comment: 8 pages, 5 figures. Prepared for the proceedings of the 12th International Workshop on Tau Lepton Physics, Sep. 17-21, 2012, Nagoya, Japan and the 10th International Conference on Confinement and the Hadron Spectrum, Oct. 6-13, 2012, Garching/Munich, German