Showing 1–50 of 55 results for author: Wood, A

Self-Assembly and Phase Behavior of Janus Rods: Competition Between Shape and Potential Anisotropy

Authors: Jared A. Wood, Laura Dal Compare, Lillian Pearse, Alicia Schuitemaker, Yawei Liu, Toby Hudson, Achille Giacometti, Asaph Widmer-Cooper

Abstract: We characterize the self-assembly and phase behavior of Janus rods over a broad range of temperatures and volume fractions, using Langevin dynamics simulations and free energy calculations. The Janus rods consist of a line of fused overlapping spheres that interact via a soft-core repulsive potential, with the addition of an attractive pseudo-square-well tail to a fraction of the spheres (the cove… ▽ More We characterize the self-assembly and phase behavior of Janus rods over a broad range of temperatures and volume fractions, using Langevin dynamics simulations and free energy calculations. The Janus rods consist of a line of fused overlapping spheres that interact via a soft-core repulsive potential, with the addition of an attractive pseudo-square-well tail to a fraction of the spheres (the coverage) ranging from 5\% to 100\% of sites. Competition between the stability of liquid crystal phases originating from shape anisotropy and assembly driven by directional interactions gives rise to a rich polymorphism that depends on the coverage. At low density near the Boyle temperature, we observe the formation of spherical and tubular micelles at low coverage, while at higher coverage randomly oriented monolayers form as the attractive parts of the rods overlap. At higher density, bilayer structures appear and merge to form smectic and crystalline lamellar phases. All of these structures gradually become unstable as the temperature is increased until eventually regular nematic and smectic phases appear, consistent with the hard rod limit. Our results indicate that the intermediate regime where shape-entropic effects compete with anisotropic attractions provided by site specificity is rich in structural possibilities, and should help guide the design of rod-like colloids for specific applications. △ Less

Submitted 23 August, 2024; originally announced August 2024.

Comments: 13 pages, 13 figures, SI not included

Muon spectroscopy investigation of anomalous dynamic magnetism in NiI$_2$

Authors: T. L. Breeze, B. M. Huddart, A. Hernández-Melían, N. P. Bentley, D. A. Mayoh, G. D. A. Wood, G. Balakrishnan, J. Wilkinson, F. L. Pratt, S. J. Clark, T. Lancaster

Abstract: We present the results of muon-spin relaxation ($μ^{+}$SR) measurements of the van der Waals magnet NiI$_2$, which probe magnetic phase transitions at $T_{\mathrm{N1}}=73$K and $T_{\mathrm{N2}}=60$K. Supporting density functional theory (DFT) calculations allow the determination of a single muon stopping site whose magnetic environment is consistent with the proposed ground-state magnetic structur… ▽ More We present the results of muon-spin relaxation ($μ^{+}$SR) measurements of the van der Waals magnet NiI$_2$, which probe magnetic phase transitions at $T_{\mathrm{N1}}=73$K and $T_{\mathrm{N2}}=60$K. Supporting density functional theory (DFT) calculations allow the determination of a single muon stopping site whose magnetic environment is consistent with the proposed ground-state magnetic structure. $μ^{+}$SR measurements of the dynamics reveal behavior consistent with spin-wave excitations below $T_{\mathrm{N2}}$, with an additional contribution from a nonmagnetic component of the material. In the region $T_{\mathrm{N1}}<T<T_{\mathrm{N2}}$ the system is magnetically ordered throughout the bulk, but the character of the dynamics changes qualitatively, resulting in an unusual region of temperature-independent fluctuations. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: 5 pages, 3 figures, submitted to PRB letters

Photo-induced charge state dynamics of the neutral and negatively charged silicon vacancy centers in room-temperature diamond

Authors: G. Garcia-Arellano, G. I. López-Morales, N. B. Manson, J. Flick, A. A. Wood, C. A. Meriles

Abstract: The silicon vacancy (SiV) center in diamond is drawing much attention due to its optical and spin properties, attractive for quantum information processing and sensing. Comparatively little is known, however, about the dynamics governing SiV charge state interconversion mainly due to challenges associated with generating, stabilizing, and characterizing all possible charge states, particularly at… ▽ More The silicon vacancy (SiV) center in diamond is drawing much attention due to its optical and spin properties, attractive for quantum information processing and sensing. Comparatively little is known, however, about the dynamics governing SiV charge state interconversion mainly due to challenges associated with generating, stabilizing, and characterizing all possible charge states, particularly at room temperature. Here, we use multi-color confocal microscopy and density functional theory to examine photo-induced SiV recombination - from neutral, to single-, to double-negatively charged - over a broad spectral window in chemical-vapor-deposition diamond under ambient conditions. For the SiV0 to SiV- transition, we find a linear growth of the photo-recombination rate with laser power at all observed wavelengths, a hallmark of single photon dynamics. Laser excitation of SiV-, on the other hand, yields only fractional recombination into SiV2-, a finding we interpret in terms of a photo-activated electron tunneling process from proximal nitrogen atoms. △ Less

Submitted 16 March, 2024; originally announced March 2024.

Journal ref: Advanced Science, 2024

Field-orientation-dependent magnetic phases in GdRu$_2$Si$_2$ probed with muon-spin spectroscopy

Authors: B. M. Huddart, A. Hernández-Melián, G. D. A. Wood, D. A. Mayoh, M. Gomilšek, Z. Guguchia, C. Wang, S. J. Blundell, G. Balakrishnan, T. Lancaster

Abstract: Centrosymmetric GdRu$_2$Si$_2$ exhibits a variety of multi-Q magnetic states as a function of temperature and applied magnetic field, including a square skyrmion-lattice phase. The material's behavior is strongly dependent on the direction of the applied field, with different phase diagrams resulting for fields applied parallel or perpendicular to the crystallographic $c$ axis. Here, we present th… ▽ More Centrosymmetric GdRu$_2$Si$_2$ exhibits a variety of multi-Q magnetic states as a function of temperature and applied magnetic field, including a square skyrmion-lattice phase. The material's behavior is strongly dependent on the direction of the applied field, with different phase diagrams resulting for fields applied parallel or perpendicular to the crystallographic $c$ axis. Here, we present the results of muon-spin relaxation ($μ^+$SR) measurements on single crystals of GdRu$_2$Si$_2$. Our analysis is based on the computation of muon stopping sites and consideration of zero-point motion effects, allowing direct comparison with the underlying spin textures in the material. Using transverse-field $μ^+$SR with fields applied along either the [001] or [100] crystallographic directions, we distinguish between the magnetic phases in this system via their distinct muon response, providing additional evidence for the skyrmion and meron-lattice phases, while also suggesting the existence of RKKY-driven muon hyperfine coupling. Zero-field $μ^+$SR provides clear evidence for a transition between two distinct magnetically-ordered phases at 39 K. △ Less

Submitted 5 August, 2024; v1 submitted 14 March, 2024; originally announced March 2024.

Comments: 16 pages, 8 figures

3D-mapping and manipulation of photocurrent in an optoelectronic diamond device

Authors: A. A. Wood, D. J. McCloskey, N. Dontschuk, A. Lozovoi, R. M. Goldblatt, T. Delord, D. A. Broadway, J. -P. Tetienne, B. C. Johnson, K. T. Mitchell, C. T. -K. Lew, C. A. Meriles, A. M. Martin

Abstract: Characterising charge transport in a material is central to the understanding of its electrical properties, and can usually only be inferred from bulk measurements of derived quantities such as current flow. Establishing connections between host material impurities and transport properties in emerging electronics materials, such as wide bandgap semiconductors, demands new diagnostic methods tailor… ▽ More Characterising charge transport in a material is central to the understanding of its electrical properties, and can usually only be inferred from bulk measurements of derived quantities such as current flow. Establishing connections between host material impurities and transport properties in emerging electronics materials, such as wide bandgap semiconductors, demands new diagnostic methods tailored to these unique systems, and the presence of optically-active defect centers in these materials offers a non-perturbative, in-situ characterisation system. Here, we combine charge-state sensitive optical microscopy and photoelectric detection of nitrogen-vacancy (NV) centres to directly image the flow of charge carriers inside a diamond optoelectronic device, in 3D and with temporal resolution. We optically control the charge state of background impurities inside the diamond on-demand, resulting in drastically different current flow such as filamentary channels nucleating from specific, defective regions of the device. We then optically engineered conducting channels that control carrier flow, key steps towards optically reconfigurable, wide bandgap designer optoelectronics. We anticipate our approach might be extended to probe other wide-bandgap semiconductors (SiC, GaN) relevant to present and emerging electronic technologies. △ Less

Submitted 10 February, 2024; originally announced February 2024.

Comments: 6 pages main text + 11 pages supplement, 3 figures, 8 supplementary figures. Comments welcome

Extreme Metastability of Diamond and its Transformation to BC8 Post-Diamond Phase of Carbon

Authors: Kien Nguyen-Cong, Jonathan T. Willman, Joseph M. Gonzalez, Ashley S. Williams, Anatoly B. Belonoshko, Stan G. Moore, Aidan P. Thompson, Mitchell A. Wood, Jon H. Eggert, Marius Millot, Luis A. Zepeda-Ruiz, Ivan I. Oleynik

Abstract: Diamond possesses exceptional physical properties due to its remarkably strong carbon-carbon bonding, leading to significant resilience to structural transformations at very high pressures and temperatures. Despite several experimental attempts, synthesis and recovery of the theoretically predicted post-diamond BC8 phase remains elusive. Through quantum accurate, multi-million atom molecular dynam… ▽ More Diamond possesses exceptional physical properties due to its remarkably strong carbon-carbon bonding, leading to significant resilience to structural transformations at very high pressures and temperatures. Despite several experimental attempts, synthesis and recovery of the theoretically predicted post-diamond BC8 phase remains elusive. Through quantum accurate, multi-million atom molecular dynamics (MD) simulations, we have uncovered the extreme metastability of diamond at very high pressures, significantly exceeding its range of thermodynamic stability. We predict the post-diamond BC8 phase to be experimentally accessible only within a narrow high pressure-temperature region of the carbon phase diagram. The diamond to BC8 transformation proceeds through pre-melting followed by BC8 nucleation and growth in the metastable carbon liquid. We propose a double-shock compression pathway to achieve BC8 synthesis, which is currently being explored in theory-inspired experiments at the National Ignition Facility. △ Less

Submitted 22 January, 2024; v1 submitted 16 January, 2024; originally announced January 2024.

Probing Iron in Earth's Core With Molecular-Spin Dynamics

Authors: Svetoslav Nikolov, Kushal Ramakrishna, Andrew Rohskopf, Mani Lokamani, Julien Tranchida, John Carpenter, Attila Cangi, Mitchell A. Wood

Abstract: Dynamic compression of iron to Earth-core conditions is one of the few ways to gather important elastic and transport properties needed to uncover key mechanisms surrounding the geodynamo effect. Herein a new machine-learned ab-initio derived molecular-spin dynamics (MSD) methodology with explicit treatment for longitudinal spin-fluctuations is utilized to probe the dynamic phase-diagram of iron.… ▽ More Dynamic compression of iron to Earth-core conditions is one of the few ways to gather important elastic and transport properties needed to uncover key mechanisms surrounding the geodynamo effect. Herein a new machine-learned ab-initio derived molecular-spin dynamics (MSD) methodology with explicit treatment for longitudinal spin-fluctuations is utilized to probe the dynamic phase-diagram of iron. This framework uniquely enables an accurate resolution of the phase-transition kinetics and Earth-core elastic properties, as highlighted by compressional wave velocity and adiabatic bulk moduli measurements. In addition, a unique coupling of MSD with time-dependent density functional theory enables gauging electronic transport properties, critically important for resolving geodynamo dynamics. △ Less

Submitted 15 November, 2023; originally announced November 2023.

Comments: 3 Figures in main document, 8 Figures in the supplemental information

Janus helices: From fully attractive to hard helices

Authors: Laura Dal Compare, Flavio Romano, Jared A. Wood, Asaph Widmer-Cooper, Achille Giacometti

Abstract: The phase diagram of hard helices differs from its hard rods counterpart by the presence of chiral "screw" phases stemming from the characteristic helical shape, in addition to the conventional liquid crystal phases also found for rod-like particles. Using extensive Monte Carlo and Molecular Dynamics simulations, we study the effect of the addition of a short-range attractive tail representing sol… ▽ More The phase diagram of hard helices differs from its hard rods counterpart by the presence of chiral "screw" phases stemming from the characteristic helical shape, in addition to the conventional liquid crystal phases also found for rod-like particles. Using extensive Monte Carlo and Molecular Dynamics simulations, we study the effect of the addition of a short-range attractive tail representing solvent-induced interactions to a fraction of the sites forming the hard helices, ranging from a single-site attraction to fully attractive helices for a specific helical shape. Different temperature regimes exist for different fractions of the attractive sites, as assessed in terms of the relative Boyle temperatures, that are found to be rather insensitive to the specific shape of the helical particle. The temperature range probed by the present study is well above the corresponding Boyle temperatures, with the phase behaviour still mainly entropically dominated and with the existence and location of the various liquid crystal phases only marginally affected. The pressure in the equation of state is found to decrease upon increasing the fraction of attractive beads and/or on lowering the temperature at fixed volume fraction, as expected on physical grounds. All screw phases are found to be stable within the considered range of temperatures with the smectic phase becoming more stable on lowering the temperature. By contrast, the location of the transition lines do not display a simple dependence on the fraction of attractive beads in the considered range of temperatures. △ Less

Submitted 9 October, 2023; originally announced October 2023.

Comments: 17 pages, 13 figures, JCP in press, SI not included

Robustness of steady state and stochastic cyclicity in generalized coalescence-fragmentation models

Authors: Brennen T. Fagan, Niall J. MacKay, A. Jamie Wood

Abstract: Processes of coalescence and fragmentation are used to understand the time-evolution of the mass distribution of various systems and may result in a steady state or in stable deterministic or stochastic cycles. Motivated by applications in insurgency warfare we investigate coalescence-fragmentation systems. We begin with a simple model of size-biased coalescence accompanied by shattering into mono… ▽ More Processes of coalescence and fragmentation are used to understand the time-evolution of the mass distribution of various systems and may result in a steady state or in stable deterministic or stochastic cycles. Motivated by applications in insurgency warfare we investigate coalescence-fragmentation systems. We begin with a simple model of size-biased coalescence accompanied by shattering into monomers. Depending on the parameters this model has an approximately power-law-distributed steady state or stochastic cycles of alternating gelation and shattering. We conduct stochastic simulations of this model and its generalizations to include different kernel types, accretion and erosion, and various distributions of non-shattering fragmentation. Our central aim is to explore the robustness of the steady state and gel-shatter cycles to these variations. We show that an approximate power-law steady state persists with the addition of accretion and erosion, and with partial rather than total shattering. However, broader distributions of fragment sizes typically vitiate both the power law steady state and gel-shatter cyclicity. This work clarifies features shown in coalescence/fragmentation model simulations and elucidates the relationship between the microscopic dynamics and observed phenomena in this widely applicable interdisciplinary model type. △ Less

Submitted 6 October, 2023; originally announced October 2023.

Comments: 29 pages, 13 Figures

Dynamic Formation of Preferentially Lattice Oriented, Self Trapped Hydrogen Clusters

Authors: M. A. Cusentino, E. L. Sikorski, M. J. McCarthy, A. P. Thompson, M. A. Wood

Abstract: A series of MD and DFT simulations were performed to investigate hydrogen self-clustering and retention in tungsten. Using a newly develop machine learned interatomic potential, spontaneous formation of hydrogen platelets was observed after implanting low-energy hydrogen into tungsten at high fluxes and temperatures. The platelets formed along low miller index orientations and neighboring tetrahed… ▽ More A series of MD and DFT simulations were performed to investigate hydrogen self-clustering and retention in tungsten. Using a newly develop machine learned interatomic potential, spontaneous formation of hydrogen platelets was observed after implanting low-energy hydrogen into tungsten at high fluxes and temperatures. The platelets formed along low miller index orientations and neighboring tetrahedral and octahedral sites and could grow to over 50 atoms in size. High temperatures above 600 K and high hydrogen concentrations were needed to observe significant platelet formation. A critical platelet size of six hydrogen atoms was needed for long term stability. Platelets smaller than this were found to be thermally unstable within a few nanoseconds. To verify these observations, characteristic platelets from the MD simulations were simulated using large-scale DFT. DFT corroborated the MD results in that large platelets were also found to be dynamically stable for five or more hydrogen atoms. The LDOS from the DFT simulated platelets indicated that hydrogen atoms, particularly at the periphery of the platelet, were found to be at least as stable as hydrogen atoms in bulk tungsten. In addition, electrons were found to be localized around hydrogen atoms in the platelet itself and that hydrogen atoms up to 4.2 Angstrom away within the platelet were found to be bonded suggesting that the hydrogen atoms are interacting across longer distances than previously suggested. These results reveal a self-clustering mechanisms for hydrogen within tungsten in the absence of radiation induced or microstructural defects that could be a precursor to blistering and potentially explain the experimentally observed high hydrogen retention particularly in the near surface region. △ Less

Submitted 19 July, 2023; originally announced July 2023.

Comments: 10 pages, 8 figures

Exploring Model Complexity in Machine Learned Potentials for Simulated Properties

Authors: Andrew Rohskopf, James Goff, Dionysios Sema, Kiarash Gordiz, Ngoc Cuong Nguyen, Asegun Henry, Aidan P. Thompson, Mitchell A. Wood

Abstract: Machine learning (ML) enables the development of interatomic potentials that promise the accuracy of first principles methods while retaining the low cost and parallel efficiency of empirical potentials. While ML potentials traditionally use atom-centered descriptors as inputs, different models such as linear regression and neural networks can map these descriptors to atomic energies and forces. T… ▽ More Machine learning (ML) enables the development of interatomic potentials that promise the accuracy of first principles methods while retaining the low cost and parallel efficiency of empirical potentials. While ML potentials traditionally use atom-centered descriptors as inputs, different models such as linear regression and neural networks can map these descriptors to atomic energies and forces. This begs the question: what is the improvement in accuracy due to model complexity irrespective of choice of descriptors? We curate three datasets to investigate this question in terms of ab initio energy and force errors: (1) solid and liquid silicon, (2) gallium nitride, and (3) the superionic conductor LGPS. We further investigate how these errors affect simulated properties with these models and verify if the improvement in fitting errors corresponds to measurable improvement in property prediction. Since linear and nonlinear regression models have different advantages and disadvantages, the results presented herein help researchers choose models for their particular application. By assessing different models, we observe correlations between fitting quantity (e.g. atomic force) error and simulated property error with respect to ab initio values. Such observations can be repeated by other researchers to determine the level of accuracy, and hence model complexity, needed for their particular systems of interest. △ Less

Submitted 4 June, 2023; originally announced June 2023.

The Double-$Q$ Ground State with Topological Charge Stripes in the Skyrmion Candidate $\text{GdRu}_{\text{2}}\text{Si}_{\text{2}}$

Authors: G. D. A. Wood, D. D. Khalyavin, D. A. Mayoh, J. Bouaziz, A. E. Hall, S. J. R. Holt, F. Orlandi, P. Manuel, S. Blügel, J. B. Staunton, O. A. Petrenko, M. R. Lees, G. Balakrishnan

Abstract: $\text{GdRu}_{\text{2}}\text{Si}_{\text{2}}$ is a centrosymmetric magnet in which a skyrmion lattice has recently been discovered. Here, we investigate the magnetic structure of the zero field ground state using neutron diffraction on single crystal and polycrystalline $^{\text{160}}\text{GdRu}_{\text{2}}\text{Si}_{\text{2}}$. In addition to observing the principal propagation vectors $\mathbf{q}_… ▽ More $\text{GdRu}_{\text{2}}\text{Si}_{\text{2}}$ is a centrosymmetric magnet in which a skyrmion lattice has recently been discovered. Here, we investigate the magnetic structure of the zero field ground state using neutron diffraction on single crystal and polycrystalline $^{\text{160}}\text{GdRu}_{\text{2}}\text{Si}_{\text{2}}$. In addition to observing the principal propagation vectors $\mathbf{q}_{1}$ and $\mathbf{q}_{2}$, we discover higher order magnetic satellites, notably $\mathbf{q}_{1} + 2\mathbf{q}_{2}$. The appearance of these satellites are explained within the framework of a new double-$Q$ constant-moment solution. Using powder diffraction we implement a quantitative refinement of this model. This structure, which contains vortexlike motifs, is shown to have a one-dimensional topological charge density. △ Less

Submitted 26 June, 2023; v1 submitted 28 April, 2023; originally announced April 2023.

Comments: Accepted PRB Letter (14.4.2023)

Transferable Interatomic Potentials for Aluminum from Ambient Conditions to Warm Dense Matter

Authors: Sandeep Kumar, Hossein Tahmasbi, Kushal Ramakrishna, Mani Lokamani, Svetoslav Nikolov, Julien Tranchida, Mitchell A. Wood, Attila Cangi

Abstract: We present a study on the transport and materials properties of aluminum spanning from ambient to warm dense matter conditions using a machine-learned interatomic potential (ML-IAP). Prior research has utilized ML-IAPs to simulate phenomena in warm dense matter, but these potentials have often been calibrated for a narrow range of temperature and pressures. In contrast, we train a single ML-IAP ov… ▽ More We present a study on the transport and materials properties of aluminum spanning from ambient to warm dense matter conditions using a machine-learned interatomic potential (ML-IAP). Prior research has utilized ML-IAPs to simulate phenomena in warm dense matter, but these potentials have often been calibrated for a narrow range of temperature and pressures. In contrast, we train a single ML-IAP over a wide range of temperatures, using density functional theory molecular dynamics (DFT-MD) data. Our approach overcomes computational limitations of DFT-MD simulations, enabling us to study transport and materials properties of matter at higher temperatures and longer time scales. We demonstrate the ML-IAP transferability across a wide range of temperatures using molecular-dynamics (MD) by examining the thermal conductivity, diffusion coefficient, viscosity, sound velocity, and ion-ion structure factor of aluminum up to about 60,000 K, where we find good agreement with previous theoretical data. △ Less

Submitted 19 April, 2023; originally announced April 2023.

Revealing contributions to conduction from transport within ordered and disordered regions in highly doped conjugated polymers through analysis of temperature-dependent Hall measurements

Authors: William A. Wood, Ian E. Jacobs, Leszek J. Spalek, Yuxuan Huang, Chen Chen, Xinglong Ren, Henning Sirringhaus

Abstract: Hall effect measurements in doped polymer semiconductors are widely reported, but are difficult to interpret due to screening of Hall voltages by carriers undergoing incoherent transport. Here, we propose a refined analysis for such Hall measurements, based on measuring the Hall coefficient as a function of temperature, and modelling carriers as existing in a regime of variable "deflectability" (i… ▽ More Hall effect measurements in doped polymer semiconductors are widely reported, but are difficult to interpret due to screening of Hall voltages by carriers undergoing incoherent transport. Here, we propose a refined analysis for such Hall measurements, based on measuring the Hall coefficient as a function of temperature, and modelling carriers as existing in a regime of variable "deflectability" (i.e. how strongly they "feel" the magnetic part of the Lorentz force). By linearly interpolating each carrier between the extremes of no deflection and full deflection, we demonstrate that it is possible to extract the (time-averaged) concentration of deflectable charge carriers, $\left<n_d\right>$, the average, temperature-dependent mobility of those carriers, $\left<μ_d\right>(T)$, as well as the ratio of conductivity that comes from such deflectable transport, $d(T)$. Our method was enabled by the construction of an improved AC Hall measurement system, as well as an improved data extraction method. We measured Hall bar devices of ion-exchange doped films of PBTTT-C$_{14}$ from 10--300 K. Our analysis provides evidence for the proportion of conductivity arising from deflectable transport, $d(T)$, increasing with doping level, ranging between 15.4% and 16.4% at room temperature. When compared to total charge-carrier-density estimates from independent methods, the values of $\left<n_d\right>$ extracted suggest that carriers spend $\sim$37% of their time of flight being deflectable in the most highly doped of the devices measured here. The extracted values of $d(T)$ being less than half this value thus suggest that the limiting factor for conductivity in such highly doped devices is carrier mobility, rather than concentration. △ Less

Submitted 24 March, 2023; originally announced March 2023.

Comments: Accepted for publication in APS Physical Review Materials

Atomic Representations of Local and Global Chemistry in Complex Alloys

Authors: Megan J. McCarthy, Jacob Startt, Rémi Dingreville, Aidan P. Thompson, Mitchell A. Wood

Abstract: The exceptional properties observed in complex concentrated alloys (CCAs) arise from the interplay between crystalline order and chemical disorder at the atomic scale, complicating a unique determination of properties. In contrast to conventional alloys, CCA properties emerge as distributions due to varying local chemical environments and the specific scale of measurement. Currently there are few… ▽ More The exceptional properties observed in complex concentrated alloys (CCAs) arise from the interplay between crystalline order and chemical disorder at the atomic scale, complicating a unique determination of properties. In contrast to conventional alloys, CCA properties emerge as distributions due to varying local chemical environments and the specific scale of measurement. Currently there are few ways to quantitatively define, track, and compare local alloy compositions (versus a global label, i.e. equiatomic) contained in a CCA. Molecular dynamics is used here to build descriptive metrics that connect a global alloy composition to the diverse local alloy compositions that define it. A machine-learned interatomic potential for MoNbTaTi is developed and we use these metrics to investigate how property distributions change with excursions in global-local composition space. Short-range order is examined through the lens of local chemistry for the equiatomic composition, demonstrating stark changes in vacancy formation energy with local chemistry evolution. △ Less

Submitted 18 April, 2023; v1 submitted 7 March, 2023; originally announced March 2023.

Comments: Version 3: updated funding statement, fixed typos, otherwise unchanged. 15 pages, 6 main figures, 1 supplemental figure

Machine Learned Interatomic Potential for Dispersion Strengthened Plasma Facing Components

Authors: E. L. Sikorski, M. A. Cusentino, M. J. McCarthy, J. Tranchida, M. A. Wood, A. P. Thompson

Abstract: Tungsten (W) is a material of choice for the divertor material due to its high melting temperature, thermal conductivity, and sputtering threshold. However, W has a very high brittle-to-ductile transition temperature and at fusion reactor temperatures ($\geq$1000K) may undergo recrystallization and grain growth. Dispersion-strengthening W with zirconium carbide (ZrC) can improve ductility and limi… ▽ More Tungsten (W) is a material of choice for the divertor material due to its high melting temperature, thermal conductivity, and sputtering threshold. However, W has a very high brittle-to-ductile transition temperature and at fusion reactor temperatures ($\geq$1000K) may undergo recrystallization and grain growth. Dispersion-strengthening W with zirconium carbide (ZrC) can improve ductility and limit grain growth, but much of the effects of the dispersoids on microstructural evolution and thermomechanical properties at high temperature are still unknown. We present a machine learned Spectral Neighbor Analysis Potential (SNAP) for W-ZrC that can now be used to study these materials. In order to construct a potential suitable for large-scale atomistic simulations at fusion reactor temperatures, it is necessary to train on ab initio data generated for a diverse set of structures, chemical environments, and temperatures. Further accuracy and stability tests of the potential were achieved using objective functions for both material properties and high temperature stability. Validation of lattice parameters, surface energies, bulk moduli, and thermal expansion is confirmed on the optimized potential. Tensile tests of W/ZrC bicrystals show that while the W(110)-ZrC(111) C-terminated bicrystal has the highest ultimate tensile strength (UTS) at room temperature, observed strength decreases with increasing temperature. At 2500K, the terminating C layer diffuses into the W, resulting in a weaker W-Zr interface. Meanwhile, the W(110)-ZrC(111) Zr-terminated bicrystal has the highest UTS at 2500K. △ Less

Submitted 2 December, 2022; originally announced December 2022.

The Thermodynamic Origins of Chiral Twist in Monolayer Assemblies of Hard Rod-like Colloids

Authors: Yawei Liu, Jared A. Wood, Achille Giacometti, Asaph Widmer-Cooper

Abstract: The propagation of chirality across scales is a common but poorly understood phenomenon in soft matter. In this work, we use computer simulations to study chiral monolayer assemblies formed by hard rod-like colloidal particles in the presence of non-adsorbing polymer and characterize the thermodynamic driving forces responsible for the twisting. Simulations show that straight (achiral) rods assemb… ▽ More The propagation of chirality across scales is a common but poorly understood phenomenon in soft matter. In this work, we use computer simulations to study chiral monolayer assemblies formed by hard rod-like colloidal particles in the presence of non-adsorbing polymer and characterize the thermodynamic driving forces responsible for the twisting. Simulations show that straight (achiral) rods assemble into monolayers with a spontaneous twist that is either left- or right-handed, while helical (chiral) rods lead to assemblies with preferential chiral features that depend on their handedness and curliness. The onset of chirality in these monolayers can be traced back to small clusters formed at the initial stage of the self-assembly. In these microscopic monolayers, entropy drives twisting in ways that differ from the assumptions on which existing continuum theory is built. Depending on the geometry of the constituent rods, the preferred chiral twist can be driven by entropy gain of the polymers, or of the rods, or both. In addition, the variation of the polymer entropy with twist depends on changes in both the surface area and the volume of the monolayer. Rod fluctuations perpendicular to the monolayer also play an important role in stabilising the twisting. △ Less

Submitted 30 November, 2022; originally announced November 2022.

Comments: 14 pages, 12 figures

Journal ref: Nanoscale 2022, 14, 16837-16844

Room-temperature photo-chromism of silicon vacancy centers in CVD diamond

Authors: Alexander Wood, Artur Lozovoi, Zi-Huai Zhang, Sachin Sharma, Gabriel I. López-Morales, Harishankar Jayakumar, Nathalie P. de Leon, Carlos A. Meriles

Abstract: The silicon-vacancy (SiV) center in diamond is typically found in three stable charge states, SiV0, SiV- and SiV2-, but studying the processes leading to their formation is challenging, especially at room temperature, due to their starkly different photo-luminescence rates. Here, we use confocal fluorescence microscopy to activate and probe charge interconversion between all three charge states un… ▽ More The silicon-vacancy (SiV) center in diamond is typically found in three stable charge states, SiV0, SiV- and SiV2-, but studying the processes leading to their formation is challenging, especially at room temperature, due to their starkly different photo-luminescence rates. Here, we use confocal fluorescence microscopy to activate and probe charge interconversion between all three charge states under ambient conditions. In particular, we witness the formation of SiV0 via the two-step capture of diffusing, photo-generated holes, a process we expose both through direct SiV0 fluorescence measurements at low temperatures and confocal microscopy observations in the presence of externally applied electric fields. Further, we show that continuous red illumination induces the converse process, first transforming SiV0 into SiV-, then into SiV2-. Our results shed light on the charge dynamics of SiV and promise opportunities for nanoscale sensing and quantum information processing. △ Less

Submitted 23 November, 2022; originally announced November 2022.

Permutation-adapted complete and independent basis for atomic cluster expansion descriptors

Authors: James M. Goff, Charles Sievers, Mitchell A. Wood, Aidan P. Thompson

Abstract: Atomic cluster expansion (ACE) methods provide a systematic way to describe particle local environments of arbitrary body order. For practical applications it is often required that the basis of cluster functions be symmetrized with respect to rotations and permutations. Existing methodologies yield sets of symmetrized functions that are over-complete. These methodologies thus require an additiona… ▽ More Atomic cluster expansion (ACE) methods provide a systematic way to describe particle local environments of arbitrary body order. For practical applications it is often required that the basis of cluster functions be symmetrized with respect to rotations and permutations. Existing methodologies yield sets of symmetrized functions that are over-complete. These methodologies thus require an additional numerical procedure, such as singular value decomposition (SVD), to eliminate redundant functions. In this work, it is shown that analytical linear relationships for subsets of cluster functions may be derived using recursion and permutation properties of generalized Wigner symbols. From these relationships, subsets (blocks) of cluster functions can be selected such that, within each block, functions are guaranteed to be linearly independent. It is conjectured that this block-wise independent set of permutation-adapted rotation and permutation invariant (PA-RPI) functions forms a complete, independent basis for ACE. Along with the first analytical proofs of block-wise linear dependence of ACE cluster functions and other theoretical arguments, numerical results are offered to demonstrate this. The utility of the method is demonstrated in the development of an ACE interatomic potential for tantalum. Using the new basis functions in combination with Bayesian compressive sensing sparse regression, some high degree descriptors are observed to persist and help achieve high-accuracy models. △ Less

Submitted 26 February, 2024; v1 submitted 2 August, 2022; originally announced August 2022.

Machine Learning Interatomic Potential for Simulations of Carbon at Extreme Conditions

Authors: Jonathan T. Willman, Kien Nguyen-Cong, Ashley S. Williams, Anatoly B. Belonoshko, Stan G. Moore, Aidan P. Thompson, Mitchell A. Wood, Ivan I. Oleynik

Abstract: A Spectral Neighbor Analysis (SNAP) machine learning interatomic potential (MLIP) has been developed for simulations of carbon at extreme pressures (up to 5 TPa) and temperatures (up to 20,000 K). This was achieved using a large database of experimentally relevant quantum molecular dynamics (QMD) data, training the SNAP potential using a robust machine learning methodology, and performing extensiv… ▽ More A Spectral Neighbor Analysis (SNAP) machine learning interatomic potential (MLIP) has been developed for simulations of carbon at extreme pressures (up to 5 TPa) and temperatures (up to 20,000 K). This was achieved using a large database of experimentally relevant quantum molecular dynamics (QMD) data, training the SNAP potential using a robust machine learning methodology, and performing extensive validation against QMD and experimental data. The resultant carbon MLIP demonstrates unprecedented accuracy and transferability in predicting the carbon phase diagram, melting curves of crystalline phases, and the shock Hugoniot, all within 3% of QMD. By achieving quantum accuracy and efficient implementation on leadership class high performance computing systems, SNAP advances frontiers of classical MD simulations by enabling atomic-scale insights at experimental time and length scales. △ Less

Submitted 2 May, 2022; originally announced May 2022.

Rheological design of thickened alcohol-based hand rubs

Authors: Andreia F Silva, Tiffany A Wood, Daniel J M Hodgson, John R Royer, Job H J Thijssen, Alex Lips, Wilson C K Poon

Abstract: The handleability and sensory perception of hand sanitisers by consumers affect the hygiene outcome. Spillage may result in under-dosing and poor sensory properties can lead to under-utilisation. We first propose four principles (low run off, spreadability, smoothness and non-stickiness) for designing the rheology of thickened alcohol-based hand rubs with acceptable handleability and hand feel. We… ▽ More The handleability and sensory perception of hand sanitisers by consumers affect the hygiene outcome. Spillage may result in under-dosing and poor sensory properties can lead to under-utilisation. We first propose four principles (low run off, spreadability, smoothness and non-stickiness) for designing the rheology of thickened alcohol-based hand rubs with acceptable handleability and hand feel. We then evaluate a commercial hand gel and a variety of simplified formulations thickened with microgels (Carbopol 974P, Carbopol Ultrez 20 and Sepimax Zen), or linear polymers (Jaguar HP 120 COS), and evaluate them against these design criteria. All four additives provide acceptable spreadability by shear thinning to $η\approx 10^{-1} Pa \, s$ at $\dotγ\sim 10^3 s^{-1}$. Either the finite yield stress conferred by the microgels ($σ_y \gtrsim 10 Pa$) or the increase in low-shear viscosity provided by the linear polymer ($η\gtrsim 1 Pa \, s$ at $\dotγ\lesssim 0.1 s^{-1}$) give rise to acceptably low run-off. However, the formulation using the linear polymer shows a filament breakage time of $τ_{\rm b} \approx 1 s$ in capillary rheology, which may result in stickiness and therefore a less than optimal hand feel. △ Less

Submitted 5 May, 2022; v1 submitted 24 January, 2022; originally announced January 2022.

Comments: 10 pages, 8 figures

Effects of Fe deficiency and Co substitution in polycrystalline and single crystals of Fe$_{3}$GeTe$_{2}$

Authors: D. A. Mayoh, G. D. A. Wood, S. J. R. Holt, G. Beckett, E. J. L. Dekker, M. R. Lees, G. Balakrishnan

Abstract: Fe$_{3}$GeTe$_{2}$ is a two-dimensional van der Waals material with a ferromagnetic ground state and a maximum transition temperature $T_{\mathrm{c}}\sim225$ K. However, when Fe$_{3}$GeTe$_{2}$ is synthesized lower values of $T_{\mathrm{c}}$ are often reported. This is attributed to a deficiency in the Fe at the 2c site in the crystal structure. Here we investigate the effect of Fe deficiency and… ▽ More Fe$_{3}$GeTe$_{2}$ is a two-dimensional van der Waals material with a ferromagnetic ground state and a maximum transition temperature $T_{\mathrm{c}}\sim225$ K. However, when Fe$_{3}$GeTe$_{2}$ is synthesized lower values of $T_{\mathrm{c}}$ are often reported. This is attributed to a deficiency in the Fe at the 2c site in the crystal structure. Here we investigate the effect of Fe deficiency and the substitution of Co for Fe on the magnetic properties of this system. We have synthesized both polycrystalline material and single crystals by chemical vapor transport and the flux method, with the largest crystals obtained using the flux method. Cobalt substitution at the Fe site is found to significantly reduce the magnetic transition temperature. Crystals of Fe$_{3}$GeTe$_{2}$ grown by chemical vapor transport with $\sim 8\%$ excess Fe in the starting materials display an optimum Fe content and magnetic transition temperature. △ Less

Submitted 7 February, 2022; v1 submitted 17 November, 2021; originally announced November 2021.

Comments: 22 pages, 8 figures, 3 tables. Published online in ACS Crystal Growth and Design

Journal ref: Cryst. Growth Des. 21, 6786 (2021)

Quantum-accurate magneto-elastic predictions with classical spin-lattice dynamics

Authors: Svetoslav Nikolov, Mitchell A. Wood, Attila Cangi, Jean-Bernard Maillet, Mihai-Cosmin Marinica, Aidan P. Thompson, Michael P. Desjarlais, Julien Tranchida

Abstract: A data-driven framework is presented for building magneto-elastic machine-learning interatomic potentials (ML-IAPs) for large-scale spin-lattice dynamics simulations. The magneto-elastic ML-IAPs are constructed by coupling a collective atomic spin model with an ML-IAP. Together they represent a potential energy surface from which the mechanical forces on the atoms and the precession dynamics of th… ▽ More A data-driven framework is presented for building magneto-elastic machine-learning interatomic potentials (ML-IAPs) for large-scale spin-lattice dynamics simulations. The magneto-elastic ML-IAPs are constructed by coupling a collective atomic spin model with an ML-IAP. Together they represent a potential energy surface from which the mechanical forces on the atoms and the precession dynamics of the atomic spins are computed. Both the atomic spin model and the ML-IAP are parametrized on data from first-principles calculations. We demonstrate the efficacy of our data-driven framework across magneto-structural phase transitions by generating a magneto-elastic ML-IAP for α-iron. The combined potential energy surface yields excellent agreement with first-principles magneto-elastic calculations and quantitative predictions of diverse materials properties including bulk modulus, magnetization, and specific heat across the ferromagnetic-paramagnetic phase transition. △ Less

Submitted 2 February, 2021; v1 submitted 18 January, 2021; originally announced January 2021.

Comments: 15 pages, 6 figures, 2 tables

Journal ref: npj Computational Materials 7, 153 (2021)

Dynamic Response of Glassy Dispersions in a Nematic Liquid Crystal

Authors: Navneeta Katyan, Tiffany Annette Wood

Abstract: We present the first systematic measurements of the dynamic response of a glassy soft-solid formed on dispersing colloids in a thermotropic nematic liquid crystal. Through applying a simple theoretical model, we reveal that the nematic elasticity dominates at low frequencies (longest timescales) and prevents flow. The Ericksen number $E_R$, determines a critical frequency $ω_c$, above which the vi… ▽ More We present the first systematic measurements of the dynamic response of a glassy soft-solid formed on dispersing colloids in a thermotropic nematic liquid crystal. Through applying a simple theoretical model, we reveal that the nematic elasticity dominates at low frequencies (longest timescales) and prevents flow. The Ericksen number $E_R$, determines a critical frequency $ω_c$, above which the viscous component $G'' \propto ω^{1/2}$. We find that the value of tan$δ$ is independent of the colloidal volume fraction ($φ$), and is determined by the ratio of lengthscales associated with the defect-line density and the size of colloid-free nematic regions trapped within the composite. △ Less

Submitted 23 December, 2019; originally announced December 2019.

Comments: 4 pages, 3 figures

Combinatorial mappings of exclusion processes

Authors: Anthony J. Wood, Richard A. Blythe, Martin R. Evans

Abstract: We review various combinatorial interpretations and mappings of stationary-state probabilities of the totally asymmetric, partially asymmetric and symmetric simple exclusion processes (TASEP, PASEP, SSEP respectively). In these steady states, the statistical weight of a configuration is determined from a matrix product, which can be written explicitly in terms of generalised ladder operators. This… ▽ More We review various combinatorial interpretations and mappings of stationary-state probabilities of the totally asymmetric, partially asymmetric and symmetric simple exclusion processes (TASEP, PASEP, SSEP respectively). In these steady states, the statistical weight of a configuration is determined from a matrix product, which can be written explicitly in terms of generalised ladder operators. This lends a natural association to the enumeration of random walks with certain properties. Specifically, there is a one-to-many mapping of steady-state configurations to a larger state space of discrete paths, which themselves map to an even larger state space of number permutations. It is often the case that the configuration weights in the extended space are of a relatively simple form (e.g., a Boltzmann-like distribution). Meanwhile, various physical properties of the nonequilibrium steady state - such as the entropy - can be interpreted in terms of how this larger state space has been partitioned. These mappings sometimes allow physical results to be derived very simply, and conversely the physical approach allows some new combinatorial problems to be solved. This work brings together results and observations scattered in the combinatorics and statistical physics literature, and also presents new results. The review is pitched at statistical physicists who, though not professional combinatorialists, are competent and enthusiastic amateurs. △ Less

Submitted 10 February, 2020; v1 submitted 2 August, 2019; originally announced August 2019.

Comments: 56 pages, 21 figures

Solvable model of a many-filament Brownian ratchet

Authors: Anthony J. Wood, Richard A. Blythe, Martin R. Evans

Abstract: We construct and exactly solve a model of an extended Brownian ratchet. The model comprises an arbitrary number of heterogeneous, growing and shrinking filaments which together move a rigid membrane by a ratchet mechanism. The model draws parallels with the dynamics of actin filament networks at the leading edge of the cell. In the model, the filaments grow and contract stochastically. The model a… ▽ More We construct and exactly solve a model of an extended Brownian ratchet. The model comprises an arbitrary number of heterogeneous, growing and shrinking filaments which together move a rigid membrane by a ratchet mechanism. The model draws parallels with the dynamics of actin filament networks at the leading edge of the cell. In the model, the filaments grow and contract stochastically. The model also includes forces which derive from a potential dependent on the separation between the filaments and the membrane. These forces serve to attract the filaments to the membrane or generate a surface tension that prevents the filaments from dispersing. We derive an N -dimensional diffusion equation for the N filament-membrane separations, which allows the steady-state probability distribution function to be calculated exactly under certain conditions. These conditions are fulfilled by the physically relevant cases of linear and quadratic interaction potentials. The exact solution of the diffusion equation furnishes expressions for the average velocity of the membrane and critical system parameters for which the system stalls and has zero net velocity. In the case of a restoring force, the membrane velocity grows as the square root of the force constant, whereas it decreases once a surface tension is introduced. △ Less

Submitted 2 October, 2019; v1 submitted 1 July, 2019; originally announced July 2019.

Journal ref: Phys. Rev. E 100, 042122 (2019)

A Performance and Cost Assessment of Machine Learning Interatomic Potentials

Authors: Yunxing Zuo, Chi Chen, Xiangguo Li, Zhi Deng, Yiming Chen, Jörg Behler, Gábor Csányi, Alexander V. Shapeev, Aidan P. Thompson, Mitchell A. Wood, Shyue Ping Ong

Abstract: Machine learning of the quantitative relationship between local environment descriptors and the potential energy surface of a system of atoms has emerged as a new frontier in the development of interatomic potentials (IAPs). Here, we present a comprehensive evaluation of ML-IAPs based on four local environment descriptors --- Behler-Parrinello symmetry functions, smooth overlap of atomic positions… ▽ More Machine learning of the quantitative relationship between local environment descriptors and the potential energy surface of a system of atoms has emerged as a new frontier in the development of interatomic potentials (IAPs). Here, we present a comprehensive evaluation of ML-IAPs based on four local environment descriptors --- Behler-Parrinello symmetry functions, smooth overlap of atomic positions (SOAP), the Spectral Neighbor Analysis Potential (SNAP) bispectrum components, and moment tensors --- using a diverse data set generated using high-throughput density functional theory (DFT) calculations. The data set comprising bcc (Li, Mo) and fcc (Cu, Ni) metals and diamond group IV semiconductors (Si, Ge) is chosen to span a range of crystal structures and bonding. All descriptors studied show excellent performance in predicting energies and forces far surpassing that of classical IAPs, as well as predicting properties such as elastic constants and phonon dispersion curves. We observe a general trade-off between accuracy and the degrees of freedom of each model, and consequently computational cost. We will discuss these trade-offs in the context of model selection for molecular dynamics and other applications. △ Less

Submitted 24 July, 2019; v1 submitted 20 June, 2019; originally announced June 2019.

Data-driven Material Models for Atomistic Simulation

Authors: Mitchell A. Wood, Mary Alice Cusentino, Brian D. Wirth, Aidan P. Thompson

Abstract: The central approximation made in classical molecular dynamics simulation of materials is the interatomic potential used to calculate the forces on the atoms. Great effort and ingenuity is required to construct viable functional forms and find accurate parameterizations for potentials using traditional approaches. Machine-learning has emerged as an effective alternative approach to develop accurat… ▽ More The central approximation made in classical molecular dynamics simulation of materials is the interatomic potential used to calculate the forces on the atoms. Great effort and ingenuity is required to construct viable functional forms and find accurate parameterizations for potentials using traditional approaches. Machine-learning has emerged as an effective alternative approach to develop accurate and robust interatomic potentials. Starting with a very general model form, the potential is learned directly from a database of electronic structure calculations and therefore can be viewed as a multiscale link between quantum and classical atomistic simulations. Risk of inaccurate extrapolation exists outside the narrow range of time- and length-scales where the two methods can be directly compared. In this work, we use the Spectral Neighbor Analysis Potential (SNAP) and show how a fit can be produced with minimal interpolation errors which is also robust in extrapolating beyond training. To demonstrate the method, we have developed a new tungsten-beryllium potential suitable for the full range of binary compositions. Subsequently, large-scale molecular dynamics simulations were performed of high energy Be atom implantation onto the (001) surface of solid tungsten. The new machine learned W-Be potential generates a population of implantation structures consistent with quantum calculations of defect formation energies. A very shallow (<2nm) average Be implantation depth is predicted which may explain ITER diverter degradation in the presence of beryllium. △ Less

Submitted 18 March, 2019; v1 submitted 25 February, 2019; originally announced February 2019.

Comments: 5 Figures, 2 Tables

Journal ref: Phys. Rev. B 99, 184305 (2019)

Diamond nano-pillar arrays for quantum microscopy of neuronal signals

Authors: Liam Hanlon, Vini Gautam, James D. A. Wood, Prithvi Reddy, Michael S. J. Barson, Marika Niihori, Alexander R. J. Silalahi, Ben Corry, Joerg Wrachtrup, Matthew J. Sellars, Vincent R. Daria, Patrick Maletinsky, Gregory J. Stuart, Marcus W. Doherty

Abstract: Modern neuroscience is currently limited in its capacity to perform long term, wide-field measurements of neuron electromagnetics with nanoscale resolution. Quantum microscopy using the nitrogen vacancy centre (NV) can provide a potential solution to this problem with electric and magnetic field sensing at nano-scale resolution and good biocompatibility. However, the performance of existing NV sen… ▽ More Modern neuroscience is currently limited in its capacity to perform long term, wide-field measurements of neuron electromagnetics with nanoscale resolution. Quantum microscopy using the nitrogen vacancy centre (NV) can provide a potential solution to this problem with electric and magnetic field sensing at nano-scale resolution and good biocompatibility. However, the performance of existing NV sensing technology does not allow for studies of small mammalian neurons yet. In this paper, we propose a solution to this problem by engineering NV quantum sensors in diamond nanopillar arrays. The pillars improve light collection efficiency by guiding excitation/emission light, which improves sensitivity. More importantly, they also improve the size of the signal at the NV by removing screening charges as well as coordinating the neuron growth to the tips of the pillars where the NV is located. Here, we provide a growth study to demonstrate coordinated neuron growth as well as the first simulation of nano-scopic neuron electric and magnetic fields to assess the enhancement provided by the nanopillar geometry. △ Less

Submitted 25 January, 2019; originally announced January 2019.

Comments: 18 pages including supplementary and references, 12 figures

Nanoparticles Binding to Lipid Membranes: from Vesicle-Based Gels to Vesicle Inversion and Destruction

Authors: Derek A. Wood, Sarah Zuraw-Weston, Ian K. Torres, YiWei Lee, Li-Sheng Wang, Ziwen Jiang, Guillermo R. Lázaro, ShiYu Wang, Avital A. Rodal, Michael F. Hagan, Vincent M. Rotello, Anthony D. Dinsmore

Abstract: Cells offer numerous inspiring examples where proteins and membranes combine to form complex structures that are key to intracellular compartmentalization, cargo transport, and specialization of cell morphology. Despite this wealth of examples, we still lack the design principles to control membrane morphology in synthetic systems. Here we show that even the relatively simple case of spherical nan… ▽ More Cells offer numerous inspiring examples where proteins and membranes combine to form complex structures that are key to intracellular compartmentalization, cargo transport, and specialization of cell morphology. Despite this wealth of examples, we still lack the design principles to control membrane morphology in synthetic systems. Here we show that even the relatively simple case of spherical nanoparticles binding to lipid-bilayer membrane vesicles results in a remarkably rich set of morphologies that can be controlled quantitatively via the particle binding energy. We find that when the binding energy is weak relative to a characteristic membrane-bending energy, the vesicles adhere to one another and form a soft solid, which could be used as a useful platform for controlled release. When the binding energy is larger, the vesicles undergo a remarkable destruction process consisting first of invaginated tubules, followed by vesicles turning inside-out, yielding a network of nanoparticle-membrane tubules. We propose that the crossover from one behavior to the other is triggered by the transition from partial to complete wrapping of nanoparticles. This model is confirmed by computer simulations and by quantitative estimates of the binding energy. These findings open the door to a new class of vesicle-based, closed-cell gels that are more than 99% water and can encapsulate and release on demand. Our results also show how to intentionally drive dramatic shape changes in vesicles as a step toward shape-responsive particles. Finally, they help us to unify the wide range of previously observed responses of vesicles and cells to added nanoparticles. △ Less

Submitted 10 May, 2018; originally announced May 2018.

Comments: 13 pages, 8 figures, and supplemental materials

Liquid-liquid displacement in slippery liquid-infused membranes (SLIMs)

Authors: Hanieh Bazyar, Pengyu Lv, Jeffery A. Wood, Slawomir Porada, Detlef Lohse, Rob G. H. Lammertink

Abstract: Liquid-infused membranes inspired by slippery liquid-infused porous surfaces (SLIPS) have been recently introduced to membrane technology. The gating mechanism of these membranes is expected to give rise to anti-fouling properties and multi-phase transport capabilities. However, the long-term retention of the infusion liquid has not yet been explored. To address this issue, we investigate the rete… ▽ More Liquid-infused membranes inspired by slippery liquid-infused porous surfaces (SLIPS) have been recently introduced to membrane technology. The gating mechanism of these membranes is expected to give rise to anti-fouling properties and multi-phase transport capabilities. However, the long-term retention of the infusion liquid has not yet been explored. To address this issue, we investigate the retention of the infusion liquid in slippery liquid-infused membranes (SLIMs) via liquid-liquid displacement porometry (LLDP) experiments combined with microscopic observations of the displacement mechanism. Our results reveal that pores will be opened corresponding to the capillary pressure, leading to preferential flow pathways for water transport. The LLDP results further suggest the presence of liquid-lined pores in SLIM. This hypothesis is analyzed theoretically using an interfacial pore flow model. We find that the displacement patterns correspond to capillary fingering in immiscible displacement in porous media. The related physics regarding two-phase flow in porous media is used to confirm the permeation mechanism appearing in SLIMs. In order to experimentally observe liquid-liquid displacement, a microfluidic chip mimicking a porous medium is designed and a highly ramified structure with trapped infusion liquid is observed. The remaining infusion liquid is retained as pools, bridges and thin films around pillar structures in the chip, which further confirms liquid-lining. Fractal dimension analysis, along with evaluation of the fluid (non-wetting phase) saturation, further confirms that the fractal patterns correspond to capillary fingering, which is consistent with an invasion percolation with trapping (IPT) model. △ Less

Submitted 11 April, 2018; originally announced April 2018.

Quantum-correlated photons from semiconductor cavity polaritons

Authors: Guillermo Muñoz-Matutano, Andrew Wood, Mattias Johnson, Xavier Vidal Asensio, Ben Baragiola, Andreas Reinhard, Aristide Lemaitre, Jaqueline Bloch, Alberto Amo, Benjamin Besga, Maxime Richard, Thomas Volz

Abstract: Over the past decade, exciton-polaritons in semiconductor microcavities have attracted a great deal of interest as a driven-dissipative quantum fluid. These systems offer themselves as a versatile platform for performing Hamiltonian simulations with light, as well as for experimentally realizing nontrivial out-of-equilibrium phase transitions. In addition, polaritons exhibit a sizeable mutual inte… ▽ More Over the past decade, exciton-polaritons in semiconductor microcavities have attracted a great deal of interest as a driven-dissipative quantum fluid. These systems offer themselves as a versatile platform for performing Hamiltonian simulations with light, as well as for experimentally realizing nontrivial out-of-equilibrium phase transitions. In addition, polaritons exhibit a sizeable mutual interaction strength that opens up a whole range of possibilities in the context of quantum state generation. While squeezed light emission from polaritons has been reported previously, the granular nature of polaritons has not been observed to date. The latter capability is particularly attractive for realizing strongly correlated many-body quantum states of light on scalable arrays of coupled cavities. Here we demonstrate that by optically confining polaritons to a very small effective mode volume, one can reach the weak blockade regime, in which the nonlinearity turns strong enough to become significant at the few particle level, and thus produce a non-negligible antibunching in the emitted photons statistics. Our results act as a door opener for accessing the newly emerging field of quantum polaritonics, and as a proof of principle that optically confined exciton-polaritons can be considered as a realistic, new strategy to generate single photons. △ Less

Submitted 15 December, 2017; originally announced December 2017.

Comments: 12 pages, 10 figures

Journal ref: Nature Materials 18, 213 (2019)

Extending the Accuracy of the SNAP Interatomic Potential Form

Authors: Mitchell A. Wood, Aidan P. Thompson

Abstract: The Spectral Neighbor Analysis Potential (SNAP) is a classical interatomic potential that expresses the energy of each atom as a linear function of selected bispectrum components of the neighbor atoms. An extension of the SNAP form is proposed that includes quadratic terms in the bispectrum components. The extension is shown to provide a large increase in accuracy relative to the linear form, whil… ▽ More The Spectral Neighbor Analysis Potential (SNAP) is a classical interatomic potential that expresses the energy of each atom as a linear function of selected bispectrum components of the neighbor atoms. An extension of the SNAP form is proposed that includes quadratic terms in the bispectrum components. The extension is shown to provide a large increase in accuracy relative to the linear form, while incurring only a modest increase in computational cost. The mathematical structure of the quadratic SNAP form is similar to the embedded atom method (EAM), with the SNAP bispectrum components serving as counterparts to the two-body density functions in EAM. The effectiveness of the new form is demonstrated using an extensive set of training data for tantalum structures. Similarly to artificial neural network potentials, the quadratic SNAP form requires substantially more training data in order to prevent overfitting. The quality of this new potential form is measured through a robust cross-validation analysis. △ Less

Submitted 20 February, 2018; v1 submitted 29 November, 2017; originally announced November 2017.

Multiscale Modeling of Shock Wave Localization in Porous Energetic Material

Authors: Mitchell A. Wood, David E. Kittell, Cole D. Yarrington, Aidan P. Thompson

Abstract: Shock wave interactions with defects, such as pores, are known to play a key role in the chemical initiation of energetic materials. The shock response of hexanitrostilbene is studied through a combination of large scale reactive molecular dynamics and mesoscale hydrodynamic simulations. In order to extend our simulation capability at the mesoscale to include weak shock conditions (< 6 GPa), atomi… ▽ More Shock wave interactions with defects, such as pores, are known to play a key role in the chemical initiation of energetic materials. The shock response of hexanitrostilbene is studied through a combination of large scale reactive molecular dynamics and mesoscale hydrodynamic simulations. In order to extend our simulation capability at the mesoscale to include weak shock conditions (< 6 GPa), atomistic simulations of pore collapse are used to define a strain rate dependent strength model. Comparing these simulation methods allows us to impose physically-reasonable constraints on the mesoscale model parameters. In doing so, we have been able to study shock waves interacting with pores as a function of this viscoplastic material response. We find that the pore collapse behavior of weak shocks is characteristically different to that of strong shocks. △ Less

Submitted 7 November, 2017; originally announced November 2017.

Comments: 9 Pages, 7 Figures + 2 pages, 2 Figures Supplemental

Report number: SAND2017-12295 J

Journal ref: Phys. Rev. B 97, 014109 (2018)

Quantum probe hyperpolarisation of molecular nuclear spins

Authors: David A. Broadway, Jean-Philippe Tetienne, Alastair Stacey, James D. A. Wood, David A. Simpson, Liam T. Hall, Lloyd C. L. Hollenberg

Abstract: The hyperpolarisation of nuclear spins within target molecules is a critical and complex challenge in magnetic resonance imaging (MRI) and nuclear magnetic resonance (NMR) spectroscopy. Hyperpolarisation offers enormous gains in signal and spatial resolution which may ultimately lead to the development of molecular MRI and NMR. At present, techniques used to polarise nuclear spins generally requir… ▽ More The hyperpolarisation of nuclear spins within target molecules is a critical and complex challenge in magnetic resonance imaging (MRI) and nuclear magnetic resonance (NMR) spectroscopy. Hyperpolarisation offers enormous gains in signal and spatial resolution which may ultimately lead to the development of molecular MRI and NMR. At present, techniques used to polarise nuclear spins generally require low temperatures and/or high magnetic fields, radio-frequency control fields, or the introduction of catalysts or free-radical mediators. The emergence of room temperature solid-state spin qubits has opened exciting new pathways to circumvent these requirements to achieve direct nuclear spin hyperpolarisation using quantum control. Employing a novel cross-relaxation induced polarisation (CRIP) protocol using a single nitrogen-vacancy (NV) centre in diamond, we demonstrate the first external nuclear spin hyperpolarisation achieved by a quantum probe, in this case of $^1$H molecular spins in poly(methyl methacrylate). In doing so, we show that a single qubit is capable of increasing the thermal polarisation of $\sim 10^6$ nuclear spins by six orders of magnitude, equivalent to an applied magnetic field of $10^5$\,T. The technique can also be tuned to multiple spin species, which we demonstrate using both \C{13} and $^1$H nuclear spin ensembles. Our results are analysed and interpreted via a detailed theoretical treatment, which is also used to describe how the system can be scaled up to a universal quantum hyperpolarisation platform for the production of macroscopic quantities of contrast agents at high polarisation levels for clinical applications. These results represent a new paradigm for nuclear spin hyperpolarisation for molecular imaging and spectroscopy, and beyond into areas such as materials science and quantum information processing. △ Less

Submitted 22 August, 2017; v1 submitted 19 August, 2017; originally announced August 2017.

Comments: 6 pages, 4 figures

Renyi entropy of the totally asymmetric exclusion process

Authors: Anthony J. Wood, Richard A. Blythe, Martin R. Evans

Abstract: The Renyi entropy is a generalisation of the Shannon entropy that is sensitive to the fine details of a probability distribution. We present results for the Renyi entropy of the totally asymmetric exclusion process (TASEP). We calculate explicitly an entropy whereby the squares of configuration probabilities are summed, using the matrix product formalism to map the problem to one involving a six d… ▽ More The Renyi entropy is a generalisation of the Shannon entropy that is sensitive to the fine details of a probability distribution. We present results for the Renyi entropy of the totally asymmetric exclusion process (TASEP). We calculate explicitly an entropy whereby the squares of configuration probabilities are summed, using the matrix product formalism to map the problem to one involving a six direction lattice walk in the upper quarter plane. We derive the generating function across the whole phase diagram, using an obstinate kernel method. This gives the leading behaviour of the Renyi entropy and corrections in all phases of the TASEP. The leading behaviour is given by the result for a Bernoulli measure and we conjecture that this holds for all Renyi entropies. Within the maximal current phase the correction to the leading behaviour is logarithmic in the system size. Finally, we remark upon a special property of equilibrium systems whereby discontinuities in the Renyi entropy arise away from phase transitions, which we refer to as secondary transitions. We find no such secondary transition for this nonequilibrium system, supporting the notion that these are specific to equilibrium cases. △ Less

Submitted 28 September, 2017; v1 submitted 1 August, 2017; originally announced August 2017.

Journal ref: A J Wood, R A Blythe and M R Evans (2017) J. Phys. A: Math. Theor. 50 475005

Magnetic pseudo-fields in a rotating electron-nuclear spin system

Authors: A. A. Wood, E. Lilette, Y. Y. Fein, V. S. Perunicic, L. C. L. Hollenberg, R. E. Scholten, A. M. Martin

Abstract: A precessing spin observed in a rotating frame of reference appears frequency-shifted, an effect analogous to the precession of a Foucault pendulum observed on the rotating Earth. This frequency shift can be understood as arising from a magnetic pseudo-field in the rotating frame that nevertheless has physically significant consequences, such as the Barnett effect. Detecting these pseudo-fields is… ▽ More A precessing spin observed in a rotating frame of reference appears frequency-shifted, an effect analogous to the precession of a Foucault pendulum observed on the rotating Earth. This frequency shift can be understood as arising from a magnetic pseudo-field in the rotating frame that nevertheless has physically significant consequences, such as the Barnett effect. Detecting these pseudo-fields is experimentally challenging, as a rotating-frame sensor is required. Previous work has realised classical rotating-frame detectors. Here we use quantum sensors, nitrogen-vacancy (NV) centres, in a rapidly rotating diamond to detect pseudo-fields in the rotating frame. While conventional magnetic fields induce precession at a rate proportional to the gyromagnetic ratio, rotation shifts the precession of all spins equally, and thus primarily affect nearby $^{13}$C nuclear spins. We are thus able to explore these effects via quantum sensing in a rapidly rotating frame, and define a new approach to quantum control using rotationally-induced nuclear spin-selective magnetic fields. This work provides an integral step towards realising precision rotation sensing and quantum spin gyroscopes. △ Less

Submitted 11 July, 2017; originally announced July 2017.

Comments: 12 pages including 3 figures

Journal ref: Nature Physics 13, 1070 (2017)

Quantum-Accurate Molecular Dynamics Potential for Tungsten

Authors: Mitchell A. Wood, Aidan P. Thompson

Abstract: The purpose of this short contribution is to report on the development of a Spectral Neighbor Analysis Potential (SNAP) for tungsten. We have focused on the characterization of elastic and defect properties of the pure material in order to support molecular dynamics simulations of plasma-facing materials in fusion reactors. A parallel genetic algorithm approach was used to efficiently search for f… ▽ More The purpose of this short contribution is to report on the development of a Spectral Neighbor Analysis Potential (SNAP) for tungsten. We have focused on the characterization of elastic and defect properties of the pure material in order to support molecular dynamics simulations of plasma-facing materials in fusion reactors. A parallel genetic algorithm approach was used to efficiently search for fitting parameters optimized against a large number of objective functions. In addition, we have shown that this many-body tungsten potential can be used in conjunction with a simple helium pair potential to produce accurate defect formation energies for the W-He binary system. △ Less

Submitted 22 February, 2017; originally announced February 2017.

Comments: 7 Pages, 3 Figures, 3 Tables

Environmentally mediated coherent control of a spin qubit in diamond

Authors: Scott E. Lillie, David A. Broadway, James D. A. Wood, David A. Simpson, Alastair Stacey, Jean-Philippe Tetienne, Lloyd C. L. Hollenberg

Abstract: The coherent control of spin qubits forms the basis of many applications in quantum information processing and nanoscale sensing, imaging and spectroscopy. Such control is conventionally achieved by direct driving of the qubit transition with a resonant global field, typically at microwave frequencies. Here we introduce an approach that relies on the resonant driving of nearby environment spins, w… ▽ More The coherent control of spin qubits forms the basis of many applications in quantum information processing and nanoscale sensing, imaging and spectroscopy. Such control is conventionally achieved by direct driving of the qubit transition with a resonant global field, typically at microwave frequencies. Here we introduce an approach that relies on the resonant driving of nearby environment spins, whose localised magnetic field in turn drives the qubit when the environmental spin Rabi frequency matches the qubit resonance. This concept of environmentally mediated resonance (EMR) is explored experimentally using a qubit based on a single nitrogen-vacancy (NV) centre in diamond, with nearby electronic spins serving as the environmental mediators. We demonstrate EMR driven coherent control of the NV spin-state, including the observation of Rabi oscillations, free induction decay, and spin-echo. This technique also provides a way to probe the nanoscale environment of spin qubits, which we illustrate by acquisition of electron spin resonance spectra of single NV centres in various settings. △ Less

Submitted 19 February, 2017; originally announced February 2017.

Comments: 12 pages including SI

Journal ref: Phys. Rev. Lett. 118, 167204 (2017)

Microwave-Free Nuclear Magnetic Resonance at Molecular Scales

Authors: James D. A. Wood, Jean-Philippe Tetienne, David A. Broadway, Liam T. Hall, David A. Simpson, Alastair Stacey, Lloyd C. L. Hollenberg

Abstract: The implementation of nuclear magnetic resonance (NMR) at the nanoscale is a major challenge, as conventional systems require relatively large ensembles of spins and limit resolution to mesoscopic scales. New approaches based on quantum spin probes, such as the nitrogen-vacancy (NV) centre in diamond, have recently achieved nano-NMR under ambient conditions. However, the measurement protocols requ… ▽ More The implementation of nuclear magnetic resonance (NMR) at the nanoscale is a major challenge, as conventional systems require relatively large ensembles of spins and limit resolution to mesoscopic scales. New approaches based on quantum spin probes, such as the nitrogen-vacancy (NV) centre in diamond, have recently achieved nano-NMR under ambient conditions. However, the measurement protocols require application of complex microwave pulse sequences of high precision and relatively high power, placing limitations on the design and scalability of these techniques. Here we demonstrate a microwave-free method for nanoscale NMR using the NV centre, which is a far less invasive, and vastly simpler measurement protocol. By utilising a carefully tuned magnetic cross-relaxation interaction between a subsurface NV spin and an external, organic environment of proton spins, we demonstrate NMR spectroscopy of $^1$H within a $\approx(10~{\rm nm})^3$ sensing volume. We also theoretically and experimentally show that the sensitivity of our approach matches that of existing microwave control-based techniques using the NV centre. Removing the requirement for coherent manipulation of either the NV or the environmental spin quantum states represents a significant step towards the development of robust, non-invasive nanoscale NMR probes. △ Less

Submitted 6 October, 2016; originally announced October 2016.

Comments: 7 pages, 4 figures

Journal ref: Nature Communications 8:15950 (2017)

Measurement and extinction of vector light shifts using interferometry of spinor condensates

Authors: A. A. Wood, L. D. Turner, R. P. Anderson

Abstract: We use differential Ramsey interferometry of ultracold atoms to characterize the vector light shift (VLS) from a far-off resonance optical dipole trap at $λ= 1064\,\mathrm{nm}$. The VLS manifests as a `fictitious' magnetic field, which we perceive as a change in the Larmor frequency of two spinor condensates exposed to different intensities of elliptically polarized light. We use our measurement s… ▽ More We use differential Ramsey interferometry of ultracold atoms to characterize the vector light shift (VLS) from a far-off resonance optical dipole trap at $λ= 1064\,\mathrm{nm}$. The VLS manifests as a `fictitious' magnetic field, which we perceive as a change in the Larmor frequency of two spinor condensates exposed to different intensities of elliptically polarized light. We use our measurement scheme to diagnose the light-induced magnetic field and suppress it to $2.1(8)\times10^{-4}$ of its maximum value, by making the trapping light linearly polarized with a quarter-wave plate in each beam. Our sensitive measurement of the VLS-induced field demonstrates high-precision, in-vacuo interferometric polarimetry of dipole trapping light and can be adapted to measure vector shifts from other lasers, advancing the application of optically trapped atoms to precision metrology. △ Less

Submitted 20 November, 2016; v1 submitted 23 July, 2016; originally announced July 2016.

Comments: 12 pages, 6 figures

Journal ref: Phys. Rev. A 94, 052503 (2016)

Multiple Silicon Atom Artificial Molecules

Authors: John A. Wood, Mohammad Rashidi, Mohammad Koleini, Jason L. Pitters, Robert A. Wolkow

Abstract: We present linear ensembles of dangling bond chains on a hydrogen terminated Si(100) surface, patterned in the closest spaced arrangement allowed by the surface lattice. Local density of states maps over a range of voltages extending spatially over the close-coupled entities reveal a rich energetic and spatial variation of electronic states. These artificial molecules exhibit collective electronic… ▽ More We present linear ensembles of dangling bond chains on a hydrogen terminated Si(100) surface, patterned in the closest spaced arrangement allowed by the surface lattice. Local density of states maps over a range of voltages extending spatially over the close-coupled entities reveal a rich energetic and spatial variation of electronic states. These artificial molecules exhibit collective electronic states resulting from covalent interaction of the constituent atoms. A pronounced electrostatic perturbation of dangling bond chain structure is induced by close placement of a negatively dangling bond. The electronic changes so induced are entirely removed, paradoxically, by addition of a second dangling bond. △ Less

Submitted 20 July, 2016; originally announced July 2016.

Spin dynamics of diamond nitrogen-vacancy centres at the ground state level anti-crossing and all-optical low frequency magnetic field sensing

Authors: David A. Broadway, James D. A. Wood, Liam T. Hall, Alastair Stacey, Matthew Markham, David A. Simpson, Jean-Philippe Tetienne, Lloyd C. L. Hollenberg

Abstract: We investigate the photo-induced spin dynamics of single nitrogen-vacancy (NV) centres in diamond near the electronic ground state level anti-crossing (GSLAC), which occurs at an axial magnetic field around 1024 G. Using optically detected magnetic resonance spectroscopy, we first find that the electron spin transition frequency can be tuned down to 100 kHz for the \NV{14} centre, while for the \N… ▽ More We investigate the photo-induced spin dynamics of single nitrogen-vacancy (NV) centres in diamond near the electronic ground state level anti-crossing (GSLAC), which occurs at an axial magnetic field around 1024 G. Using optically detected magnetic resonance spectroscopy, we first find that the electron spin transition frequency can be tuned down to 100 kHz for the \NV{14} centre, while for the \NV{15} centre the transition strength vanishes for frequencies below about 2 MHz owing to the GSLAC level structure. Using optical pulses to prepare and readout the spin state, we observe coherent spin oscillations at 1024 G for the \NV{14}, which originate from spin mixing induced by residual transverse magnetic fields. This effect is responsible for limiting the smallest observable transition frequency, which can span two orders of magnitude from 100 kHz to tens of MHz depending on the local magnetic noise. A similar feature is observed for the \NV{15} centre at 1024 G. As an application of these findings, we demonstrate all-optical detection and spectroscopy of externally-generated fluctuating magnetic fields at frequencies from 8 MHz down to 500 kHz, using a \NV{14} centre. Since the Larmor frequency of most nuclear spin species lies within this frequency range near the GSLAC, these results pave the way towards all-optical, nanoscale nuclear magnetic resonance spectroscopy, using longitudinal spin cross-relaxation. △ Less

Submitted 14 July, 2016; originally announced July 2016.

Comments: 10 pages, 6 figures

Journal ref: Phys. Rev. Applied 6, 064001 (2016)

Persistent superfluid flow arising from the He-McKellar-Wilkens effect in molecular dipolar condensates

Authors: A. A. Wood, B. H. J. McKellar, A. M. Martin

Abstract: We show that the He-McKellar-Wilkens effect can induce a persistent flow in a Bose-Einstein condensate of polar molecules confined in a toroidal trap, with the dipolar interaction mediated via an electric dipole moment. For Bose-Einstein condensates of atoms with a magnetic dipole moment, we show that although it is theoretically possible to induce persistent flow via the Aharonov-Casher effect, t… ▽ More We show that the He-McKellar-Wilkens effect can induce a persistent flow in a Bose-Einstein condensate of polar molecules confined in a toroidal trap, with the dipolar interaction mediated via an electric dipole moment. For Bose-Einstein condensates of atoms with a magnetic dipole moment, we show that although it is theoretically possible to induce persistent flow via the Aharonov-Casher effect, the strength of electric field required is prohibitive. We also outline an experimental geometry tailored specifically for observing the He-McKellar-Wilkens effect in toroidally-trapped condensates. △ Less

Submitted 21 June, 2016; v1 submitted 18 April, 2016; originally announced April 2016.

Comments: 5 pages 2 figures

Journal ref: Phys. Rev. Lett. 116, 250403 (2016)

Momentum sharing in imbalanced Fermi systems

Authors: O. Hen, M. Sargsian, L. B. Weinstein, E. Piasetzky, H. Hakobyan, D. W. Higinbotham, M. Braverman, W. K. Brooks, S. Gilad, K. P. Adhikari, J. Arrington, G. Asryan, H. Avakian, J. Ball, N. A. Baltzell, M. Battaglieri, A. Beck, S. May-Tal Beck, I. Bedlinskiy, W. Bertozzi, A. Biselli, V. D. Burkert, T. Cao, D. S. Carman, A. Celentano , et al. (116 additional authors not shown)

Abstract: The atomic nucleus is composed of two different kinds of fermions, protons and neutrons. If the protons and neutrons did not interact, the Pauli exclusion principle would force the majority fermions (usually neutrons) to have a higher average momentum. Our high-energy electron scattering measurements using 12C, 27Al, 56Fe and 208Pb targets show that, even in heavy neutron-rich nuclei, short-range… ▽ More The atomic nucleus is composed of two different kinds of fermions, protons and neutrons. If the protons and neutrons did not interact, the Pauli exclusion principle would force the majority fermions (usually neutrons) to have a higher average momentum. Our high-energy electron scattering measurements using 12C, 27Al, 56Fe and 208Pb targets show that, even in heavy neutron-rich nuclei, short-range interactions between the fermions form correlated high-momentum neutron-proton pairs. Thus, in neutron-rich nuclei, protons have a greater probability than neutrons to have momentum greater than the Fermi momentum. This finding has implications ranging from nuclear few body systems to neutron stars and may also be observable experimentally in two-spin state, ultra-cold atomic gas systems. △ Less

Submitted 29 November, 2014; originally announced December 2014.

Comments: Published in Science. 10 pages, 3 figures

Journal ref: Science 346, 614 (2014)

Magnetic tensor gradiometry using Ramsey interferometry of spinor condensates

Authors: A. A. Wood, L. M. Bennie, A. Duong, M. Jasperse, L. D. Turner, R. P. Anderson

Abstract: We have realized a magnetic tensor gradiometer by interferometrically measuring the relative phase between two spatially separated Bose-Einstein condensates (BECs). We perform simultaneous Ramsey interferometry of the proximate $^{87}$Rb spin-1 condensates in freefall and infer their relative Larmor phase -- and thus the differential magnetic field strength -- with a common-mode phase noise suppre… ▽ More We have realized a magnetic tensor gradiometer by interferometrically measuring the relative phase between two spatially separated Bose-Einstein condensates (BECs). We perform simultaneous Ramsey interferometry of the proximate $^{87}$Rb spin-1 condensates in freefall and infer their relative Larmor phase -- and thus the differential magnetic field strength -- with a common-mode phase noise suppression exceeding $50\,\mathrm{dB}$. By appropriately biasing the magnetic field and separating the BECs along orthogonal directions, we measure the magnetic field gradient tensor of ambient and applied magnetic fields with a nominal precision of $30\,\mathrm{μG\,cm^{-1}}$ and a sensor volume of $2\times10^{-5}\,\mathrm{mm}^3$. We predict a spin-projection noise limited magnetic energy resolution of order $\hbar$ for typical Zeeman coherence times of trapped condensates with this scheme, even with the low measurement duty cycle inherent to BEC experiments. △ Less

Submitted 18 May, 2015; v1 submitted 5 August, 2014; originally announced August 2014.

Comments: 7 pages, 4 figures

Journal ref: Phys. Rev. A 92, 053604 (2015)

Self-Assembly on a Cylinder: A Model System for Understanding the Constraint of Commensurability

Authors: D. A. Wood, C. D. Santangelo, A. D. Dinsmore

Abstract: A crystal lattice, when confined to the surface of a cylinder, must have a periodic structure that is commensurate with the cylinder circumference. This constraint can frustrate the system, leading to oblique crystal lattices or to structures with a chiral seam known as a "line slip" phase, neither of which are stable for isotropic particles in equilibrium on flat surfaces. In this study, we use m… ▽ More A crystal lattice, when confined to the surface of a cylinder, must have a periodic structure that is commensurate with the cylinder circumference. This constraint can frustrate the system, leading to oblique crystal lattices or to structures with a chiral seam known as a "line slip" phase, neither of which are stable for isotropic particles in equilibrium on flat surfaces. In this study, we use molecular dynamics simulations to find the steady-state structure of spherical particles with short-range repulsion and long-range attraction far below the melting temperature. We vary the range of attraction using the Lennard-Jones and Morse potentials and find that a shorter-range attraction favors the line-slip. We develop a simple model based only on geometry and bond energy to predict when the crystal or line-slip phases should appear, and find reasonable agreement with the simulations. The simplicity of this model allows us to understand the influence of the commensurability constraint, an understanding that might be extended into the more general problem of self-assembling particles in strongly confined spaces. △ Less

Submitted 2 July, 2013; originally announced July 2013.

Comments: 12 pages, 9 figures. Submitted for publication, 2013

Direct measurement of the effective charge in nonpolar suspensions by optical tracking of single particles

Authors: G. Seth Roberts, Tiffany A. Wood, William J. Frith, Paul Bartlett

Abstract: We demonstrate a novel technique for the measurement of the charge carried by a colloidal particle. The technique uses the phenomenon of the resonance of a particle held in an optical tweezers trap and driven by a sinusoidal electric field. The trapped particle forms a strongly damped harmonic oscillator whose fluctuations are a function of $γ$, the ratio of the root-mean square average of the e… ▽ More We demonstrate a novel technique for the measurement of the charge carried by a colloidal particle. The technique uses the phenomenon of the resonance of a particle held in an optical tweezers trap and driven by a sinusoidal electric field. The trapped particle forms a strongly damped harmonic oscillator whose fluctuations are a function of $γ$, the ratio of the root-mean square average of the electric and thermal forces on the particle. At low applied fields, where $γ\ll 1$, the particle is confined to the optical axis while at high fields ($γ\gg 1$) the probability distribution of the particle is double-peaked. The periodically-modulated thermal fluctuations are measured with nanometer sensitivity using an interferometric position detector. Charges, as low as a few elementary charges, can be measured with an uncertainty of about 0.25 $e$. This is significantly better than previous techniques and opens up new possibilities for the study of nonpolar suspensions. △ Less

Submitted 26 January, 2007; originally announced January 2007.

Computer Algebra calculation of XRMS polarisation dependence in the non spherical case

Authors: Thomas A. Wood, Alessandro Mirone

Abstract: Simple analytical formulae, directly relating the experimental geometry and sample orientation to the measured R(M)XS scattered intensity are very useful to design experiments and analyse data. Such formulae can be obtained by the contraction of an expression containing the polarisations and crystal field tensors, and where the magnetisation vector acts as a rotation derivative\cite{mirone}. The… ▽ More Simple analytical formulae, directly relating the experimental geometry and sample orientation to the measured R(M)XS scattered intensity are very useful to design experiments and analyse data. Such formulae can be obtained by the contraction of an expression containing the polarisations and crystal field tensors, and where the magnetisation vector acts as a rotation derivative\cite{mirone}. The result of a contraction contains a scalar product of (rotated) polarisation vectors and the crystal field axis. The contraction rules give rise to combinatorial algorithms which can be efficiently treated by computers. In this work we provide and discuss a concise Mathematica code along with a few example applications to non-centrosymmetric magnetic systems. △ Less

Submitted 13 November, 2006; originally announced November 2006.

Comments: 25 pages, 1 figure. Trainee stage report

Tunable dipolar magnetism in high-spin molecular clusters

Authors: M. Evangelisti, A. Candini, A. Ghirri, M. Affronte, G. W. Powell, I. A. Gass, P. A. Wood, S. Parsons, E. K. Brechin, D. Collison, S. L. Heath

Abstract: We report on the Fe17 high-spin molecular cluster and show that this system is an exemplification of nanostructured dipolar magnetism. Each Fe17 molecule, with spin S=35/2 and axial anisotropy as small as D=-0.02K, is the magnetic unit that can be chemically arranged in different packing crystals whilst preserving both spin ground-state and anisotropy. For every configuration, molecular spins ar… ▽ More We report on the Fe17 high-spin molecular cluster and show that this system is an exemplification of nanostructured dipolar magnetism. Each Fe17 molecule, with spin S=35/2 and axial anisotropy as small as D=-0.02K, is the magnetic unit that can be chemically arranged in different packing crystals whilst preserving both spin ground-state and anisotropy. For every configuration, molecular spins are correlated only by dipolar interactions. The ensuing interplay between dipolar energy and anisotropy gives rise to macroscopic behaviors ranging from superparamagnetism to long-range magnetic order at temperatures below 1K. △ Less

Submitted 20 September, 2006; v1 submitted 7 February, 2006; originally announced February 2006.

Comments: Replaced with version accepted for publication in Physical Review Letters

Journal ref: Phys. Rev. Lett. 97, 167202 (2006)