Case study of the uniaxial anisotropic spin-1 bilinear-biquadratic Heisenberg model on a triangular lattice

We study the spin-1 model on a triangular lattice in the presence of a uniaxial anisotropy field using a cluster mean-field (CMF) approach. The interplay among antiferromagnetic exchange, lattice geometry, and anisotropy forces Gutzwiller mean-field approaches to fail in a certain region of the phase diagram. There, the CMF method yields two supersolid phases compatible with those present in the spin-1/2 XXZ model onto which the spin-1 system maps. Between these two supersolid phases, the three-sublattice order is broken and the results of the CMF approach depend heavily on the geometry and size of the cluster. We discuss the possible presence of a spin liquid in this region.

Current-driven magnetic rearrangements in spin-polarized point contacts

A method for investigating the dynamics of atomic magnetic moments in current-carrying magnetic point contacts under bias is presented. This combines the nonequilibrium Green's function (NEGF) method for evaluating the current and the charge density with a description of the dynamics of the magnetization in terms of quasistatic thermally activated transitions between stationary configurations. This method is then implemented in a tight-binding (TB) model with parameters chosen to simulate the main features of the electronic structures of magnetic transition metals. We investigate the domain wall (DW) migration in magnetic monoatomic chains sandwiched between magnetic leads, and for realistic parameters find that collinear arrangement of the magnetic moments of the chain is always favorable. Several stationary magnetic configurations are identified, corresponding to a different number of Bloch walls in the chain and to a different current. The relative stability of these configurations depends on the geometrical details of the junction and on the bias; however, we predict transitions between different configurations with activation barriers of the order of a few tens of meV. Since different magnetic configurations are associated with different resistances, this suggests an intrinsic random telegraph noise at microwave frequencies in the I-V curves of magnetic atomic point contacts at room temperature. Finally, we investigate whether or not current-induced torques are conservative.

Spontaneous spin polarization and electron localization in constrained geometries: The Wigner transition in nanowires

The Wigner transition in a jellium model of cylindrical nanowires has been investigated by density-functional computations using the local spin-density approximation. A wide range of background densities rho(b) has been explored from the nearly ideal metallic regime (r(s)=[3/4 pi rho(b)](1/3)=1) to the high correlation limit (r(s)=100). Computations have been performed using an unconstrained plane wave expansion for the Kohn-Sham orbitals and a large simulation cell with up to 480 electrons. The electron and spin distributions retain the cylindrical symmetry of the Hamiltonian at high density, while electron localization and spin polarization arise nearly simultaneously in low-density wires (r(s)similar to 30). At sufficiently low density (r(s)>= 40), the ground-state electron distribution is the superposition of well defined and nearly disjoint droplets, whose charge and spin densities integrate almost exactly to one electron and 1/2 mu(B), respectively. Droplets are arranged on radial shells and define a distorted lattice whose structure is intermediate between bcc and fcc. Dislocations and grain boundaries are apparent in the droplets' configuration found by our simulations. Our computations aim at modeling the behavior of experimental low-carried density systems made of lightly doped semiconductor nanostructures or conducting polymers.

Nonlinear theory of solitary waves associated with longitudinal particle motion in lattices - Application to longitudinal grain oscillations in a dust crystal

The nonlinear aspects of longitudinal motion of interacting point masses in a lattice are revisited, with emphasis on the paradigm of charged dust grains in a dusty plasma (DP) crystal. Different types of localized excitations, predicted by nonlinear wave theories, are reviewed and conditions for their occurrence (and characteristics) in DP crystals are discussed. Making use of a general formulation, allowing for an arbitrary (e.g. the Debye electrostatic or else) analytic potential form phi(r) and arbitrarily long site-to-site range of interactions, it is shown that dust-crystals support nonlinear kink-shaped localized excitations propagating at velocities above the characteristic DP lattice sound speed v(0). Both compressive and rarefactive kink-type excitations are predicted, depending on the physical parameter values, which represent pulse- (shock-)like coherent structures for the dust grain relative displacement. Furthermore, the existence of breather-type localized oscillations, envelope-modulated wavepackets and shocks is established. The relation to previous results on atomic chains as well as to experimental results on strongly-coupled dust layers in gas discharge plasmas is discussed.

Probing magnetic order in ultracold lattice gases

A forthcoming challenge in ultracold lattice gases is the simulation of quantum magnetism. That involves both the preparation of the lattice atomic gas in the desired spin state and the probing of the state. Here we demonstrate how a probing scheme based on atom-light interfaces gives access to the order parameters of nontrivial quantum magnetic phases, allowing us to characterize univocally strongly correlated magnetic systems produced in ultracold gases. This method, which is also nondemolishing, yields spatially resolved spin correlations and can be applied to bosons or fermions. As a proof of principle, we apply this method to detect the complete phase diagram displayed by a chain of (rotationally invariant) spin-1 bosons.

Phase diagram of spin-1 bosons on one-dimensional lattices

