Algorithmic differentiation and the calculation of forces by quantum Monte Carlo

We describe an efficient algorithm to compute forces in quantum Monte Carlo using adjoint algorithmic differentiation. This allows us to apply the space warp coordinate transformation in differential form, and compute all the 3M force components of a system with M atoms with a computational effort comparable with the one to obtain the total energy. Few examples illustrating the method for an electronic system containing several water molecules are presented. With the present technique, the calculation of finite-temperature thermodynamic properties of materials with quantum Monte Carlo will be feasible in the near future.Comment: 32 pages, 4 figure, to appear in The Journal of Chemical Physic

Ising transition in the two-dimensional quantum J1−J2J_1-J_2 Heisenberg model

We study the thermodynamics of the spin-$S$ two-dimensional quantum Heisenberg antiferromagnet on the square lattice with nearest ($J_1$) and next-nearest ($J_2$) neighbor couplings in its collinear phase ($J_2/J_1>0.5$), using the pure-quantum self-consistent harmonic approximation. Our results show the persistence of a finite-temperature Ising phase transition for every value of the spin, provided that the ratio $J_2/J_1$ is greater than a critical value corresponding to the onset of collinear long-range order at zero temperature. We also calculate the spin- and temperature-dependence of the collinear susceptibility and correlation length, and we discuss our results in light of the experiments on Li$_2$VOSiO$_4$ and related compounds.Comment: 4 page, 4 figure

Spin-lattice coupling in frustrated antiferromagnets

We review the mechanism of spin-lattice coupling in relieving the geometrical frustration of pyrochlore antiferromagnets, in particular spinel oxides. The tetrahedral unit, which is the building block of the pyrochlore lattice, undergoes a spin-driven Jahn-Teller instability when lattice degrees of freedom are coupled to the antiferromagnetism. By restricting our considerations to distortions which preserve the translational symmetries of the lattice, we present a general theory of the collective spin-Jahn-Teller effect in the pyrochlore lattice. One of the predicted lattice distortions breaks the inversion symmetry and gives rise to a chiral pyrochlore lattice, in which frustrated bonds form helices with a definite handedness. The chirality is transferred to the spin system through spin-orbit coupling, resulting in a long-period spiral state, as observed in spinel CdCr2O4. We discuss explicit models of spin-lattice coupling using local phonon modes, and their applications in other frustrated magnets.Comment: 23 pages, 6 figures. Lecture notes for Trieste Summer School, August 2007. To appear as a chapter in "Highly Frustrated Magnetism", Eds. C. Lacroix, P. Mendels, F. Mil

Reentrant behavior of the phase stiffness in Josephson junction arrays

The phase diagram of a 2D Josephson junction array with large substrate resistance, described by a quantum XY model, is studied by means of Fourier path-integral Monte Carlo. A genuine Berezinskii-Kosterlitz-Thouless transition is found up to a threshold value g* of the quantum coupling, beyond which no phase coherence is established. Slightly below g* the phase stiffness shows a reentrant behavior with temperature, in connection with a low-temperature disappearance of the superconducting phase, driven by strong nonlinear quantum fluctuations.Comment: 4 pages, 7 figures, to appear in Phys.Rev.Let

Suppression of Dimer Correlations in the Two-Dimensional J1J_1-J2J_2 Heisenberg Model: an Exact Diagonalization Study

We present an exact diagonalization study of the ground state of the spin-half $J_1{-}J_2$ model. Dimer correlation functions and the susceptibility associated to the breaking of the translational invariance are calculated for the $4\times 4$ and the $6\times 6$ clusters. These results -- especially when compared to the one dimensional case, where the occurrence of a dimerized phase for large enough frustration is well established -- suggest either a homogeneous spin liquid or, possibly, a dimerized state with a rather small order parameter

Thermodynamics of the quantum easy-plane antiferromagnet on the triangular lattice

The classical XXZ triangular-lattice antiferromagnet (TAF) shows both an Ising and a BKT transition, related to the chirality and the in-plane spin components, respectively. In this paper the quantum effects on the thermodynamic quantities are evaluated by means of the pure-quantum self-consistent harmonic approximation (PQSCHA), that allows one to deal with any spin value through classical MC simulations. We report the internal energy, the specific heat, and the in-plane correlation length of the quantum XX0 TAF, for S=1/2, 1, 5/2. The quantum transition temperatures turn out to be smaller the smaller the spin, and agree with the few available theoretical and numerical estimates.Comment: 4 pages,3 postscript figure

Real-Time Risk Management: An AAD-PDE Approach

We apply adjoint algorithmic differentiation (AAD) to the risk management of securities when their price dynamics are given by partial differential equations (PDE). We show how AAD can be applied to forward and backward PDEs in a straightforward manner. In the context of one-factor models for interest rates or default intensities, we show how price sensitivities are computed reliably and orders of magnitude faster than with a standard finite-difference approach. This significantly increased efficiency is obtained by combining (i) the adjoint forward PDE for calibrating model parameters, (ii) the adjoint backward PDE for derivatives pricing, and (iii) the implicit function theorem to avoid iterating the calibration procedure

Spiral order by disorder and lattice nematic order in a frustrated Heisenberg antiferromagnet on the honeycomb lattice

Motivated by recent experiments on Bi$_3$Mn$_4$O$_{12}$(NO$_3$), we study a frustrated $J_1$-$J_2$ Heisenberg model on the two dimensional (2D) honeycomb lattice. The classical $J_1$-$J_2$ Heisenberg model on the two dimensional (2D) honeycomb lattice has N\'eel order for $J_2 J_1/6$, it exhibits a one-parameter family of degenerate incommensurate spin spiral ground states where the spiral wave vector can point in any direction. Spin wave fluctuations at leading order lift this accidental degeneracy in favor of specific wave vectors, leading to spiral order by disorder. For spin $S=1/2$, quantum fluctuations are, however, likely to be strong enough to melt the spiral order parameter over a wide range of $J_2/J_1$. Over a part of this range, we argue that the resulting state is a valence bond solid (VBS) with staggered dimer order - this VBS is a nematic which breaks lattice rotational symmetry. Our arguments are supported by comparing the spin wave energy with the energy of the dimer solid obtained using a bond operator formalism. Turning to the effect of thermal fluctuations on the spiral ordered state, any nonzero temperature destroys the magnetic order, but the discrete rotational symmetry of the lattice remains broken resulting in a thermal analogue of the nematic VBS. We present arguments, supported by classical Monte Carlo simulations, that this nematic transforms into the high temperature symmetric paramagnet via a thermal phase transition which is in the universality class of the classical 3-state Potts (clock) model in 2D. We discuss the possible relevance of our results for honeycomb magnets, such as Bi$_3$M$_4$O$_{12}$(NO$_3$) (with M=Mn,V,Cr), and bilayer triangular lattice magnets.Comment: Slightly revise