-
Stationary particle currents in sedimenting active matter wetting a wall
Authors:
Matthieu Mangeat,
Shauri Chakraborty,
Adam Wysocki,
Heiko Rieger
Abstract:
Recently it was predicted, on the basis of a lattice gas model, that scalar active matter in a gravitational field would rise against gravity up a confining wall or inside a thin capillary - in spite of repulsive particle-wall interactions [Phys. Rev. Lett. 124, 048001 (2020)]. In this paper we confirm this prediction with sedimenting active Brownian particles (ABPs) in a box numerically and eluci…
▽ More
Recently it was predicted, on the basis of a lattice gas model, that scalar active matter in a gravitational field would rise against gravity up a confining wall or inside a thin capillary - in spite of repulsive particle-wall interactions [Phys. Rev. Lett. 124, 048001 (2020)]. In this paper we confirm this prediction with sedimenting active Brownian particles (ABPs) in a box numerically and elucidate the mechanism leading to the formation of a meniscus rising above the bulk of the sedimentation region. The height of the meniscus increases with the activity of the system, algebraically with the Péclet number. The formation of the meniscus is determined by a stationary circular particle current, a vortex, centered at the base of the meniscus, whose size and strength increase with the ABP activity. The origin of these vortices can be traced back to the confinement of the ABPs in a box: already the stationary state of ideal (non-interacting) ABPs without gravitation displays circular currents that arrange in a highly symmetric way in the eight octants of the box. Gravitation distorts this vortex configuration downward, leaving two major vortices at the two side walls, with a strong downward flow along the walls. Repulsive interactions between the ABPs change this situation only as soon as motility induced phase separation (MIPS) sets in and forms a dense, sedimented liquid region at the bottom, which pushes the center of the vortex upwards towards the liquid-gas interface. Self-propelled particles therefore represent an impressive realization of scalar active matter that forms stationary particle currents being able to perform visible work against gravity or any other external field, which we predict to be observable experimentally in active colloids under gravitation.
△ Less
Submitted 15 January, 2024; v1 submitted 18 September, 2023;
originally announced September 2023.
-
Relativistic Douglas-Kroll-Hess Calculations of Hyperfine Interactions within First Principles Multireference Methods
Authors:
Aleksander L. Wysocki,
Kyungwha Park
Abstract:
Relativistic magnetic hyperfine interaction Hamiltonian based on the Douglas-Kroll-Hess (DKH) theory up to the second order is implemented within the ab initio multireference methods including spin-orbit coupling in the Molcas/OpenMolcas package. This implementation is applied to calculate relativistic hyperfine coupling (HFC) parameters for atomic systems and diatomic radicals with valence s or d…
▽ More
Relativistic magnetic hyperfine interaction Hamiltonian based on the Douglas-Kroll-Hess (DKH) theory up to the second order is implemented within the ab initio multireference methods including spin-orbit coupling in the Molcas/OpenMolcas package. This implementation is applied to calculate relativistic hyperfine coupling (HFC) parameters for atomic systems and diatomic radicals with valence s or d orbitals by systematically varying active space size in the restricted active space self-consistent field (RASSCF) formalism with restricted active space state interaction (RASSI) for spin-orbit coupling. The DKH relativistic treatment of the hyperfine interaction reduces the Fermi contact contribution to the HFC due to the presence of kinetic factors that regularize the singularity of the Dirac delta function in the nonrelativitic Fermi contact operator. This effect is more prominent for heavier nuclei. As the active space size increases, the relativistic correction of the Fermi contact contribution converges well to the experimental data for light and moderately heavy nuclei. The relativistic correction, however, does not significantly affect the spin-dipole contribution to the hyperfine interaction. In addition to the atomic and molecular systems, the implementation is applied to calculate the relativistic HFC parameters for large trivalent and divalent Tb-based single-molecule magnets (SMMs) such as Tb(III)Pc$_2$ and Tb(II)(Cp$^\text{iPr5}$)$_2$ without ligand truncation using well-converged basis sets. In particular, for the divalent SMM which has an unpaired valence 6s/5d hybrid orbital, the relativistic treatment of HFC is crucial for a proper description of the Fermi contact contribution. Even with the relativistic hyperfine Hamiltonian, the divalent SMM is shown to exhibit strong tunability of HFC via an external electric field (i.e., strong hyperfine Stark effect).
△ Less
Submitted 17 September, 2023;
originally announced September 2023.
-
How to steer active colloids up a vertical wall
Authors:
Adérito Fins Carreira,
Adam Wysocki,
Christophe Ybert,
Mathieu Leocmach,
Heiko Rieger,
Cécile Cottin-Bizonne
Abstract:
An important challenge in active matter lies in harnessing useful global work from entities that produce work locally, e.g., via self-propulsion. We investigate here the active matter version of a classical capillary rise effect, by considering a non-phase separated sediment of self-propelled Janus colloids in contact with a vertical wall. We provide experimental evidence of an unexpected and dyna…
▽ More
An important challenge in active matter lies in harnessing useful global work from entities that produce work locally, e.g., via self-propulsion. We investigate here the active matter version of a classical capillary rise effect, by considering a non-phase separated sediment of self-propelled Janus colloids in contact with a vertical wall. We provide experimental evidence of an unexpected and dynamic adsorption layer at the wall. Additionally, we develop a complementary numerical model that recapitulates the experimental observations. We show that an adhesive and aligning wall enhances the pre-existing polarity heterogeneity within the bulk, enabling polar active particles to climb up a wall against gravity, effectively powering a global flux. Such steady-state flux has no equivalent in a passive wetting layer.
△ Less
Submitted 6 July, 2023;
originally announced July 2023.
-
Computational Insights into Electronic Excitations, Spin-Orbit Coupling Effects, and Spin Decoherence in Cr(IV)-based Molecular Qubits
Authors:
Karolina Janicka,
Aleksander L. Wysocki,
Kyungwha Park
Abstract:
The great success of point defects and dopants in semiconductors for quantum information processing has invigorated a search for molecules with analogous properties. Flexibility and tunability of desired properties in a large chemical space have great advantages over solid-state systems. The properties analogous to point defects were demonstrated in Cr(IV)-based molecular family, Cr(IV)(aryl)$_4$,…
▽ More
The great success of point defects and dopants in semiconductors for quantum information processing has invigorated a search for molecules with analogous properties. Flexibility and tunability of desired properties in a large chemical space have great advantages over solid-state systems. The properties analogous to point defects were demonstrated in Cr(IV)-based molecular family, Cr(IV)(aryl)$_4$, where the electronic spin states were optically initialized, read out, and controlled. Despite this kick-start, there is still a large room for enhancing properties crucial for molecular qubits. Here we provide computational insights into key properties of the Cr(IV)-based molecules aimed at assisting chemical design of efficient molecular qubits. Using the multireference ab-initio methods, we investigate the electronic states of Cr(IV)(aryl)$_4$ molecules with slightly different ligands, showing that the zero-phonon line energies agree with the experiment, and that the excited spin-triplet and spin-singlet states are highly sensitive to small chemical perturbations. By adding spin-orbit interaction, we find that the sign of the uniaxial zero-field splitting (ZFS) parameter is negative for all considered molecules, and discuss optically-induced spin initialization via non-radiative intersystem crossing. We quantify (super)hyperfine coupling to the $^{53}$Cr nuclear spin and to the $^{13}$C and $^1$H nuclear spins, and we discuss electron spin decoherence. We show that the splitting or broadening of the electronic spin sub-levels due to superhyperfine interaction with $^1$H nuclear spins decreases by an order of magnitude when the molecules have a substantial transverse ZFS parameter.
△ Less
Submitted 13 October, 2022; v1 submitted 30 April, 2022;
originally announced May 2022.
-
Interacting particles in an activity landscape
Authors:
Adam Wysocki,
Anil K. Dasanna,
Heiko Rieger
Abstract:
We study interacting active Brownian particles (ABPs) with a space-dependent swim velocity via simulation and theory. We find that, although an equation of state exists, a mechanical equilibrium does not apply to ABPs in activity landscapes. The pressure difference originates in the flux of polar order and the gradient of swim velocity across the interface between regions of different activity. In…
▽ More
We study interacting active Brownian particles (ABPs) with a space-dependent swim velocity via simulation and theory. We find that, although an equation of state exists, a mechanical equilibrium does not apply to ABPs in activity landscapes. The pressure difference originates in the flux of polar order and the gradient of swim velocity across the interface between regions of different activity. In contrast to motility-induced phase separation of ABPs with a homogeneous swim velocity, a critical point does not exist for an active-passive patch system, which continuously splits into a dense and a dilute phase with increasing activity. However, if the global density is so high that not all particles can be packed onto the inactive patch, then MIPS-like behavior is restored and the pressure is balanced again.
△ Less
Submitted 3 April, 2022;
originally announced April 2022.
-
Striped Patterns in Radially Driven Suspensions with Open Boundaries
Authors:
Mahdieh Mohammadi,
Maniya Maleki,
Adam Wysocki,
M. Reza Shaebani
Abstract:
We study the motion of radially driven fluid-immersed particles in a novel Hele-Shaw cell with open boundaries. The initially uniform suspension forms a striped pattern within a specific range of horizontal oscillation frequencies and for sufficiently large amplitudes. We observe that the initial coarsening dynamics of the stripes gradually slows down and the pattern reaches a steady state after a…
▽ More
We study the motion of radially driven fluid-immersed particles in a novel Hele-Shaw cell with open boundaries. The initially uniform suspension forms a striped pattern within a specific range of horizontal oscillation frequencies and for sufficiently large amplitudes. We observe that the initial coarsening dynamics of the stripes gradually slows down and the pattern reaches a steady state after a few minutes. The distance between the stripes in the steady state exhibits an exponentially saturating increase with increased oscillation amplitude or frequency. The width of the stripes decreases as a power-law with the frequency while its amplitude dependence follows a logistic function. We propose a mechanism -- based on the interplay between shear stress, hydrodynamic interactions, and frictional forces -- to link the structural characteristics of the stripes to the properties of the oscillatory external drive.
△ Less
Submitted 30 October, 2021; v1 submitted 14 July, 2021;
originally announced July 2021.
-
A multiconfigurational study of the negatively charged nitrogen-vacancy center in diamond
Authors:
Churna Bhandari,
Aleksander L. Wysocki,
Sophia E. Economou,
Pratibha Dev,
Kyungwha Park
Abstract:
Deep defects in wide band gap semiconductors have emerged as leading qubit candidates for realizing quantum sensing and information applications. Due to the spatial localization of the defect states, these deep defects can be considered as artificial atoms/molecules in a solid state matrix. Here we show that unlike single-particle treatments, the multiconfigurational quantum chemistry methods, tra…
▽ More
Deep defects in wide band gap semiconductors have emerged as leading qubit candidates for realizing quantum sensing and information applications. Due to the spatial localization of the defect states, these deep defects can be considered as artificial atoms/molecules in a solid state matrix. Here we show that unlike single-particle treatments, the multiconfigurational quantum chemistry methods, traditionally reserved for atoms/molecules, accurately describe the many-body characteristics of the electronic states of these defect centers and correctly predict properties that single-particle treatments fail to obtain. We choose the negatively charged nitrogen-vacancy (NV$^-$) center in diamond as the prototype defect to study with these techniques due to its importance for quantum information applications and because its properties are well-known, which makes it an ideal benchmark system. By properly accounting for electron correlations and including spin-orbit coupling and dipolar spin-spin coupling in the quantum chemistry calculations, for the NV$^-$ center in diamond clusters, we are able to: (i) show the correct splitting of the ground (first-excited) triplet state into two levels (four levels), (ii) calculate zero-field splitting values of the ground and excited triplet states, in good agreement with experiment, and (iii) calculate the energy differences between ground and exited spin-triplet and spin-singlet states, as well as their ordering, which are also found to be in good agreement with recent experimental data. The numerical procedure we have developed is general and it can screen other color centers whose properties are not well known but promising for applications.
△ Less
Submitted 11 October, 2020; v1 submitted 23 August, 2020;
originally announced August 2020.
-
Electrically tuned hyperfine spectrum in neutral Tb(II)(Cp$^{\rm{iPr5}}$)$_2$ single-molecule magnet
Authors:
Robert L. Smith,
Aleksander L. Wysocki,
Kyungwha Park
Abstract:
Molecular spin qubits with long spin coherence time as well as non-invasive operation methods on such qubits are in high demand. It was shown that both molecular electronic and nuclear spin levels can be used as qubits. In solid state systems with dopants, an electric field was shown to effectively change the spacing between the nuclear spin qubit levels when the electron spin density is high at t…
▽ More
Molecular spin qubits with long spin coherence time as well as non-invasive operation methods on such qubits are in high demand. It was shown that both molecular electronic and nuclear spin levels can be used as qubits. In solid state systems with dopants, an electric field was shown to effectively change the spacing between the nuclear spin qubit levels when the electron spin density is high at the nucleus of the dopant. Inspired by such solid-state systems, we propose that divalent lanthanide (Ln) complexes with an unusual electronic configuration of Ln$^{2+}$ have a strong interaction between the Ln nuclear spin and the electronic degrees of freedom, which renders electrical tuning of the interaction. As an example, we study electronic structure and hyperfine interaction of the $^{159}$Tb nucleus in a neutral Tb(II)(Cp$^{\rm{iPr5}}$)$_2$ single-molecule magnet (SMM) using the complete active space self-consistent field method with spin-orbit interaction included within the restricted active space state interaction. Our calculations show that the low-energy states arise from $4f^8(6s,5d_{z^2})^1$, 4$f^8$(5$d_{x^2-y^2}$)$^1$, and $4f^8(5d_{xy})^1$ configurations. We compute the hyperfine interaction parameters and the electronic-nuclear spectrum within our multiconfigurational approach. We find that the hyperfine interaction is about one order of magnitude greater than that for Tb(III)Pc$_2$ SMMs. This stems from the strong Fermi contact interaction between the Tb nuclear spin and the electron spin density at the nucleus that originates from the occupation of the $(6s,5d)$ orbitals. We also uncover that the response of the Fermi contact term to electric field results in electrical tuning of the electronic-nuclear level separations. This hyperfine Stark effect may be useful for applications of molecular nuclear spins for quantum computing.
△ Less
Submitted 30 July, 2020;
originally announced July 2020.
-
Hyperfine and quadrupole interactions for Dy isotopes in DyPc$_2$ molecules
Authors:
Aleksander L. Wysocki,
Kyungwha Park
Abstract:
Nuclear spin levels play an important role in understanding magnetization dynamics and implementation and control of quantum bits in lanthanide-based single-molecule magnets. We investigate the hyperfine and nuclear quadrupole interactions for $^{161}$Dy and $^{163}$Dy nucleus in anionic DyPc$_2$ (Pc=phthalocyanine) single-molecule magnets, using multiconfigurational ab-initio methods (beyond dens…
▽ More
Nuclear spin levels play an important role in understanding magnetization dynamics and implementation and control of quantum bits in lanthanide-based single-molecule magnets. We investigate the hyperfine and nuclear quadrupole interactions for $^{161}$Dy and $^{163}$Dy nucleus in anionic DyPc$_2$ (Pc=phthalocyanine) single-molecule magnets, using multiconfigurational ab-initio methods (beyond density-functional theory) including spin-orbit interaction. The two isotopes of Dy are chosen because the others have zero nuclear spin. Both isotopes have the nuclear spin $I=5/2$, although the magnitude and sign of the nuclear magnetic moment differ from each other. The large energy gap between the electronic ground and first-excited Kramers doublets, allows us to map the microscopic hyperfine and quadrupole interaction Hamiltonian onto an effective Hamiltonian with an electronic pseudo-spin $S_{\rm eff}=1/2$ that corresponds to the ground Kramers doublet. Our ab-initio calculations show that the coupling between the nuclear spin and electronic orbital angular momentum contributes the most to the hyperfine interaction and that both the hyperfine and nuclear quadrupole interactions for $^{161}$Dy and $^{163}$Dy nucleus are much smaller than those for $^{159}$Tb nucleus in TbPc$_2$ single-molecule magnets. The calculated separations of the electronic-nuclear levels are comparable to experimental data reported for $^{163}$DyPc$_2$. We demonstrate that hyperfine interaction for Dy Kramers ion leads to tunnel splitting (or quantum tunneling of magnetization) at zero field. This effect does not occur for TbPc$_2$ single-molecule magnets. The magnetic field values of the avoided level crossings for $^{161}$DyPc$_2$ and $^{163}$DyPc$_2$ are found to be noticeably different, which can be observed from experiment.
△ Less
Submitted 12 February, 2020;
originally announced February 2020.
-
Nature of hyperfine interactions in TbPc$_2$ single-molecule magnets: Multireference ab-initio study
Authors:
Aleksander L. Wysocki,
Kyungwha Park
Abstract:
Lanthanide-based single-ion magnetic molecules can have large magnetic hyperfine interactions as well as large magnetic anisotropy. Recent experimental studies reported tunability of these properties by changes of chemical environments or by application of external stimuli for device applications. In order to provide insight onto the origin and mechanism of such tunability, here we investigate the…
▽ More
Lanthanide-based single-ion magnetic molecules can have large magnetic hyperfine interactions as well as large magnetic anisotropy. Recent experimental studies reported tunability of these properties by changes of chemical environments or by application of external stimuli for device applications. In order to provide insight onto the origin and mechanism of such tunability, here we investigate the magnetic hyperfine and nuclear quadrupole interactions for $^{159}$Tb nucleus in TbPc$_2$ (Pc=phthalocyanine) single-molecule magnets using multireference ab-initio methods including spin-orbit interaction. Since the electronic ground and first-excited (quasi)doublets are well separated in energy, the microscopic Hamiltonian can be mapped onto an effective Hamiltonian with an electronic pseudo-spin $S=1/2$. From the ab-initio-calculated parameters, we find that the magnetic hyperfine coupling is dominated by the interaction of the Tb nuclear spin with electronic orbital angular momentum. The asymmetric $4f$-like electronic charge distribution leads to a strong nuclear quadrupole interaction with significant non-axial terms for the molecule with low symmetry. The ab-initio calculated electronic-nuclear spectrum including the magnetic hyperfine and quadrupole interactions is in excellent agreement with experiment. We further find that the non-axial quadrupole interactions significantly influence the avoided level crossings in magnetization dynamics and that the molecular distortions affect mostly the Fermi contact terms as well as the non-axial quadrupole interactions.
△ Less
Submitted 25 October, 2019;
originally announced October 2019.
-
Computational models for active matter
Authors:
M Reza Shaebani,
Adam Wysocki,
Roland G Winkler,
Gerhard Gompper,
Heiko Rieger
Abstract:
A variety of computational models have been developed to describe active matter at different length and time scales. The diversity of the methods and the challenges in modeling active matter---ranging from molecular motors and cytoskeletal filaments over artificial and biological swimmers on microscopic to groups of animals on macroscopic scales---mainly originate from their out-of-equilibrium cha…
▽ More
A variety of computational models have been developed to describe active matter at different length and time scales. The diversity of the methods and the challenges in modeling active matter---ranging from molecular motors and cytoskeletal filaments over artificial and biological swimmers on microscopic to groups of animals on macroscopic scales---mainly originate from their out-of-equilibrium character, multiscale nature, nonlinearity, and multibody interactions. In the present review, various modeling approaches and numerical techniques are addressed, compared, and differentiated to illuminate the innovations and current challenges in understanding active matter. The complexity increases from minimal microscopic models of dry active matter toward microscopic models of active matter in fluids. Complementary, coarse-grained descriptions and continuum models are elucidated. Microscopic details are often relevant and strongly affect collective behaviors, which implies that the selection of a proper level of modeling is a delicate choice, with simple models emphasizing universal properties and detailed models capturing specific features. Finally, current approaches to further advance the existing models and techniques to cope with real-world applications, such as complex media and biological environments, are discussed.
△ Less
Submitted 13 October, 2019; v1 submitted 6 October, 2019;
originally announced October 2019.
-
Capillary action in scalar active matter
Authors:
Adam Wysocki,
Heiko Rieger
Abstract:
We study the capacity of active matter to rise in thin tubes against gravity and other related phenomena, like, wetting of vertical plates and spontaneous imbibition, where a wetting liquid is drawn into a porous medium. This capillary action or capillarity is well known in classical fluids and originates from attractive interactions between the liquid molecules and the container walls, and from t…
▽ More
We study the capacity of active matter to rise in thin tubes against gravity and other related phenomena, like, wetting of vertical plates and spontaneous imbibition, where a wetting liquid is drawn into a porous medium. This capillary action or capillarity is well known in classical fluids and originates from attractive interactions between the liquid molecules and the container walls, and from the attraction of the liquid molecules among each other. We observe capillarity in a minimal model for scalar active matter with purely repulsive interactions, where an effective attraction emerges due to slowdown during collisions between active particles and between active particles and walls. Simulations indicate that the capillary rise in thin tubes is approximately proportional to the active sedimentation length $λ$ and that the wetting height of a vertical plate grows superlinear with $λ$. In a disordered porous medium the imbibition height scales as $\langle h\rangle\proptoλφ_m$, where $φ_m$ is its packing fraction.
△ Less
Submitted 9 August, 2019;
originally announced August 2019.
-
Concentration tuned tetragonal strain in alloys: application to magnetic anisotropy of FeNi$_{1-x}$Co$_x$
Authors:
Aleksander L. Wysocki,
Manh Cuong Nguyen,
Cai-Zhuang Wang,
Kai-Ming Ho,
Andrey V. Postnikov,
Vladimir P. Antropov
Abstract:
We explore an opportunity to induce and control tetragonal distortion in materials. The idea involves formation of a binary alloy from parent compounds having body-centered and face-centered symmetries. The concept is illustrated in the case of FeNi$_{1-x}$Co$_x$ magnetic alloy formed by substitutional doping of the L1$_0$ FeNi magnet with Co. Using electronic structure calculations we demonstrate…
▽ More
We explore an opportunity to induce and control tetragonal distortion in materials. The idea involves formation of a binary alloy from parent compounds having body-centered and face-centered symmetries. The concept is illustrated in the case of FeNi$_{1-x}$Co$_x$ magnetic alloy formed by substitutional doping of the L1$_0$ FeNi magnet with Co. Using electronic structure calculations we demonstrate that the tetragonal strain in this system can be controlled by concentration and it reaches maximum for $x=0.5$. This finding is then applied to create a large magnetocrystalline anisotropy (MAE) in FeNi$_{1-x}$Co$_x$ system by considering an interplay of the tetragonal distortion with electronic concentration and chemical anisotropy. In particular, we identify a new ordered FeNi$_{0.5}$Co$_{0.5}$ system with MAE larger by a factor 4.5 from the L1$_0$ FeNi magnet.
△ Less
Submitted 26 June, 2019;
originally announced June 2019.
-
Multireference Ab Initio Studies of Magnetic Properties of Terbium-Based Single-Molecule Magnets
Authors:
Ryan Pederson,
Aleksander L. Wysocki,
Nicholas Mayhall,
Kyungwha Park
Abstract:
We investigate how different chemical environment influences magnetic properties of terbium(III) (Tb)-based single-molecule magnets (SMMs), using first-principles relativistic multireference methods. Recent experiments showed that Tb-based SMMs can have exceptionally large magnetic anisotropy and that they can be used for experimental realization of quantum information applications, with a judicio…
▽ More
We investigate how different chemical environment influences magnetic properties of terbium(III) (Tb)-based single-molecule magnets (SMMs), using first-principles relativistic multireference methods. Recent experiments showed that Tb-based SMMs can have exceptionally large magnetic anisotropy and that they can be used for experimental realization of quantum information applications, with a judicious choice of chemical environment. Here, we perform complete active space self-consistent field (CASSCF) calculations including relativistic spin-orbit interaction (SOI) for representative Tb-based SMMs such as TbPc$_2$ and TbPcNc in three charge states. We calculate low-energy electronic structure from which we compute the Tb crystal-field parameters and construct an effective pseudospin Hamiltonian. Our calculations show that ligand type and fine points of molecular geometry do not affect the zero-field splitting, while the latter varies weakly with oxidation number. On the other hand, higher-energy levels have a strong dependence on all these characteristics. For neutral TbPc$_2$ and TbPcNc molecules, the Tb magnetic moment and the ligand spin are parallel to each other and the coupling strength between them does not depend much on ligand type and details of atomic structure. However, ligand distortion and molecular symmetry play a crucial role in transverse crystal-field parameters which lead to tunnel splitting. The tunnel splitting induces quantum tunneling of magnetization by itself or by combining with other processes. Our results provide insight into mechanisms of magnetization relaxation in the representative Tb-based SMMs.
△ Less
Submitted 25 May, 2019;
originally announced May 2019.
-
Toward long-range entanglement between electrically driven single-molecule magnets
Authors:
Khadijeh Najafi,
Alexander Wysocki,
Kyungwha Park,
Sophia E. Economou,
Edwin Barnes
Abstract:
Over the past two decades, several molecules have been explored as possible building blocks of a quantum computer, a device that would provide exponential speedups for a number of problems, including the simulation of large, strongly correlated chemical systems. Achieving strong interactions and entanglement between molecular qubits remains an outstanding challenge. Here, we show that the TbPc…
▽ More
Over the past two decades, several molecules have been explored as possible building blocks of a quantum computer, a device that would provide exponential speedups for a number of problems, including the simulation of large, strongly correlated chemical systems. Achieving strong interactions and entanglement between molecular qubits remains an outstanding challenge. Here, we show that the TbPc$_2$ single-molecule magnet has the potential to overcome this obstacle due to its sensitivity to electric fields stemming from the hyperfine Stark effect. We show how this feature can be leveraged to achieve long-range entanglement between pairs of molecules using a superconducting resonator as a mediator. Our results suggest that the molecule-resonator interaction is near the edge of the strong-coupling regime and could potentially pass into it given a more detailed, quantitative understanding of the TbPc$_2$ molecule.
△ Less
Submitted 7 February, 2020; v1 submitted 22 May, 2019;
originally announced May 2019.
-
Single Crystal Permanent Magnet: Extraordinary Magnetic Behavior in the ta, Cu and Fe Substituted CeCo5 System
Authors:
Tej N. Lamichhane,
Michael T. Onyszczak,
Olena Palasyuk,
Saba Sharikadze,
Tae-Hoon Kim,
Matthew J. Kramer,
R. W. MacCallum,
Aleksander L. Wysocki,
Manh Cuong Nguyen,
Vladimir P. Antropov,
Tribhuvan Pandey,
David Parker,
Sergey L. Bud'ko,
Paul C. Canfield,
Andriy Palasyuk
Abstract:
To reduce material and processing costs of commercial permanent magnets and to attempt to fill the empty niche of energy products, 10 - 20 MGOe, between low-flux (ferrites, alnico) and high-flux (Nd2Fe14B- and SmCo5-type) magnets, we report synthesis, structure, magnetic properties and modeling of Ta, Cu and Fe substituted CeCo5. Using a self-flux technique, we grew single crystals of I - Ce15.1Ta…
▽ More
To reduce material and processing costs of commercial permanent magnets and to attempt to fill the empty niche of energy products, 10 - 20 MGOe, between low-flux (ferrites, alnico) and high-flux (Nd2Fe14B- and SmCo5-type) magnets, we report synthesis, structure, magnetic properties and modeling of Ta, Cu and Fe substituted CeCo5. Using a self-flux technique, we grew single crystals of I - Ce15.1Ta1.0Co74.4Cu9.5, II - Ce16.3Ta0.6Co68.9Cu14.2, III - Ce15.7Ta0.6Co67.8Cu15.9, IV - Ce16.3Ta0.3Co61.7Cu21.7 and V - Ce14.3Ta1.0Co62.0Fe12.3Cu10.4. X-ray diffraction analysis (XRD) showed that these materials retain a CaCu5 substructure and incorporate small amounts of Ta in the form of \dumb-bells", filling the 2e crystallographic sites within the 1D hexagonal channel with the 1a Ce site, whereas Co, Cu and Fe are statistically distributed among the 2c and 3g crystallographic sites. Scanning electron microscopy, energy dispersive X-ray spectroscopy (SEM-EDS) and scanning transmission electron microscopy (STEM) examinations provided strong evidence of the single-phase nature of the as-grown crystals, even though they readily exhibited significant magnetic coercivitie of ~1.6 - ~1.8 kOe caused by Co-enriched, nano-sized, structural defects and faults that can serve as pinning sites. Formation of the "composite crystal" during the heat treatment creates a 3D array of extended defects within a primarily single grain single crystal, which greatly improves its magnetic characteristics. Possible causes of the formation of the "composite crystal" may be associated with Ta atoms leaving matrix interstices at lower temperatures and/or matrix degradation induced by decreased miscibility at lower temperatures. Fe strongly improves both the Curie temperature and magnetization of the system resulting in (BH)max:~13 MGOe at room temperature.
△ Less
Submitted 16 May, 2018;
originally announced May 2018.
-
Spin-density fluctuations and the fluctuation-dissipation theorem in 3d ferromagnetic metals
Authors:
A. L. Wysocki,
V. N. Valmispild,
A. Kutepov,
S. Sharma,
J. K. Dewhurst,
E. K. U. Gross,
A. I. Lichtenstein,
V. P. Antropov
Abstract:
Spatial and time scales of spin density fluctuations (SDF) were analyzed in 3d ferromagnets using ab initio linear response calculations of complete wavevector and energy dependence of the dynamic spin susceptibility tensor. We demonstrate that SDF are spread continuously over the entire Brillouin zone and while majority of them reside within the 3d bandwidth, a significant amount comes from much…
▽ More
Spatial and time scales of spin density fluctuations (SDF) were analyzed in 3d ferromagnets using ab initio linear response calculations of complete wavevector and energy dependence of the dynamic spin susceptibility tensor. We demonstrate that SDF are spread continuously over the entire Brillouin zone and while majority of them reside within the 3d bandwidth, a significant amount comes from much higher energies. A validity of the adiabatic approximation in spin dynamics is discussed. The SDF spectrum is shown to have two main constituents: a minor low-energy spin wave contribution and a much larger high-energy component from more localized excitations. Using the fluctuation-dissipation theorem (FDT), the on-site spin correlator (SC) and the related effective fluctuating moment were properly evaluated and their universal dependence on the 3d band population is further discussed.
△ Less
Submitted 15 November, 2017; v1 submitted 1 June, 2017;
originally announced June 2017.
-
Strength and scales of itinerant spin fluctuations in 3d paramagnetic metals
Authors:
Aleksander L. Wysocki,
Andrey Kutepov,
Vladimir P. Antropov
Abstract:
The full spin density fluctuations (SDF) spectra in 3d paramagnetic metals are analyzed from first principles using the linear response technique. Using the calculated complete wavevector and energy dependence of the dynamic spin susceptibility, we obtain the most important, but elusive, characteristic of SDF in solids: on-site spin correlator (SC). We demonstrate that the SDF have a mixed charact…
▽ More
The full spin density fluctuations (SDF) spectra in 3d paramagnetic metals are analyzed from first principles using the linear response technique. Using the calculated complete wavevector and energy dependence of the dynamic spin susceptibility, we obtain the most important, but elusive, characteristic of SDF in solids: on-site spin correlator (SC). We demonstrate that the SDF have a mixed character consisting of interacting collective and single-particle excitations of similar strength spreading continuously over the entire Brillouin zone and a wide energy range up to femtosecond time scales. These excitations cannot be adiabatically separated and their intrinsically multiscale nature should be always taken into account for a proper description of metallic systems. Overall, in all studied systems, despite the lack of local moment, we found a very large SC resulting in an effective fluctuating moment of the order of several Bohr magnetons.
△ Less
Submitted 22 September, 2016;
originally announced September 2016.
-
Traveling fronts in active-passive particle mixtures
Authors:
Adam Wysocki,
Roland G. Winkler,
Gerhard Gompper
Abstract:
The emergent dynamics in phase-separated mixtures of isometric active and passive Brownian particles is studied numerically in two dimensions. A novel steady-state of well-defined traveling fronts is observed, where the interface between the dense and the dilute phase propagates and the bulk of both phases is (nearly) at rest. Two kind of interfaces, advancing and receding, are formed by spontaneo…
▽ More
The emergent dynamics in phase-separated mixtures of isometric active and passive Brownian particles is studied numerically in two dimensions. A novel steady-state of well-defined traveling fronts is observed, where the interface between the dense and the dilute phase propagates and the bulk of both phases is (nearly) at rest. Two kind of interfaces, advancing and receding, are formed by spontaneous symmetry breaking, induced by an instability of a planar interface due to the formation of localized vortices. The propagation arises due to flux imbalance at the interface, strongly resembling traveling fronts in reaction-diffusion systems. Above a threshold, the interface velocity decreases linearly with increasing fraction of active particles.
△ Less
Submitted 5 January, 2016;
originally announced January 2016.
-
Micromagnetic simulations with periodic boundary conditions: Hard-soft nanocomposites
Authors:
Aleksander L. Wysocki,
Vladimir P. Antropov
Abstract:
We developed a micromagnetic method for modeling magnetic systems with periodic boundary conditions along an arbitrary number of dimensions. The main feature is an adaptation of the Ewald summation technique for evaluation of long-range dipolar interactions. The method was applied to investigate the hysteresis process in hard-soft magnetic nanocomposites with various geometries. The dependence of…
▽ More
We developed a micromagnetic method for modeling magnetic systems with periodic boundary conditions along an arbitrary number of dimensions. The main feature is an adaptation of the Ewald summation technique for evaluation of long-range dipolar interactions. The method was applied to investigate the hysteresis process in hard-soft magnetic nanocomposites with various geometries. The dependence of the results on different micromagnetic parameters was studied. We found that for layered structures with an out-of-plane hard phase easy axis the hysteretic properties are very sensitive to the strength of the interlayer exchange coupling, as long as the spontaneous magnetization for the hard phase is significantly smaller than for the soft phase. The origin of this behavior was discussed. Additionally, we investigated the soft phase size optimizing the energy product of hard-soft nanocomposites.
△ Less
Submitted 28 October, 2015;
originally announced October 2015.
-
Magnetically induced phonon splitting in ACr$_2$O$_4$ spinels from first principles
Authors:
Aleksander L. Wysocki,
Turan Birol
Abstract:
We study the magnetically-induced phonon splitting in cubic ACr$_2$O$_4$ (A=Mg, Zn, Cd, Hg) spinels from first principles, and demonstrate that the sign of the splitting, which is experimentally observed to be opposite in CdCr$_2$O$_4$ compared to ZnCr$_2$O$_4$ and MgCr$_2$O$_4$, is determined solely by the particular magnetic ordering pattern observed in these compounds. We further show that this…
▽ More
We study the magnetically-induced phonon splitting in cubic ACr$_2$O$_4$ (A=Mg, Zn, Cd, Hg) spinels from first principles, and demonstrate that the sign of the splitting, which is experimentally observed to be opposite in CdCr$_2$O$_4$ compared to ZnCr$_2$O$_4$ and MgCr$_2$O$_4$, is determined solely by the particular magnetic ordering pattern observed in these compounds. We further show that this interaction between magnetism and phonon frequencies can be fully described by the previously proposed spin-phonon coupling model that includes only the nearest neighbor exchange. Finally, using this model with materials specific parameters calculated from first principles, we provide additional insights into the physics of spin-phonon coupling in this intriguing family of compounds.
△ Less
Submitted 4 August, 2015;
originally announced August 2015.
-
Virial pressure in systems of active Brownian particles
Authors:
Roland G. Winkler,
Adam Wysocki,
Gerhard Gompper
Abstract:
The pressure of suspensions of self-propelled objects is studied theoretically and by simulation of spherical active Brownian particles (ABP). We show that for certain geometries, the mechanical pressure as force/area of a confined systems can equally be expressed by bulk properties, which implies the existence of an nonequilibrium equation of state. Exploiting the virial theorem, we derive expres…
▽ More
The pressure of suspensions of self-propelled objects is studied theoretically and by simulation of spherical active Brownian particles (ABP). We show that for certain geometries, the mechanical pressure as force/area of a confined systems can equally be expressed by bulk properties, which implies the existence of an nonequilibrium equation of state. Exploiting the virial theorem, we derive expressions for the pressure of ABPs confined by solid walls or exposed to periodic boundary conditions. In both cases, the pressure comprises three contributions: the ideal-gas pressure due to white-noise random forces, an activity-induce pressure (swim pressure), which can be expressed in terms of a product of the bare and a mean effective propulsion velocity, and the contribution by interparticle forces. We find that the pressure of spherical ABPs in confined systems explicitly depends on the presence of the confining walls and the particle-wall interactions, which has no correspondence in systems with periodic boundary conditions. Our simulations of three-dimensional APBs in systems with periodic boundary conditions reveal a pressure-concentration dependence that becomes increasingly nonmonotonic with increasing activity. Above a critical activity and ABP concentration, a phase transition occurs, which is reflected in a rapid and steep change of the pressure. We present and discuss the pressure for various activities and analyse the contributions of the individual pressure components.
△ Less
Submitted 12 June, 2015;
originally announced June 2015.
-
Giant adsorption of microswimmers: duality of shape asymmetry and wall curvature
Authors:
Adam Wysocki,
Jens Elgeti,
Gerhard Gompper
Abstract:
The effect of shape asymmetry of microswimmers on their adsorption capacity at confining channel walls is studied by a simple dumbbell model. For a shape polarity of a forward-swimming cone, like the stroke-averaged shape of a sperm, extremely long wall retention times are found, caused by a non-vanishing component of the propulsion force pointing steadily into the wall, which grows exponentially…
▽ More
The effect of shape asymmetry of microswimmers on their adsorption capacity at confining channel walls is studied by a simple dumbbell model. For a shape polarity of a forward-swimming cone, like the stroke-averaged shape of a sperm, extremely long wall retention times are found, caused by a non-vanishing component of the propulsion force pointing steadily into the wall, which grows exponentially with the self-propulsion velocity and the shape asymmetry. A direct duality relation between shape asymmetry and wall curvature is proposed and verified. Our results are relevant for the design microswimmer with controlled wall-adhesion properties. In addition, we confirm that pressure in active systems is strongly sensitive to the details of the particle-wall interactions.
△ Less
Submitted 11 May, 2015; v1 submitted 10 March, 2015;
originally announced March 2015.
-
First-principles microscopic model of exchange-driven magnetoelectric response with application to Cr$_2$O$_3$
Authors:
Sai Mu,
A. L. Wysocki,
K. D. Belashchenko
Abstract:
The exchange-driven contribution to the magnetoelectric susceptibility $\hatα$ is formulated using a microscopic model Hamiltonian coupling the spin degrees of freedom to lattice displacements and electric field, which may be constructed from first-principles data. Electronic and ionic contributions are sorted out, and the latter is resolved into a sum of contributions from different normal modes.…
▽ More
The exchange-driven contribution to the magnetoelectric susceptibility $\hatα$ is formulated using a microscopic model Hamiltonian coupling the spin degrees of freedom to lattice displacements and electric field, which may be constructed from first-principles data. Electronic and ionic contributions are sorted out, and the latter is resolved into a sum of contributions from different normal modes. If intrasublattice spin correlations can be neglected, the longitudinal component $α_\parallel$ becomes proportional to the product of magnetic susceptibility and sublattice magnetization in accordance with Rado's phenomenological model. As an illustration, the method is applied to analyze the temperature dependence of the longitudinal magnetoelectric susceptibility of Cr$_2$O$_3$ using first-principles calculations and the quantum pair cluster approximation for magnetic thermodynamics. A substantial electronic contribution is found, which is opposite to the ionic part. The sensitivity of the results to the Hubbard $U$ parameter and the sources of error are studied. It is also found that non-Heisenberg interactions are too weak to account for the sign change of $α_\parallel$ in Cr$_2$O$_3$.
△ Less
Submitted 23 April, 2014;
originally announced April 2014.
-
Hydrodynamic simulations of self-phoretic microswimmers
Authors:
Mingcheng Yang,
Adam Wysocki,
Marisol Ripoll
Abstract:
A mesoscopic hydrodynamic model to simulate synthetic self-propelled Janus particles which is thermophoretically or diffusiophoretically driven is here developed. We first propose a model for a passive colloidal sphere which reproduces the correct rotational dynamics together with strong phoretic effect. This colloid solution model employs a multiparticle collision dynamics description of the solv…
▽ More
A mesoscopic hydrodynamic model to simulate synthetic self-propelled Janus particles which is thermophoretically or diffusiophoretically driven is here developed. We first propose a model for a passive colloidal sphere which reproduces the correct rotational dynamics together with strong phoretic effect. This colloid solution model employs a multiparticle collision dynamics description of the solvent, and combines potential interactions with the solvent, with stick boundary conditions. Asymmetric and specific colloidal surface is introduced to produce the properties of self-phoretic Janus particles. A comparative study of Janus and microdimer phoretic swimmers is performed in terms of their swimming velocities and induced flow behavior. Self-phoretic microdimers display long range hydrodynamic interactions and can be characterized as pullers or pushers. In contrast, Janus particles are characterized by short range hydrodynamic interactions and behave as neutral swimmers. Our model nicely mimics those recent experimental realization of the self-phoretic Janus particles.
△ Less
Submitted 3 March, 2014;
originally announced March 2014.
-
Glass-Transition Properties from Hard Spheres to Charged Point Particles
Authors:
Anoosheh Yazdi,
Alexei Ivlev,
Sergei Khrapak,
Hubertus Thomas,
Gregor E. Morfill,
Hartmut Löwen,
Adam Wysocki,
Matthias Sperl
Abstract:
The glass transition is investigated in three dimensions for single and double Yukawa potentials for the full range of control parameters. For vanishing screening parameter, the limit of the one-component plasma is obtained; for large screening parameters and high coupling strengths, the glass-transition properties crossover to the hard-sphere system. Between the two limits, the entire transition…
▽ More
The glass transition is investigated in three dimensions for single and double Yukawa potentials for the full range of control parameters. For vanishing screening parameter, the limit of the one-component plasma is obtained; for large screening parameters and high coupling strengths, the glass-transition properties crossover to the hard-sphere system. Between the two limits, the entire transition diagram can be described by analytical functions. Different from other potentials, the glass-transition and melting lines for Yukawa potentials are found to follow shifted but otherwise identical curves in control-parameter space.
△ Less
Submitted 28 January, 2014;
originally announced January 2014.
-
Cooperative Motion of Active Brownian Spheres in Three-Dimensional Dense Suspensions
Authors:
Adam Wysocki,
Roland G. Winkler,
Gerhard Gompper
Abstract:
The structural and dynamical properties of suspensions of self-propelled Brownian particles of spherical shape are investigated in three spatial dimensions. Our simulations reveal a phase separation into a dilute and a dense phase, above a certain density and strength of self-propulsion. The packing fraction of the dense phase approaches random close packing at high activity, yet the system remain…
▽ More
The structural and dynamical properties of suspensions of self-propelled Brownian particles of spherical shape are investigated in three spatial dimensions. Our simulations reveal a phase separation into a dilute and a dense phase, above a certain density and strength of self-propulsion. The packing fraction of the dense phase approaches random close packing at high activity, yet the system remains fluid. Although no alignment mechanism exists in this model, we find long-lived cooperative motion of the particles in the dense regime. This behavior is probably due to an interface-induced sorting process. Spatial displacement correlation functions are nearly scale-free for systems with densities close to or above the glass transition density of passive systems.
△ Less
Submitted 11 May, 2015; v1 submitted 29 August, 2013;
originally announced August 2013.
-
Bulk magnetoelectricity in the hexagonal manganites and ferrites
Authors:
Hena Das,
Aleksander L. Wysocki,
Yanan Geng,
Weida Wu,
Craig J. Fennie
Abstract:
Improper ferroelectricity (trimerization) in the hexagonal manganites RMnO$_3$ leads to a network of coupled structural and magnetic vortices that induce domain wall magnetoelectricity and magnetization neither of which, however, occurs in the bulk. Here we combined first-principles calculations, group-theoretic techniques, and microscopic spin models to show how the trimerization not only induces…
▽ More
Improper ferroelectricity (trimerization) in the hexagonal manganites RMnO$_3$ leads to a network of coupled structural and magnetic vortices that induce domain wall magnetoelectricity and magnetization neither of which, however, occurs in the bulk. Here we combined first-principles calculations, group-theoretic techniques, and microscopic spin models to show how the trimerization not only induces a polarization but also a bulk magnetization and bulk magnetoelectric (ME) effect. This results in the existence of a bulk linear ME vortex structure or a bulk ME coupling such that if direction of polarization reverses so does magnetization. To measure the predicted ME vortex, we suggest RMnO$_3$ under large magnetic field. We suggest a family of materials, the hexagonal RFeO$_3$ ferrites, also display the predicted phenomena in their ground state.
△ Less
Submitted 1 August, 2013; v1 submitted 5 February, 2013;
originally announced February 2013.
-
Spin injection from a half-metal at finite temperatures
Authors:
K. D. Belashchenko,
J. K. Glasbrenner,
A. L. Wysocki
Abstract:
Spin injection from a half-metallic electrode in the presence of thermal spin disorder is analyzed using a combination of random matrix theory, spin-diffusion theory, and explicit simulations for the tight-binding s-d model. It is shown that efficient spin injection from a half-metal is possible as long as the effective resistance of the normal metal does not exceed a characteristic value, which d…
▽ More
Spin injection from a half-metallic electrode in the presence of thermal spin disorder is analyzed using a combination of random matrix theory, spin-diffusion theory, and explicit simulations for the tight-binding s-d model. It is shown that efficient spin injection from a half-metal is possible as long as the effective resistance of the normal metal does not exceed a characteristic value, which does not depend on the resistance of the half-metallic electrode, but is rather controlled by spin-flip scattering at the interface. This condition can be formulated as α<(l/L)/T, where αis the relative deviation of the magnetization from saturation, l and L the mean-free path and the spin-diffusion length in the non-magnetic channel, and T the transparency of the tunnel barrier at the interface (if present). The general conclusions are confirmed by tight-binding s-d model calculations. A rough estimate suggests that efficient spin injection from true half-metallic ferromagnets into silicon or copper may be possible at room temperature across a transparent interface.
△ Less
Submitted 27 November, 2012; v1 submitted 26 October, 2012;
originally announced October 2012.
-
Microscopic origin of the structural phase transitions at the Cr2O3 (0001) surface
Authors:
A. L. Wysocki,
Siqi Shi,
K. D. Belashchenko
Abstract:
The surface of a Cr2O3 (0001) film epitaxially grown on Cr undergoes an unusual reentrant sequence of structural phase transitions. In order to understand the underlying microscopic mechanisms, the structural and magnetic properties of the Cr2O3 (0001) surface are here studied using first-principles electronic structure calculations. Two competing surface Cr sites are identified. The energetics of…
▽ More
The surface of a Cr2O3 (0001) film epitaxially grown on Cr undergoes an unusual reentrant sequence of structural phase transitions. In order to understand the underlying microscopic mechanisms, the structural and magnetic properties of the Cr2O3 (0001) surface are here studied using first-principles electronic structure calculations. Two competing surface Cr sites are identified. The energetics of the surface is described by a configurational Hamiltonian with parameters determined using total energy calculations for several surface supercells. Effects of epitaxial strain and magnetic ordering on configurational interaction are also included. The thermodynamics of the system is studied using Monte Carlo simulations. At zero strain the surface undergoes an ordering phase transition at 165 K. Tensile epitaxial strain together with antiferromagnetic ordering drive the system toward strong configurational frustration, suggesting the mechanism for the disordering phase transition at lower temperatures.
△ Less
Submitted 19 October, 2012;
originally announced October 2012.
-
The Magnetoelectric Effect in Transition Metal Oxides: Insights and the Rational Design of New Materials from First Principles
Authors:
Turan Birol,
Nicole A. Benedek,
Hena Das,
Aleksander L. Wysocki,
Andrew T. Mulder,
Brian M. Abbett,
Eva H. Smith,
Saurabh Ghosh,
Craig J. Fennie
Abstract:
The search for materials displaying a large magnetoelectric effect has occupied researchers for many decades. The rewards could include not only advanced electronics technologies, but also fundamental insights concerning the dielectric and magnetic properties of condensed matter. In this article, we focus on the magnetoelectric effect in transition metal oxides and review the manner in which first…
▽ More
The search for materials displaying a large magnetoelectric effect has occupied researchers for many decades. The rewards could include not only advanced electronics technologies, but also fundamental insights concerning the dielectric and magnetic properties of condensed matter. In this article, we focus on the magnetoelectric effect in transition metal oxides and review the manner in which first-principles calculations have helped guide the search for (and increasingly, predicted) new materials and shed light on the microscopic mechanisms responsible for magnetoelectric phenomena.
△ Less
Submitted 16 October, 2012; v1 submitted 20 July, 2012;
originally announced July 2012.
-
Consistent model of magnetism in ferropnictides
Authors:
A. L. Wysocki,
K. D. Belashchenko,
V. P. Antropov
Abstract:
The discovery of superconductivity in LaFeAsO introduced the ferropnictides as a major new class of superconducting compounds with critical temperatures second only to cuprates. The presence of magnetic iron makes ferropnictides radically different from cuprates. Antiferromagnetism of the parent compounds strongly suggests that superconductivity and magnetism are closely related. However, the char…
▽ More
The discovery of superconductivity in LaFeAsO introduced the ferropnictides as a major new class of superconducting compounds with critical temperatures second only to cuprates. The presence of magnetic iron makes ferropnictides radically different from cuprates. Antiferromagnetism of the parent compounds strongly suggests that superconductivity and magnetism are closely related. However, the character of magnetic interactions and spin fluctuations in ferropnictides, in spite of vigorous efforts, has until now resisted understanding within any conventional model of magnetism. Here we show that the most puzzling features can be naturally reconciled within a rather simple effective spin model with biquadratic interactions, which is consistent with electronic structure calculations. By going beyond the Heisenberg model, this description explains numerous experimentally observed properties, including the peculiarities of the spin wave spectrum, thin domain walls, crossover from first to second order phase transition under doping in some compounds, and offers new insight in the occurrence of the nematic phase above the antiferromagnetic phase transition.
△ Less
Submitted 17 November, 2010; v1 submitted 8 November, 2010;
originally announced November 2010.
-
First-principles analysis of spin-disorder resistivity of Fe and Ni
Authors:
A. L. Wysocki,
R. F. Sabirianov,
M. van Schilfgaarde,
K. D. Belashchenko
Abstract:
Spin-disorder resistivity of Fe and Ni and its temperature dependence are analyzed using noncollinear density functional calculations within the supercell method. Different models of thermal spin disorder are considered, including the mean-field approximation and the nearest-neighbor Heisenberg model. Spin-disorder resistivity is found to depend weakly on magnetic short-range order. If the local…
▽ More
Spin-disorder resistivity of Fe and Ni and its temperature dependence are analyzed using noncollinear density functional calculations within the supercell method. Different models of thermal spin disorder are considered, including the mean-field approximation and the nearest-neighbor Heisenberg model. Spin-disorder resistivity is found to depend weakly on magnetic short-range order. If the local moments are kept frozen at their zero-temperature values, very good agreement with experiment is obtained for Fe, but for Ni the resistivity at elevated temperatures is significantly overestimated. Agreement with experiment for Fe is improved if the local moments are iterated to self-consistency. The overestimation of the resistivity for paramagnetic Ni is attributed to the reduction of the local moments down to 0.35 Bohr magnetons. Overall, the results suggest that low-energy spin fluctuations in Fe and Ni are better viewed as classical rotations of local moments rather than quantized spin fluctuations that would require an (S+1)/S correction.
△ Less
Submitted 15 May, 2009;
originally announced May 2009.
-
Oscillatory driven colloidal binary mixtures: axial segregation versus laning
Authors:
Adam Wysocki,
Hartmut Löwen
Abstract:
Using Brownian dynamics computer simulations we show that binary mixtures of colloids driven in opposite directions by an oscillating external field exhibit axial segregation in sheets perpendicular to the drive direction. The segregation effect is stable only in a finite window of oscillation frequencies and driving strengths and is taken over by lane formation in direction of the driving field…
▽ More
Using Brownian dynamics computer simulations we show that binary mixtures of colloids driven in opposite directions by an oscillating external field exhibit axial segregation in sheets perpendicular to the drive direction. The segregation effect is stable only in a finite window of oscillation frequencies and driving strengths and is taken over by lane formation in direction of the driving field if the driving force is increased. In the crossover regime, bands tilted relative to the drive direction are observed. Possible experiments to verify the axial segregation are discussed.
△ Less
Submitted 23 December, 2008;
originally announced December 2008.
-
Dynamics of lane formation in driven binary complex plasmas
Authors:
K. R. Sütterlin,
A. Wysocki,
A. V. Ivlev,
C. Räth,
H. M. Thomas,
M. Rubin-Zuzic,
W. J. Goedheer,
V. E. Fortov,
A. M. Lipaev,
V. I. Molotkov,
O. F. Petrov,
G. E. Morfill,
H. Löwen
Abstract:
The dynamical onset of lane formation is studied in experiments with binary complex plasmas under microgravity conditions. Small microparticles are driven and penetrate into a cloud of big particles, revealing a strong tendency towards lane formation. The observed time-resolved lane formation process is in good agreement with computer simulations of a binary Yukawa model with Langevin dynamics.…
▽ More
The dynamical onset of lane formation is studied in experiments with binary complex plasmas under microgravity conditions. Small microparticles are driven and penetrate into a cloud of big particles, revealing a strong tendency towards lane formation. The observed time-resolved lane formation process is in good agreement with computer simulations of a binary Yukawa model with Langevin dynamics. The laning is quantified in terms of the anisotropic scaling index, leading to a universal order parameter for driven systems.
△ Less
Submitted 16 December, 2008;
originally announced December 2008.
-
Thermodynamics of itinerant magnets in a classical spin fluctuation model
Authors:
A. L. Wysocki,
J. K. Glasbrenner,
K. D. Belashchenko
Abstract:
Thermodynamics of itinerant magnets is studied using a classical model with one parameter characterizing the degree of itinerancy. Monte Carlo simulations for bcc and fcc lattices are compared with the mean-field approximation and with the Onsager cavity field approximation extended to itinerant systems. The qualitative features of thermodynamics are similar to the known results of the functiona…
▽ More
Thermodynamics of itinerant magnets is studied using a classical model with one parameter characterizing the degree of itinerancy. Monte Carlo simulations for bcc and fcc lattices are compared with the mean-field approximation and with the Onsager cavity field approximation extended to itinerant systems. The qualitative features of thermodynamics are similar to the known results of the functional integral method. It is found that magnetic short-range order is weak and almost independent on the degree of itinerancy, and the mean-field approximation describes the thermodynamics reasonably well. Ambiguity of the phase space measure for classical models is emphasized. The Onsager cavity field method is extended to itinerant systems, which involves the renormalization of both the Weiss field and the on-site exchange interaction. The predictions of this approximation are in excellent agreement with Monte Carlo results.
△ Less
Submitted 28 October, 2008;
originally announced October 2008.
-
Magnetism of chromia
Authors:
Siqi Shi,
A. L. Wysocki,
K. D. Belashchenko
Abstract:
The electronic structure and magnetism of chromia (corundum-type Cr2O3) are studied using full-potential first-principles calculations. The electronic correlations are included within the LSDA+U method. The energies of different magnetic configurations are very well fitted by the Heisenberg Hamiltonian with strong exchange interaction with two nearest neighbors and additional weak interaction up…
▽ More
The electronic structure and magnetism of chromia (corundum-type Cr2O3) are studied using full-potential first-principles calculations. The electronic correlations are included within the LSDA+U method. The energies of different magnetic configurations are very well fitted by the Heisenberg Hamiltonian with strong exchange interaction with two nearest neighbors and additional weak interaction up to the fifth neighbor shell. These energies are insensitive to the position of the oxygen states, indicating that magnetism in Cr2O3 is dominated by direct exchange. The Neel temperature is calculated using the pair-cluster approximation for localized quantum spins of magnitude 3/2. Very good agreement with experiment is found for all properties including the equilibrium volume, spectral density, local magnetic moment, band gap, and the Neel temperature for the values of U and J that are close to those obtained within the constrained occupation method. The band gap is of the Mott-Hubbard type.
△ Less
Submitted 21 October, 2008;
originally announced October 2008.
-
Direct observation of hydrodynamic instabilities in driven non-uniform colloidal dispersions
Authors:
Adam Wysocki,
C. Patrick Royall,
Roland G. Winkler,
Gerhard Gompper,
Hajime Tanaka,
Alfons van Blaaderen,
Hartmut Löwen
Abstract:
A Rayleigh-Taylor-like instability of a dense colloidal layer under gravity in a capillary of microfluidic dimensions is considered. We access all relevant lengthscales with particle-level microscopy and computer simulations which incorporate long-range hydrodynamic interactions between the particles. By tuning the gravitational driving force, we reveal a mechanism whose growth is connected to t…
▽ More
A Rayleigh-Taylor-like instability of a dense colloidal layer under gravity in a capillary of microfluidic dimensions is considered. We access all relevant lengthscales with particle-level microscopy and computer simulations which incorporate long-range hydrodynamic interactions between the particles. By tuning the gravitational driving force, we reveal a mechanism whose growth is connected to the fluctuations of specific wavelengths, non-linear pattern formation and subsequent diffusion-dominated relaxation. Our linear stability theory captures the initial regime and thus predicts mixing conditions, with important implications for fields ranging from biology to nanotechnology.
△ Less
Submitted 7 October, 2008;
originally announced October 2008.
-
Calculations of spin-disorder resistivity from first principles
Authors:
A. L. Wysocki,
K. D. Belashchenko,
J. P. Velev,
M. van Schilfgaarde
Abstract:
Spin-disorder resistivity of Fe and Ni is studied using the noncollinear density functional theory. The Landauer conductance is averaged over random disorder configurations and fitted to Ohm's law. The distribution function is approximated by the mean-field theory. The dependence of spin-disorder resistivity on magnetization in Fe is found to be in excellent agreement with the results for the is…
▽ More
Spin-disorder resistivity of Fe and Ni is studied using the noncollinear density functional theory. The Landauer conductance is averaged over random disorder configurations and fitted to Ohm's law. The distribution function is approximated by the mean-field theory. The dependence of spin-disorder resistivity on magnetization in Fe is found to be in excellent agreement with the results for the isotropic s-d model. In the fully disordered state, spin-disorder resistivity for Fe is close to experiment, while for fcc Ni it exceeds the experimental value by a factor of 2.3. This result indicates strong magnetic short-range order in Ni at the Curie temperature.
△ Less
Submitted 20 November, 2006; v1 submitted 14 November, 2006;
originally announced November 2006.
-
Stick boundary conditions and rotational velocity auto-correlation functions for colloidal particles in a coarse-grained representation of the solvent
Authors:
J. T. Padding,
H. Löwen,
A. Wysocki,
A. A. Louis
Abstract:
We show how to implement stick boundary conditions for a spherical colloid in a solvent that is coarse-grained by the method of stochastic rotation dynamics. This allows us to measure colloidal rotational velocity auto-correlation functions by direct computer simulation. We find quantitative agreement with Enskog theory for short times and with hydrodynamic mode-coupling theory for longer times.…
▽ More
We show how to implement stick boundary conditions for a spherical colloid in a solvent that is coarse-grained by the method of stochastic rotation dynamics. This allows us to measure colloidal rotational velocity auto-correlation functions by direct computer simulation. We find quantitative agreement with Enskog theory for short times and with hydrodynamic mode-coupling theory for longer times. For aqueous colloidal suspensions, the Enskog contribution to the rotational friction is larger than the hydrodynamic one when the colloidal radius drops below 35nm.
△ Less
Submitted 29 October, 2005; v1 submitted 29 June, 2005;
originally announced June 2005.
-
Molecular view of the Rayleigh-Taylor instability in compressible Brownian fluids
Authors:
A. Wysocki,
H. Löwen
Abstract:
The onset of the Rayleigh-Taylor instability is studied a compressible Brownian Yukawa fluid mixture on the ``molecular'' length and time scales of the individual particles. As a model, a two-dimensional phase-separated symmetric binary mixture of colloidal particles of type $A$ and $B$ with a fluid-fluid interface separating an $A$-rich phase from a $B$-rich phase is investigated by Brownian co…
▽ More
The onset of the Rayleigh-Taylor instability is studied a compressible Brownian Yukawa fluid mixture on the ``molecular'' length and time scales of the individual particles. As a model, a two-dimensional phase-separated symmetric binary mixture of colloidal particles of type $A$ and $B$ with a fluid-fluid interface separating an $A$-rich phase from a $B$-rich phase is investigated by Brownian computer simulations when brought into non-equilibrium via a constant external driving field which acts differently on the different particles and perpendicular to the interface. Two different scenarios are observed which occur either for high or for low interfacial free energies as compared to the driving force. In the first scenario for high interfacial tension, the critical wavelength $λ_c$ of the unstable interface modes is in good agreement with the classical Rayleigh-Taylor formula provided dynamically rescaled values for the interfacial tension are used. The wavelength $λ_{c}$ increases with time representing a self-healing effect of the interface due to a local density increase near the interface. The Rayleigh-Taylor formula is confirmed even if $λ_c$ is of the order of a molecular correlation length. In the second scenario for very large driving forces as compared to the interfacial line tensions, on the other hand, the particle penetrate easily the interface by the driving field and form microscopic lanes with a width different from the predictions of the classical Rayleigh-Taylor formula.
△ Less
Submitted 11 February, 2004;
originally announced February 2004.