Spinor Bose condensates loaded in optical lattices have a rich phase diagram characterized by different magnetic order. Here we apply the density matrix renormalization group to accurately determine the phase diagram for spin-1 bosons loaded on a one-dimensional lattice. The Mott lobes present an even or odd asymmetry associated to the boson filling. We show that for odd fillings the insulating phase is always in a dimerized state. The results obtained in this work are also relevant for the determination of the ground state phase diagram of the S=1 Heisenberg model with biquadratic interaction.

Nonlinear excitations in ferromagnetic chains with nearest and next nearest neighbor exchange interactions

Weakly nonlinear excitations in one-dimensional isotropic Heisenberg ferromagnetic chains with nearest- and next-nearest-neighbor exchange interactions are considered. Based on the properties of modulational stability of corresponding linear spin waves, the existence regions of bright and dark magnetic solitons of the system are discussed in the whole Brillouin zone. The antidark soliton mode which is convex soliton super-imposed with a plane wave component is obtained near the zero-dispersion points of the spin wave frequency spectrum. The analytical results are checked by numerical simulations. [S0163;1829(98)01838-4].

Detection of Entanglement in Ultracold Lattice Gases

Probing non trivial magnetic ordering in quantum magnets realized with ultracold lattice gases demands detection methods with some spatial resolution built on it. Here we demonstrate that the Faraday matter-light interface provides an experimentally feasible tool to distinguish indubitably different quantum phases of a given many-body system in a non-demolishing way. We illustrate our approach by focussing on the Heisenberg chain for spin-1 bosons in the presence of a SU(2) symmetry breaking field. We explain how using the light signal obtained via homodyne detection one can reconstruct the phase diagram of the model. Further we show that the very same technique that provides a direct experimentally measurable signal of different order parameters can be extended to detect also the presence of multipartite entanglement in such systems.

Lorentz TEM investigation of the switching behavior in artificial spin ice building blocks

Frustration – the inability to simultaneously satisfy all interactions – occurs in a wide range of systems including neural networks, water ice and magnetic systems. An example of the latter is the so called spin-ice in pyrochlore materials [1] which have attracted a lot of interest not least due to the emergence of magnetic monopole defects when the ‘ice rules’ governing the local ordering breaks down [2]. However it is not possible to directly measure the frustrated property – the direction of the magnetic moments – in such spin ice systems with current experimental techniques. This problem can be solved by instead studying artificial spin-ice systems where the molecular magnetic moments are replaced by nanoscale ferromagnetic islands [3-8]. Two different arrangements of the ferromagnetic islands have been shown to exhibit spin ice behaviour: a square lattice maintaining four moments at each vertex [3,8] and the Kagome lattice which has only three moments per vertex but equivalent interactions between them [4-7]. Magnetic monopole defects have been observed in both types of lattices [7-8]. One of the challenges when studying these artificial spin-ice systems is that it is difficult to arrive at the fully demagnetised ground-state [6-8].
Here we present a study of the switching behaviour of building blocks of the Kagome lattice influenced by the termination of the lattice. Ferromagnetic islands of nominal size 1000 nm by 100 nm were fabricated in five island blocks using electron-beam lithography and lift-off techniques of evaporated 18 nm Permalloy (Ni80Fe20) films. Each block consists of a central island with four arms terminated by a different number and placement of ‘injection pads’, see Figure 1. The islands are single domain and magnetised along their long axis. The structures were grown on a 50 nm thick electron transparent silicon nitride membrane to allow TEM observation, which was back-coated with a 5 nm film of Au to prevent charge build-up during the TEM experiments.
To study the switching behaviour the sample was subjected to a magnetic field strong enough to magnetise all the blocks in one direction, see Figure 1. Each block obeys the Kagome lattice ‘ice-rules’ of “2-in, 1-out” or “1-in, 2-out” in this fully magnetised state. Fresnel mode Lorentz TEM images of the sample were then recorded as a magnetic field of increasing magnitude was applied in the opposite direction. While the Fresnel mode is normally used to image magnetic domain structures [9] for these types of samples it is possible to deduce the direction of the magnetisation from the Lorentz contrast [5]. All images were recorded at the same over-focus judged to give good Lorentz contrast.
The magnetisation was found to switch at different magnitudes of the applied field for nominally identical blocks. However, trends could still be identified: all the blocks with any injection pads, regardless of placement and number, switched the direction of the magnetisation of their central island at significantly smaller magnitudes of the applied magnetic field than the blocks without injection pads. It can therefore be concluded that the addition of an injection pad lowers the energy barrier to switching the connected island, acting as a nucleation site for monopole defects. In these five island blocks the defects immediately propagate through to the other side, but in a larger lattice the monopoles could potentially become trapped at a vertex and observed [10].
References

[1] M J Harris et al, Phys Rev Lett 79 (1997) p.2554.
[2] C Castelnovo, R Moessner and S L Sondhi, Nature 451 (2008) p. 42.
[3] R F Wang et al, Nature 439 (2006) 303.
[4] M Tanaka et al, Phys Rev B 73 (2006) 052411.
[5] Y Qi, T Brintlinger and J Cumings, Phys Rev B 77 (2008) 094418.
[6] E Mengotti et al, Phys Rev B 78 (2008) 144402.
[7] S Ladak et al, Nature Phys 6 (2010) 359.
[8] C Phatak et al, Phys Rev B 83 (2011) 174431.
[9] J N Chapman, J Phys D 17 (1984) 623.
[10] The authors gratefully acknowledge funding from the EPSRC under grant number EP/D063329/1.

