-
Competition between Neel, Haldane nematic, plaquette valence bond solid, and $(π,π)$ valence bond solid phases in SU(N) analogs of $S=1$ square-lattice antiferromagnets
Authors:
Souvik Kundu,
Nisheeta Desai,
Kedar Damle
Abstract:
We use stochastic series expansion (SSE) quantum Monte Carlo (QMC) methods to study the phases and transitions displayed by a class of sign-free designer Hamiltonians for SU($N$) analogs of spin $S=1$ quantum antiferromagnets on the square lattice. The SU($N$) spins are generators of the single-row two-column representation (complex conjugate of single-row two-column representation) on $A$ ($B$) s…
▽ More
We use stochastic series expansion (SSE) quantum Monte Carlo (QMC) methods to study the phases and transitions displayed by a class of sign-free designer Hamiltonians for SU($N$) analogs of spin $S=1$ quantum antiferromagnets on the square lattice. The SU($N$) spins are generators of the single-row two-column representation (complex conjugate of single-row two-column representation) on $A$ ($B$) sublattices of the square lattice, and the Hamiltonian is designed to explore the competition between the nearest neighbour antiferromagnetic exchange couplings $J$ and four-spin interactions $Q$ that favor a plaquette-ordered valence bond solid (p-VBS) ground state. We find that this state is indeed established at large $Q/J$ for all $N > 3$. For $3< N \leq 9$, the ground state exhibits a direct first order quantum phase transition from a small-$Q/J$ Néel ordered antiferromagnetic state to this large-$Q/J$ p-VBS state. The ground state at $Q/J=0$ for $N \geq 10$ has been previously reported to be a valence bond nematic state, dubbed the Haldane nematic in recent literature. For small nonzero $Q/J$ and $N \geq 10$, we additionally find an unusual intermediate state in which the bond energy has a Bragg peak at wavevector $(π,π)$ with no accompanying Bragg peaks at wavevectors $(π,0)$ and $(0, π)$. This $(π, π)$ state also appears to be metastable at $Q/J =0$, as evidenced by low temperature histograms of the Haldane nematic and $(π, π)$ order parameters. Deep in the p-VBS phase of the ground state phase diagram, we find the temperature-driven melting of the p-VBS order is in the Ashkin-Teller universality class. In this regime, we identify an interesting signature of Ashkin-Teller criticality in bond correlations at wavevector $(π, π)$; this is in addition to the expected critical fluctuations of the conventional p-VBS order parameter at wavevectors $(π,0)$ and $(0,π)$.
△ Less
Submitted 21 September, 2023;
originally announced September 2023.
-
Deconfined pseudocriticality in a model spin-1 quantum antiferromagnet
Authors:
Vikas Vijigiri,
Sumiran Pujari,
Nisheeta Desai
Abstract:
Berry phase interference arguments that underlie the theory of deconfined quantum criticality (DQC) for spin-$\frac{1}{2}$ antiferromagnets have also been invoked to allow for continuous transitions in spin-1 magnets including a Néel to (columnar) valence bond solid (cVBS) transition. We provide a microscopic model realization of this transition on the square lattice consisting of Heisenberg excha…
▽ More
Berry phase interference arguments that underlie the theory of deconfined quantum criticality (DQC) for spin-$\frac{1}{2}$ antiferromagnets have also been invoked to allow for continuous transitions in spin-1 magnets including a Néel to (columnar) valence bond solid (cVBS) transition. We provide a microscopic model realization of this transition on the square lattice consisting of Heisenberg exchange ($J_H$) and biquadratic exchange ($J_B$) that favor a Néel phase, and a designed $Q$-term ($Q_B$) interaction which favors a cVBS through large-scale quantum Monte Carlo (QMC) simulations. For $J_H=0$, this model is equivalent to the $SU(3)$ $JQ$ model with a Néel-cVBS transition that has been argued to be DQC through QMC. Upon turning on $J_H$ which brings down the symmetry to $SU(2)$, we find multiple signatures -- a single critical point, high quality collapse of correlation ratios and order parameters, "$U(1)$-symmetric" cVBS histograms and lack of double-peak in order parameter histograms for largest sizes studied near the critical point -- that are highly suggestive of a continuous transition scenario. However, Binder analysis finds negative dips that grow sub-extensively that we interpret as these transitions rather being pseudocritical. This along with recent results on spin-$\frac{1}{2}$ models suggests that deconfined pseudocriticality is the more generic scenario.
△ Less
Submitted 13 June, 2024; v1 submitted 26 July, 2023;
originally announced July 2023.
-
Notes on resummation-based quantum Monte Carlo vis-à-vis sign-problematic Heisenberg models on canonical geometrically frustrated lattices
Authors:
Nisheeta Desai,
Sumiran Pujari
Abstract:
We show here that a direct application of resummation-based quantum Monte Carlo (QMC) -- implemented recently for sign-problem-free SU(2)-symmetric Hamiltonians in the stochastic series expansion (SSE) framework -- does not reduce the sign problem for frustrated SU(2)-symmetric Heisenberg antiferromagnets on canonical geometrically frustrated lattices composed of triangular motifs such as the tria…
▽ More
We show here that a direct application of resummation-based quantum Monte Carlo (QMC) -- implemented recently for sign-problem-free SU(2)-symmetric Hamiltonians in the stochastic series expansion (SSE) framework -- does not reduce the sign problem for frustrated SU(2)-symmetric Heisenberg antiferromagnets on canonical geometrically frustrated lattices composed of triangular motifs such as the triangular lattice. In the process, we demonstrate that resummation-based updates do provide an ergodic sampling of the SSE-based QMC configurations which can be an issue when using the standard SSE updates, however, severely limited by the sign problem as previously mentioned. The notions laid out in these notes may be useful in the design of better algorithms for geometrically frustrated magnets.
△ Less
Submitted 20 June, 2022;
originally announced June 2022.
-
Instability and self-propulsion of active droplets along a wall
Authors:
Nikhil Desai,
Sebastien Michelin
Abstract:
Active droplets can swim spontaneously in viscous flows as a result of the non-linear convective transport of a chemical solute produced at their surface by the Marangoni and/or phoretic flows generated by this solute's inhomogeneous distribution, provided the ratio of convective-to-diffusive solute transport, or Péclet number $\mbox{Pe}$ is large enough. As the result of their net buoyancy, activ…
▽ More
Active droplets can swim spontaneously in viscous flows as a result of the non-linear convective transport of a chemical solute produced at their surface by the Marangoni and/or phoretic flows generated by this solute's inhomogeneous distribution, provided the ratio of convective-to-diffusive solute transport, or Péclet number $\mbox{Pe}$ is large enough. As the result of their net buoyancy, active drops typically evolve at a small finite distance $d$ from rigid boundaries. Yet, existing models systematically focus on unbounded flows, ignoring the effect of the wall proximity on the intrisically-nonlinear nature of their propulsion mechanism. In contrast, we obtain here a critical insight on the propulsion of active drops near walls by analysing their stability to non-axisymmetric perturbations and the resulting emergence of self-propulsion along the wall with no limiting assumption on the wall-distance $d$. Dipolar or quadrupolar axisymmetric (levitating) base states are identified depending on $d$ and $\mbox{Pe}$. Perhaps counter-intuitively, a reduction in the drop-wall separation $d$ is observed to destabilize these modes and to promote self-propulsion, as a result from the confinement-induced localisation of the chemical gradients driving the motion. In addition, quadrupolar states are more unstable than their dipolar counterparts due to the redistribution of the chemical perturbation by the base flow, favouring the emergence of stronger slip forcing on the drop surface.
△ Less
Submitted 2 November, 2021;
originally announced November 2021.
-
Melting of three-sublattice order in triangular lattice Ising antiferromagnets: Power-law order, $Z_6$ parafermionic multicriticality, and weakly first order transitions
Authors:
G. Rakala,
N. Desai,
S. Shivam,
K. Damle
Abstract:
The nature of the thermal melting process by which triangular-lattice Ising antiferromagnets lose their low-temperature ferrimagnetic three-sublattice order depends on the range of the interactions: It changes character when second and third neighbour ferromagnetic interactions become comparable to the nearest-neighbour antiferromagnetic coupling. We present a detailed numerical characterization o…
▽ More
The nature of the thermal melting process by which triangular-lattice Ising antiferromagnets lose their low-temperature ferrimagnetic three-sublattice order depends on the range of the interactions: It changes character when second and third neighbour ferromagnetic interactions become comparable to the nearest-neighbour antiferromagnetic coupling. We present a detailed numerical characterization of the corresponding threshold at which two-step melting of three-sublattice order gives way to a direct first-order transition at which this order is lost. The multicritical behaviour at this threshold is argued to be in the universality class of the $Z_6$ parafermion conformal field theory with central charge $c=5/4$. The presence of this multicritical threshold influences the melting behaviour and long-wavelength properties over a fairly large range of parameters, and at temperatures that are of the same order as the exchange interactions. It is therefore of potential experimental relevance in the context of easy-axis triangular lattice antiferromagnets that display such low temperature ordering.
△ Less
Submitted 7 September, 2021;
originally announced September 2021.
-
Resummation-based updates for Stochastic Series Expansion Quantum Monte Carlo
Authors:
Nisheeta Desai,
Sumiran Pujari
Abstract:
For spin rotational symmetric models with a positive-definite high-temperature expansion of the partition function, a stochastic sampling of the series expansion upon partial resummation becomes logically equivalent to sampling an uncoloured closely-packed loop-gas model in one higher dimension. Based on this, we devise quantum Monte Carlo updates that importance-sample loop configurations for gen…
▽ More
For spin rotational symmetric models with a positive-definite high-temperature expansion of the partition function, a stochastic sampling of the series expansion upon partial resummation becomes logically equivalent to sampling an uncoloured closely-packed loop-gas model in one higher dimension. Based on this, we devise quantum Monte Carlo updates that importance-sample loop configurations for general $SU(N)$ in fundamental and higher-symmetric representations. The algorithmic performance systematically improves with increase in (continuous) $N$ allowing efficient simulation of quantum paramagnets. The underlying reason for the increased efficacy is the correspondence of quantum paramagnetic phases like valence bond solids to short-loop phases on the loop-gas side rather than the particular value of $N$. This also gives a connection between Sandvik's $JQ$ model class and classical loop-gas models in the deconfined universality class.
△ Less
Submitted 10 August, 2021; v1 submitted 14 April, 2021;
originally announced April 2021.
-
Bilayer Coulomb phase of two dimensional dimer models: Absence of power-law columnar order
Authors:
Nisheeta Desai,
Sumiran Pujari,
K. Damle
Abstract:
We study the fully-packed dimer model on the bilayer square lattice with fugacity equal to $z$ ($1$) for inter-layer (intra-layer) dimers, and intra-layer interaction $V$ between neighbouring parallel dimers on any elementary plaquette in either layer. For a range of not-too-large $z> 0$ and repulsive interactions $0< V < V_s$ (with $V_s \approx 2.1$), we demonstrate the existence of a {\em bilaye…
▽ More
We study the fully-packed dimer model on the bilayer square lattice with fugacity equal to $z$ ($1$) for inter-layer (intra-layer) dimers, and intra-layer interaction $V$ between neighbouring parallel dimers on any elementary plaquette in either layer. For a range of not-too-large $z> 0$ and repulsive interactions $0< V < V_s$ (with $V_s \approx 2.1$), we demonstrate the existence of a {\em bilayer Coulomb phase} with purely dipolar two-point functions, {\em i.e.}, without the power-law columnar order that characterizes the usual Coulomb phase of square and honeycomb lattice dimer models. The transition line $z_{c}(V)$ separating this bilayer Coulomb phase from a large-$z$ disordered phase is argued to be in the inverted Kosterlitz-Thouless universality class. Additionally, we argue for the possibility of a tricritical point at which the bilayer Coulomb phase, the large-$z$ disordered phase and the large-$V$ staggered phase meet in the large-$z$, large-$V$ part of the phase diagram. In contrast, for the attractive case with $ V_{cb} < V \leq 0$ ($V_{cb} \approx -1.2$), we argue that any $z > 0$ destroys the power-law correlations of the $z=0$ decoupled layers, and leads immediately to a short-range correlated state, albeit with a slow crossover for small $|V|$. For $V_{c} < V < V_{cb}$ ($V_{c} \approx -1.55$), we predict that any small nonzero $z$ immediately gives rise to long-range {\em bilayer columnar order} although the $z=0$ decoupled layers remain power-law correlated in this regime; this implies a non-monotonic $z$ dependence of the columnar order parameter for fixed $V$ in this regime. This bilayer columnar ordered state is separated from the large-$z$ disordered state by a line of Ashkin-Teller transitions $z_{\rm AT}(V)$.
△ Less
Submitted 29 April, 2021; v1 submitted 9 November, 2020;
originally announced November 2020.
-
Duality and domain wall dynamics in a twisted Kitaev chain
Authors:
C. M. Morris,
Nisheeta Desai,
J. Viirok,
D. Huvonen,
U. Nagel,
T. Room,
J. W. Krizan,
R. J. Cava,
T. M. McQueen,
S. M. Koohpayeh,
Ribhu K. Kaul,
N. P. Armitage
Abstract:
The Ising chain in transverse field is a paradigmatic model for a host of physical phenomena, including spontaneous symmetry breaking, topological defects, quantum criticality, and duality. Although the quasi-1D ferromagnet CoNb$_2$O$_6$ has been put forward as the best material example of the transverse field Ising model, it exhibits significant deviations from ideality. Through a combination of…
▽ More
The Ising chain in transverse field is a paradigmatic model for a host of physical phenomena, including spontaneous symmetry breaking, topological defects, quantum criticality, and duality. Although the quasi-1D ferromagnet CoNb$_2$O$_6$ has been put forward as the best material example of the transverse field Ising model, it exhibits significant deviations from ideality. Through a combination of THz spectroscopy and theory, we show that CoNb$_2$O$_6$ in fact is well described by a different model with strong bond dependent interactions, which we dub the {\it twisted Kitaev chain}, as these interactions share a close resemblance to a one-dimensional version of the intensely studied honeycomb Kitaev model. In this model the ferromagnetic ground state of CoNb$_2$O$_6$ arises from the compromise between two distinct alternating axes rather than a single easy axis. Due to this frustration, even at zero applied field domain-wall excitations have quantum motion that is described by the celebrated Su-Schriefer-Heeger model of polyacetylene. This leads to rich behavior as a function of field. Despite the anomalous domain wall dynamics, close to a critical transverse field the twisted Kitaev chain enters a universal regime in the Ising universality class. This is reflected by the observation that the excitation gap in CoNb$_2$O$_6$ in the ferromagnetic regime closes at a rate precisely twice that of the paramagnet. This originates in the duality between domain walls and spin-flips and the topological conservation of domain wall parity. We measure this universal ratio `2' to high accuracy -- the first direct evidence for the Kramers-Wannier duality in nature.
△ Less
Submitted 29 September, 2020;
originally announced September 2020.
-
First order phase transitions in the square lattice "easy-plane" J-Q model
Authors:
Nisheeta Desai,
Ribhu K. Kaul
Abstract:
We study the quantum phase transition between the superfluid and valence bond solid in "easy-plane" J-Q models on the square lattice. The Hamiltonian we study is a linear combination of two model Hamiltonians: (1) an SU(2) symmetric model, which is the well known J-Q model that does not show any direct signs of a discontinuous transition on the largest lattices and is presumed continuous, and (2)…
▽ More
We study the quantum phase transition between the superfluid and valence bond solid in "easy-plane" J-Q models on the square lattice. The Hamiltonian we study is a linear combination of two model Hamiltonians: (1) an SU(2) symmetric model, which is the well known J-Q model that does not show any direct signs of a discontinuous transition on the largest lattices and is presumed continuous, and (2) an easy plane version of the J-Q model, which shows clear evidence for a first order transition even on rather small lattices of size $L\approx16$. A parameter $0\leqλ\leq 1$ ($λ=0$ being the easy-plane model and $λ=1$ being the SU(2) symmetric J-Q model) allows us to smoothly interpolate between these two limiting models. We use stochastic series expansion (SSE) quantum Monte Carlo (QMC) to investigate the nature of this transition as $λ$ is varied - here we present studies for $λ=0,0.5,0.75,0.85,0.95$ and $1$. While we find that the first order transition weakens as $λ$ is increased from 0 to 1, we find no evidence that the transition becomes continuous until the SU(2) symmetric point, $λ=1$. We thus conclude that the square lattice superfluid-VBS transition in the two-component easy-plane model is generically first order.
△ Less
Submitted 29 November, 2020; v1 submitted 26 September, 2019;
originally announced September 2019.
-
Spin-$S$ designer hamiltonians and the square lattice $S=1$ Haldane nematic
Authors:
Nisheeta Desai,
Ribhu K. Kaul
Abstract:
We introduce a strategy to write down lattice models of spin rotational symmetric Hamiltonians with arbitrary spin-$S$ that are Marshall positive and can be simulated efficiently using world line Monte Carlo methods. As an application of our approach we consider a square lattice $S=1$ model for which we design a $3\times 3$ - spin plaquette interaction. By numerical simulations we establish that o…
▽ More
We introduce a strategy to write down lattice models of spin rotational symmetric Hamiltonians with arbitrary spin-$S$ that are Marshall positive and can be simulated efficiently using world line Monte Carlo methods. As an application of our approach we consider a square lattice $S=1$ model for which we design a $3\times 3$ - spin plaquette interaction. By numerical simulations we establish that our model realizes a novel "Haldane nematic" phase that breaks lattice rotational symmetry by the spontaneous formation of Haldane chains, while preserving spin rotations, time reversal and lattice translations. By supplementing our model with a two-spin Heisenberg interaction, we present a study of the transition between Néel and Haldane nematic phase, which we find to be of first order.
△ Less
Submitted 27 September, 2019; v1 submitted 21 April, 2019;
originally announced April 2019.
-
First-order Néel-cVBS transition in a model square lattice $S=1$ antiferromagnet
Authors:
Julia Wildeboer,
Nisheeta Desai,
Jonathan D'Emidio,
Ribhu K. Kaul
Abstract:
We study the Néel to four-fold columnar valence bond solid (cVBS) quantum phase transition in a sign free $S=1$ square lattice model. This is the same kind of transition that for $S=1/2$ has been argued to realize the prototypical deconfined critical point. Extensive numerical simulations of the square lattice $S=1/2$ Néel-VBS transition have found consistency with the DCP scenario with no direct…
▽ More
We study the Néel to four-fold columnar valence bond solid (cVBS) quantum phase transition in a sign free $S=1$ square lattice model. This is the same kind of transition that for $S=1/2$ has been argued to realize the prototypical deconfined critical point. Extensive numerical simulations of the square lattice $S=1/2$ Néel-VBS transition have found consistency with the DCP scenario with no direct evidence for first order behavior. In contrast to the $S=1/2$ case, in our quantum Monte Carlo simulations for the $S=1$ model, we present unambiguous evidence for a direct conventional first-order quantum phase transition. Classic signs for a first order transition demonstrating co-existence including double peaked histograms and switching behavior are observed. The sharp contrast from the $S=1/2$ case is remarkable, and is a striking demonstration of the role of the size of the quantum spin in the phase diagram of two dimensional lattice models.
△ Less
Submitted 4 September, 2019; v1 submitted 14 August, 2018;
originally announced August 2018.
-
Three Dimensional Polarimetric Neutron Tomography of Magnetic Fields
Authors:
Morten Sales,
Markus Strobl,
Takenao Shinohara,
Anton Tremsin,
Luise Theil Kuhn,
William R. B. Lionheart,
Naeem M. Desai,
Anders Bjorholm Dahl,
Søren Schmidt
Abstract:
Through the use of Time-of-Flight Three Dimensional Polarimetric Neutron Tomography (ToF 3DPNT) we have for the first time successfully demonstrated a technique capable of measuring and reconstructing three dimensional magnetic field strengths and directions unobtrusively and non-destructively with the potential to probe the interior of bulk samples which is not amenable otherwise.
Using a pione…
▽ More
Through the use of Time-of-Flight Three Dimensional Polarimetric Neutron Tomography (ToF 3DPNT) we have for the first time successfully demonstrated a technique capable of measuring and reconstructing three dimensional magnetic field strengths and directions unobtrusively and non-destructively with the potential to probe the interior of bulk samples which is not amenable otherwise.
Using a pioneering polarimetric set-up for ToF neutron instrumentation in combination with a newly developed tailored reconstruction algorithm, the magnetic field generated by a current carrying solenoid has been measured and reconstructed, thereby providing the proof-of-principle of a technique able to reveal hitherto unobtainable information on the magnetic fields in the bulk of materials and devices, due to a high degree of penetration into many materials, including metals, and the sensitivity of neutron polarisation to magnetic fields. The technique puts the potential of the ToF time structure of pulsed neutron sources to full use in order to optimise the recorded information quality and reduce measurement time.
△ Less
Submitted 2 February, 2018; v1 submitted 17 April, 2017;
originally announced April 2017.
-
TiO2 based Nanostructured Memristor for RRAM and Neuromorphic Applications: A Simulation Approach
Authors:
T. D. Dongale,
P. J. Patil,
N. K. Desai,
P. P. Chougule,
S. M. Kumbhar,
P. P. Waifalkar,
P. B. Patil,
R. S. Vhatkar,
M. V. Takale,
P. K. Gaikwad,
R. K. Kamat
Abstract:
We report simulation of nanostructured memristor device using piecewise linear and nonlinear window functions for RRAM and neuromorphic applications. The linear drift model of memristor has been exploited for the simulation purpose with the linear and non-linear window function as the mathematical and scripting basis. The results evidences that the piecewise linear window function can aptly simula…
▽ More
We report simulation of nanostructured memristor device using piecewise linear and nonlinear window functions for RRAM and neuromorphic applications. The linear drift model of memristor has been exploited for the simulation purpose with the linear and non-linear window function as the mathematical and scripting basis. The results evidences that the piecewise linear window function can aptly simulate the memristor characteristics pertaining to RRAM application. However, the nonlinear window function could exhibit the nonlinear phenomenon in simulation only at the lower magnitude of control parameter. This has motivated us to propose a new nonlinear window function for emulating the simulation model of the memristor. Interestingly, the proposed window function is scalable up to f(x)=1 and exhibits the nonlinear behavior at higher magnitude of control parameter. Moreover, the simulation results of proposed nonlinear window function are encouraging and reveals the smooth nonlinear change from LRS to HRS and vice versa and therefore useful for the neuromorphic applications.
△ Less
Submitted 25 January, 2016;
originally announced January 2016.
-
Fabrication of Ag/Tl-Ba-Ca-CuO/CdSe nanostructure by electrodeposition technique
Authors:
P. M. Shirage,
D. D. Shivagan,
L. A. Ekal,
N. V. Desai,
S. B. Mane,
S. H. Pawar
Abstract:
Junction between metal, superconductor and semiconductor materials seems to be promising for their applications as the electronic devices, microelectronic devices and high power transmission. E-mail shpawar_phy@unishivaji.ac.in.
Junction between metal, superconductor and semiconductor materials seems to be promising for their applications as the electronic devices, microelectronic devices and high power transmission. E-mail shpawar_phy@unishivaji.ac.in.
△ Less
Submitted 16 May, 2003;
originally announced May 2003.