Distillable entanglement and area laws in spin and harmonic-oscillator systems

We address the presence of nondistillable (bound) entanglement in natural many-body systems. In particular, we consider standard harmonic and spin-1/2 chains, at thermal equilibrium and characterized by few interaction parameters. The existence of bound entanglement is addressed by calculating explicitly the negativity of entanglement for different partitions. This allows us to individuate a range of temperatures for which no entanglement can be distilled by means of local operations, despite the system being globally entangled. We discuss how the appearance of bound entanglement can be linked to entanglement-area laws, typical of these systems. Various types of interactions are explored, showing that the presence of bound entanglement is an intrinsic feature of these systems. In the harmonic case, we analytically prove that thermal bound entanglement persists for systems composed by an arbitrary number of particles. Our results strongly suggest the existence of bound entangled states in the macroscopic limit also for spin-1/2 systems.

Quantum annealing of an ising spin-glass by Green's function Monte Carlo

We present an implementation of quantum annealing (QA) via lattice Green's function Monte Carlo (GFMC), focusing on its application to the Ising spin glass in transverse field. In particular, we study whether or not such a method is more effective than the path-integral Monte Carlo- (PIMC) based QA, as well as classical simulated annealing (CA), previously tested on the same optimization problem. We identify the issue of importance sampling, i.e., the necessity of possessing reasonably good (variational) trial wave functions, as the key point of the algorithm. We performed GFMC-QA runs using such a Boltzmann-type trial wave function, finding results for the residual energies that are qualitatively similar to those of CA (but at a much larger computational cost), and definitely worse than PIMC-QA. We conclude that, at present, without a serious effort in constructing reliable importance sampling variational wave functions for a quantum glass, GFMC-QA is not a true competitor of PIMC-QA.

Entanglement properties of spin models in triangular lattices

The different quantum phases appearing in strongly correlated systems as well as their transitions are closely related to the entanglement shared between their constituents. In 1D systems, it is well established that the entanglement spectrum is linked to the symmetries that protect the different quantum phases. This relation extends even further at the phase transitions where a direct link associates the entanglement spectrum to the conformal field theory describing the former. For 2D systems much less is known. The lattice geometry becomes a crucial aspect to consider when studying entanglement and phase transitions. Here, we analyze the entanglement properties of triangular spin lattice models by also considering concepts borrowed from quantum information theory such as geometric entanglement.

Thermometry precision in strongly correlated ultracold lattice gases

The precise knowledge of the temperature of an ultracold lattice gas simulating a strongly correlated
system is a question of both fundamental and technological importance. Here, we address such
question by combining tools from quantum metrology together with the study of the quantum
correlations embedded in the system at finite temperatures. Within this frame we examine the spin-
1 2 XY chain, first estimating, by means of the quantum Fisher information, the lowest attainable
bound on the temperature precision. We then address the estimation of the temperature of the sample
from the analysis of correlations using a quantum non demolishing Faraday spectroscopy method.
Remarkably, our results show that the collective quantum correlations can become optimal
observables to accurately estimate the temperature of our model in a given range of temperatures.

The Rossiter-McLaughlin effect reloaded: Probing the 3D spin-orbit geometry, differential stellar rotation, and the spatially-resolved stellar spectrum of star-planet systems

When a planet transits its host star, it blocks regions of the stellar surface from view; this causes a distortion of the spectral lines and a change in the line-of-sight (LOS) velocities, known as the Rossiter-McLaughlin (RM) effect. Since the LOS velocities depend, in part, on the stellar rotation, the RM waveform is sensitive to the star-planet alignment (which provides information on the system’s dynamical history). We present a new RM modelling technique that directly measures the spatially-resolved stellar spectrum behind the planet. This is done by scaling the continuum flux of the (HARPS) spectra by the transit light curve, and then subtracting the infrom the out-of-transit spectra to isolate the starlight behind the planet. This technique does not assume any shape for the intrinsic local profiles. In it, we also allow for differential stellar rotation and centre-to-limb variations in the convective blueshift. We apply this technique to HD 189733 and compare to 3D magnetohydrodynamic (MHD) simulations. We reject rigid body rotation with high confidence (>99% probability), which allows us to determine the occulted stellar latitudes and measure the stellar inclination. In turn, we determine both the sky-projected (λ ≈ −0.4 ± 0.2◦) and true 3D obliquity (ψ ≈ 7+12 −4 ◦ ). We also find good agreement with the MHD simulations, with no significant centre-to-limb variations detectable in the local profiles. Hence, this technique provides a new powerful tool that can probe stellar photospheres, differential rotation, determine 3D obliquities, and remove sky-projection biases in planet migration theories. This technique can be implemented with existing instrumentation, but will become even more powerful with the next generation of high-precision radial velocity spectrographs.