Showing 1–50 of 162 results for author: Archer, A

Dynamics of particle aggregation in dewetting films of complex liquids

Authors: J. Zhang, D. N. Sibley, D. Tseluiko, A. J. Archer

Abstract: We consider the dynamic wetting and dewetting processes of films and droplets of complex liquids on planar surfaces, focusing on the case of colloidal suspensions, where the particle interactions can be sufficiently attractive to cause agglomeration of the colloids within the film. This leads to an interesting array of dynamic behaviours within the liquid and of the liquid-air interface. Incorpora… ▽ More We consider the dynamic wetting and dewetting processes of films and droplets of complex liquids on planar surfaces, focusing on the case of colloidal suspensions, where the particle interactions can be sufficiently attractive to cause agglomeration of the colloids within the film. This leads to an interesting array of dynamic behaviours within the liquid and of the liquid-air interface. Incorporating concepts from thermodynamics and using the thin-film approximation, we construct a model consisting of a pair of coupled partial differential equations that represent the evolution of the liquid film and the effective colloidal height profiles. We determine the relevant phase behaviour of the uniform system, including finding associated binodal and spinodal curves, helping to uncover how the emerging behaviour depends on the particle interactions. Performing a linear stability analysis of our system enables us to identify parameter regimes where agglomerates form, which we independently confirm through numerical simulations and continuation of steady states, to construct bifurcation diagrams. We obtain various dynamics such as uniform colloidal profiles in an unstable situation evolving into agglomerates and thus elucidate the interplay between dewetting and particle aggregation in complex liquids on surfaces. △ Less

Submitted 30 April, 2024; originally announced April 2024.

Comments: 39 pages, 18 figures

Experimental and Theoretical Bulk Phase Diagram and Interfacial Tension of Ouzo

Authors: Andrew J. Archer, Benjamin D. Goddard, David N. Sibley, James T. Rawlings, Ross Broadhurst, Fouzia F. Ouali, David J. Fairhurst

Abstract: Ouzo is a well-known drink in Mediterranean countries, with ingredients water, alcohol and trans-anethole oil. The oil is insoluble in water, but completely soluble in alcohol, so when water is added to the spirit, the available alcohol is depleted and the mixture exhibits spontaneous emulsification. This process is commonly known as the louche or Ouzo effect. Although the phase boundaries of this… ▽ More Ouzo is a well-known drink in Mediterranean countries, with ingredients water, alcohol and trans-anethole oil. The oil is insoluble in water, but completely soluble in alcohol, so when water is added to the spirit, the available alcohol is depleted and the mixture exhibits spontaneous emulsification. This process is commonly known as the louche or Ouzo effect. Although the phase boundaries of this archetypal ternary mixture are well known, the properties of coexisting phases have not previously been studied. Here, we present a detailed experimental investigation into the phase behaviour, including tie-lines connecting coexisting phases, determination of the critical point (also called the plait point in ternary systems) and measurements of the surface tension and density for varying alcohol concentrations. Additionally, we present a theory for the thermodynamics and phase diagram of the system. With suitable selection of the interaction parameters, the theory captures nearly all features of the experimental work. This simple model can be used to determine both bulk and non-uniform (e.g. interfacial) properties, paving the way for a wide range of future applications of the model to ternary mixtures in general. We show how our accurate equilibrium phase diagram can be used to provide improved understanding of non-equilibrium phenomena. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: 15 pages, 9 figures

Dynamic Density Functional Theory with Inertia and Background Flow

Authors: Rory D. Mills-Williams, Benjamin D. Goddard, Andrew J. Archer

Abstract: We present dynamic density functional theory (DDFT) incorporating general inhomogeneous, incompressible, time dependent background flows and inertia, describing externally driven passive colloidal systems out of equilibrium. We start by considering the underlying nonequilibrium Langevin dynamics, including the effect of the local velocity of the surrounding liquid bath, to obtain the nonlinear, no… ▽ More We present dynamic density functional theory (DDFT) incorporating general inhomogeneous, incompressible, time dependent background flows and inertia, describing externally driven passive colloidal systems out of equilibrium. We start by considering the underlying nonequilibrium Langevin dynamics, including the effect of the local velocity of the surrounding liquid bath, to obtain the nonlinear, nonlocal partial differential equations governing the evolution of the (coarse--grained) density and velocity fields describing the dynamics of colloids. Additionally, we show both with heuristic arguments, and by numerical solution, that our equations and solutions agree with existing DDFTs in the overdamped (high friction) limit. We provide numerical solutions that model the flow of hard spheres, in both unbounded and confined domains, and compare to previously--derived DDFTs with and without the background flow. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: 13 pages, 3 figures

MSC Class: 82-10; 82C05; 35R09

Binding potential and wetting behaviour of binary liquid mixtures on surfaces

Authors: Mounirah Areshi, Dmitri Tseluiko, Uwe Thiele, Benjamin D. Goddard, Andrew J. Archer

Abstract: We present a theory for the interfacial wetting phase behaviour of binary liquid mixtures on rigid solid substrates, applicable to both miscible and immiscible mixtures. In particular, we calculate the binding potential as a function of the adsorptions, i.e. the excess amounts of each of the two liquids at the substrate. The binding potential fully describes the corresponding interfacial thermodyn… ▽ More We present a theory for the interfacial wetting phase behaviour of binary liquid mixtures on rigid solid substrates, applicable to both miscible and immiscible mixtures. In particular, we calculate the binding potential as a function of the adsorptions, i.e. the excess amounts of each of the two liquids at the substrate. The binding potential fully describes the corresponding interfacial thermodynamics. Our approach is based on classical density functional theory. Binary liquid mixtures can exhibit complex bulk phase behaviour, including both liquid-liquid and vapour-liquid phase separation, depending on the nature of the interactions between all the particles of the two different liquids, the temperature and the chemical potentials. Here we show that the interplay between the bulk phase behaviour of the mixture and the properties of the interactions with the substrate gives rise to a wide variety of interfacial phase behaviours, including mixing and demixing situations. We find situations where the final state is a coexistence of up to three different phases. We determine how the liquid density profiles close to the substrate change as the interaction parameters are varied and how these determine the form of the binding potential, which in certain cases can be a multi-valued function of the adsorptions. We also present profiles for sessile droplets of both miscible and immiscible binary liquids. △ Less

Submitted 9 January, 2024; originally announced January 2024.

Comments: 22 pages, 23 figures

Journal ref: Phys. Rev. E 109, 024801 (2024)

Load is not what you should balance: Introducing Prequal

Authors: Bartek Wydrowski, Robert Kleinberg, Stephen M. Rumble, Aaron Archer

Abstract: We present Prequal (Probing to Reduce Queuing and Latency), a load balancer for distributed multi-tenant systems. Prequal aims to minimize real-time request latency in the presence of heterogeneous server capacities and non-uniform, time-varying antagonist load. It actively probes server load to leverage the power-of-d-choices paradigm, extending it with asynchronous and reusable probes. Cutting a… ▽ More We present Prequal (Probing to Reduce Queuing and Latency), a load balancer for distributed multi-tenant systems. Prequal aims to minimize real-time request latency in the presence of heterogeneous server capacities and non-uniform, time-varying antagonist load. It actively probes server load to leverage the power-of-d-choices paradigm, extending it with asynchronous and reusable probes. Cutting against received wisdom, Prequal does not balance CPU load, but instead selects servers according to estimated latency and active requests-in-flight (RIF). We explore its major design features on a testbed system and evaluate it on YouTube, where it has been deployed for more than two years. Prequal has dramatically decreased tail latency, error rates, and resource use, enabling YouTube and other production systems at Google to run at much higher utilization. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Comments: Accepted to USENIX NSDI 2024

ACM Class: C.2.4

VERITAS contributions to the 38th International Cosmic Ray Conference

Authors: A. Acharyya, C. B. Adams, A. Archer, P. Bangale, J. T. Bartkoske, P. Batista, W. Benbow, J. L. Christiansen, A. J. Chromey, A. Duerr, M. Errando, Q. Feng, G. M. Foote, L. Fortson, A. Furniss, W. Hanlon, O. Hervet, C. E. Hinrichs, J. Hoang, J. Holder, Z. Hughes, T. B. Humensky, W. Jin, M. N. Johnson, M. Kertzman , et al. (39 additional authors not shown)

Abstract: Compilation of papers presented by the VERITAS Collaboration at the 38th International Cosmic Ray Conference (ICRC), held July 26 through August 3, 2023 in Nagoya, Japan. Compilation of papers presented by the VERITAS Collaboration at the 38th International Cosmic Ray Conference (ICRC), held July 26 through August 3, 2023 in Nagoya, Japan. △ Less

Submitted 12 December, 2023; originally announced December 2023.

Comments: html page. ICRC 2023, Nagoya, Japan

Stability of nanoparticle laden aerosol liquid droplets

Authors: A. J. Archer, B. D. Goddard, R. Roth

Abstract: We develop a model for the thermodynamics and evaporation dynamics of aerosol droplets of a liquid such as water, surrounded by the gas. When the temperature and the chemical potential (or equivalently the humidity) are such that the vapour phase is the thermodynamic equilibrium state, then of course droplets of the pure liquid evaporate over a relatively short time. However, if the droplets also… ▽ More We develop a model for the thermodynamics and evaporation dynamics of aerosol droplets of a liquid such as water, surrounded by the gas. When the temperature and the chemical potential (or equivalently the humidity) are such that the vapour phase is the thermodynamic equilibrium state, then of course droplets of the pure liquid evaporate over a relatively short time. However, if the droplets also contain nanoparticles or any other non-volatile solute, then the droplets can become thermodynamically stable. We show that the equilibrium droplet size depends strongly on the amount and solubility of the nanoparticles within, i.e. on the nature of the particle interactions with the liquid, and of course also on the vapour temperature and chemical potential. We develop a simple thermodynamic model for such droplets and compare predictions with results from a lattice density functional theory that takes as input the same particle interaction properties, finding very good agreement. We also use dynamical density functional theory to study the evaporation/condensation dynamics of liquid from/to droplets as they equilibrate with the vapour, thereby demonstrating droplet stability. △ Less

Submitted 22 November, 2023; originally announced November 2023.

Comments: 25 pages, 9 figures

Journal ref: J. Chem. Phys. 159, 194503 (2023)

Practical Performance Guarantees for Pipelined DNN Inference

Authors: Aaron Archer, Matthew Fahrbach, Kuikui Liu, Prakash Prabhu

Abstract: We optimize pipeline parallelism for deep neural network (DNN) inference by partitioning model graphs into $k$ stages and minimizing the running time of the bottleneck stage, including communication. We give practical and effective algorithms for this NP-hard problem, but our emphasis is on tackling the practitioner's dilemma of deciding when a solution is good enough. To this end, we design novel… ▽ More We optimize pipeline parallelism for deep neural network (DNN) inference by partitioning model graphs into $k$ stages and minimizing the running time of the bottleneck stage, including communication. We give practical and effective algorithms for this NP-hard problem, but our emphasis is on tackling the practitioner's dilemma of deciding when a solution is good enough. To this end, we design novel mixed-integer programming (MIP) relaxations for proving lower bounds. Applying these methods to a diverse testbed of 369 production models, for $k \in \{2, 4, 8, 16, 32, 64\}$, we empirically show that these lower bounds are strong enough to be useful in practice. Our lower bounds are substantially stronger than standard combinatorial bounds. For example, evaluated via geometric means across a production testbed with $k = 16$ pipeline stages, our MIP formulations raise the lower bound from 0.4598 to 0.9452, expressed as a fraction of the best partition found. In other words, our improved lower bounds close the optimality gap by a factor of 9.855x. △ Less

Submitted 4 June, 2024; v1 submitted 6 November, 2023; originally announced November 2023.

Comments: 17 pages, 5 figures

Multiwavelength Observations of the Blazar PKS 0735+178 in Spatial and Temporal Coincidence with an Astrophysical Neutrino Candidate IceCube-211208A

Authors: A. Acharyya, C. B. Adams, A. Archer, P. Bangale, J. T. Bartkoske, P. Batista, W. Benbow, A. Brill, J. H. Buckley, J. L. Christiansen, A. J. Chromey, M. Errando, A. Falcone, Q. Feng, G. M. Foote, L. Fortson, A. Furniss, G. Gallagher, W. Hanlon, D. Hanna, O. Hervet, C. E. Hinrichs, J. Hoang, J. Holder, T. B. Humensky , et al. (185 additional authors not shown)

Abstract: We report on multiwavelength target-of-opportunity observations of the blazar PKS 0735+178, located 2.2$^\circ$ away from the best-fit position of the IceCube neutrino event IceCube-211208A detected on December 8, 2021. The source was in a high-flux state in the optical, ultraviolet, X-ray, and GeV gamma-ray bands around the time of the neutrino event, exhibiting daily variability in the soft X-ra… ▽ More We report on multiwavelength target-of-opportunity observations of the blazar PKS 0735+178, located 2.2$^\circ$ away from the best-fit position of the IceCube neutrino event IceCube-211208A detected on December 8, 2021. The source was in a high-flux state in the optical, ultraviolet, X-ray, and GeV gamma-ray bands around the time of the neutrino event, exhibiting daily variability in the soft X-ray flux. The X-ray data from Swift-XRT and NuSTAR characterize the transition between the low-energy and high-energy components of the broadband spectral energy distribution (SED), and the gamma-ray data from Fermi -LAT, VERITAS, and H.E.S.S. require a spectral cut-off near 100 GeV. Both X-ray and gamma-ray measurements provide strong constraints on the leptonic and hadronic models. We analytically explore a synchrotron self-Compton model, an external Compton model, and a lepto-hadronic model. Models that are entirely based on internal photon fields face serious difficulties in matching the observed SED. The existence of an external photon field in the source would instead explain the observed gamma-ray spectral cut-off in both leptonic and lepto-hadronic models and allow a proton jet power that marginally agrees with the Eddington limit in the lepto-hadronic model. We show a numerical lepto-hadronic model with external target photons that reproduces the observed SED and is reasonably consistent with the neutrino event despite requiring a high jet power. △ Less

Submitted 30 June, 2023; originally announced June 2023.

Comments: 21 pages, 3 figures, accepted by ApJ

A VERITAS/Breakthrough Listen Search for Optical Technosignatures

Authors: Atreya Acharyya, Colin Adams, Avery Archer, Priyadarshini Bangale, Pedro Batista, Wystan Benbow, Aryeh Brill, M Capasso, Manel Errando, Abraham Falcone, Qi Feng, John Finley, Gregory Foote, Lucy Fortson, Amy Furniss, Sean Griffin, William Hanlon, David Hanna, Olivier Hervet, Claire Hinrichs, John Hoang, Jamie Holder, T. Humensky, Weidong Jin, Philip Kaaret , et al. (43 additional authors not shown)

Abstract: The Breakthrough Listen Initiative is conducting a program using multiple telescopes around the world to search for "technosignatures": artificial transmitters of extraterrestrial origin from beyond our solar system. The VERITAS Collaboration joined this program in 2018, and provides the capability to search for one particular technosignature: optical pulses of a few nanoseconds duration detectabl… ▽ More The Breakthrough Listen Initiative is conducting a program using multiple telescopes around the world to search for "technosignatures": artificial transmitters of extraterrestrial origin from beyond our solar system. The VERITAS Collaboration joined this program in 2018, and provides the capability to search for one particular technosignature: optical pulses of a few nanoseconds duration detectable over interstellar distances. We report here on the analysis and results of dedicated VERITAS observations of Breakthrough Listen targets conducted in 2019 and 2020 and of archival VERITAS data collected since 2012. Thirty hours of dedicated observations of 136 targets and 249 archival observations of 140 targets were analyzed and did not reveal any signals consistent with a technosignature. The results are used to place limits on the fraction of stars hosting transmitting civilizations. We also discuss the minimum-pulse sensitivity of our observations and present VERITAS observations of CALIOP: a space-based pulsed laser onboard the CALIPSO satellite. The detection of these pulses with VERITAS, using the analysis techniques developed for our technosignature search, allows a test of our analysis efficiency and serves as an important proof-of-principle. △ Less

Submitted 30 June, 2023; originally announced June 2023.

Comments: 15 pages, 7 figures

VERITAS discovery of very high energy gamma-ray emission from S3 1227+25 and multiwavelength observations

Authors: Atreya Acharyya, Colin Adams, Avery Archer, Priyadarshini Bangale, Wystan Benbow, Aryeh Brill, Jodi Christiansen, Alisha Chromey, Manel Errando, Abe Falcone, Qi Feng, John Finley, Gregory Foote, Lucy Fortson, Amy Furniss, Greg Gallagher, William Hanlon, David Hanna, Olivier Hervet, Claire Hinrichs, John Hoang, Jamie Holder, Weidong Jin, Madalyn Johnson, Philip Kaaret , et al. (46 additional authors not shown)

Abstract: We report the detection of very high energy gamma-ray emission from the blazar S3 1227+25 (VER J1230+253) with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). VERITAS observations of the source were triggered by the detection of a hard-spectrum GeV flare on May 15, 2015 with the Fermi-Large Area Telescope (LAT). A combined five-hour VERITAS exposure on May 16th and May 18th… ▽ More We report the detection of very high energy gamma-ray emission from the blazar S3 1227+25 (VER J1230+253) with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). VERITAS observations of the source were triggered by the detection of a hard-spectrum GeV flare on May 15, 2015 with the Fermi-Large Area Telescope (LAT). A combined five-hour VERITAS exposure on May 16th and May 18th resulted in a strong 13$σ$ detection with a differential photon spectral index, $Γ$ = 3.8 $\pm$ 0.4, and a flux level at 9% of the Crab Nebula above 120 GeV. This also triggered target of opportunity observations with Swift, optical photometry, polarimetry and radio measurements, also presented in this work, in addition to the VERITAS and Fermi-LAT data. A temporal analysis of the gamma-ray flux during this period finds evidence of a shortest variability timescale of $τ_{obs}$ = 6.2 $\pm$ 0.9 hours, indicating emission from compact regions within the jet, and the combined gamma-ray spectrum shows no strong evidence of a spectral cut-off. An investigation into correlations between the multiwavelength observations found evidence of optical and gamma-ray correlations, suggesting a single-zone model of emission. Finally, the multiwavelength spectral energy distribution is well described by a simple one-zone leptonic synchrotron self-Compton radiation model. △ Less

Submitted 4 May, 2023; originally announced May 2023.

Comments: 18 pages, 6 figures. Accepted for publication in the Astrophysical Journal (ApJ)

Moss-like growth of metal electrodes: On the role of competing faradaic reactions and fast-charging

Authors: J. X. Kent Zheng, Jiefu Yin, Tian Tang, Lynden A. Archer

Abstract: Uncontrolled crystal growth during electroreduction of reactive metals in liquid electrolytes produces porous, low-density, mossy metal deposits that grow primarily along the surface normal vector to a planar electrode substrate. The moss-like deposits are fragile and cause premature failure of batteries by chemical, physical, and mechanical pathways. Here we use electroanalytical Rotating-Disk El… ▽ More Uncontrolled crystal growth during electroreduction of reactive metals in liquid electrolytes produces porous, low-density, mossy metal deposits that grow primarily along the surface normal vector to a planar electrode substrate. The moss-like deposits are fragile and cause premature failure of batteries by chemical, physical, and mechanical pathways. Here we use electroanalytical Rotating-Disk Electrode (RDE) studies in a three-electrode electrochemical cell to elucidate the fundamental origin of moss-like growth of metals. We report that competing Faradaic reactions occurring on the electrode surface is the source of the phenomenon. On this basis, we conclude that a moss-like growth regime can be accessed during electrodeposition of any metal by subtle shifts in electrolyte chemistry and deposition rate. Specifically, for Zn, a metal that conventionally is not known to form moss-like electrodeposits, obvious moss-like deposition patterns emerge at low-current densities in strongly-alkaline electrolytes that undergo electroreduction to form an interphase on the electrodeposited Zn. Conversely, we find that under conditions where the rate of metal electroplating is large relative to that of other competing Faradaic reactions, it is possible to eliminate the mossy-like growth regime for Zn. Taken together, our findings open up a new approach for simultaneously achieving favorable metal deposition morphology and fast charging in next-generation batteries using metal anodes. △ Less

Submitted 21 January, 2023; originally announced February 2023.

Search for Ultraheavy Dark Matter from Observations of Dwarf Spheroidal Galaxies with VERITAS

Authors: A. Acharyya, A. Archer, P. Bangale, J. T. Bartkoske, P. Batista, M. Baumgart, W. Benbow, J. H. Buckley, A. Falcone, Q. Feng, J. P. Finley, G. M. Foote, L. Fortson, A. Furniss, G. Gallagher, W. F. Hanlon, O. Hervet, J. Hoang, J. Holder, T. B. Humensky, W. Jin, P. Kaaret, M. Kertzman, M. Kherlakian, D. Kieda , et al. (29 additional authors not shown)

Abstract: Dark matter is a key piece of the current cosmological scenario, with weakly interacting massive particles (WIMPs) a leading dark matter candidate. WIMPs have not been detected in their conventional parameter space (100 GeV $\lesssim M_χ \lesssim$ 100 TeV), a mass range accessible with current Imaging Atmospheric Cherenkov Telescopes. As ultraheavy dark matter (UHDM; $M_χ \gtrsim$ 100 TeV) has bee… ▽ More Dark matter is a key piece of the current cosmological scenario, with weakly interacting massive particles (WIMPs) a leading dark matter candidate. WIMPs have not been detected in their conventional parameter space (100 GeV $\lesssim M_χ \lesssim$ 100 TeV), a mass range accessible with current Imaging Atmospheric Cherenkov Telescopes. As ultraheavy dark matter (UHDM; $M_χ \gtrsim$ 100 TeV) has been suggested as an under-explored alternative to the WIMP paradigm, we search for an indirect dark matter annihilation signal in a higher mass range (up to 30 PeV) with the VERITAS gamma-ray observatory. With 216 hours of observations of four dwarf spheroidal galaxies, we perform an unbinned likelihood analysis. We find no evidence of a $γ$-ray signal from UHDM annihilation above the background fluctuation for any individual dwarf galaxy nor for a joint-fit analysis, and consequently constrain the velocity-weighted annihilation cross section of UHDM for dark matter particle masses between 1 TeV and 30 PeV. We additionally set constraints on the allowed radius of a composite UHDM particle. △ Less

Submitted 24 February, 2023; v1 submitted 17 February, 2023; originally announced February 2023.

Comments: 10 pages, 7 figures

VTSCat: The VERITAS Catalog of Gamma-Ray Observations

Authors: A. Acharyya, C. B. Adams, A. Archer, P. Bangale, J. T. Bartkoske, P. Batista, W. Benbow, J. H. Buckley, A. Brill, M. Capasso, J. L. Christiansen, A. J. Chromey, M. K. Daniel, M. Errando, A. Falcone, K. A Farrell, Q. Feng, J. P. Finley, G. M Foote, L. Fortson, A. Furniss, G. Gallagher, A. Gent, C. Giuri, O. Gueta , et al. (64 additional authors not shown)

Abstract: The ground-based gamma-ray observatory VERITAS (Very Energetic Radiation Imaging Telescope Array System) is sensitive to photons of astrophysical origin with energies in the range between $\approx 85$ GeV to $\approx 30$ TeV. The instrument consists of four 12-m diameter imaging Cherenkov telescopes operating at the Fred Lawrence Whipple Observatory (FLWO) in southern Arizona. VERITAS started four… ▽ More The ground-based gamma-ray observatory VERITAS (Very Energetic Radiation Imaging Telescope Array System) is sensitive to photons of astrophysical origin with energies in the range between $\approx 85$ GeV to $\approx 30$ TeV. The instrument consists of four 12-m diameter imaging Cherenkov telescopes operating at the Fred Lawrence Whipple Observatory (FLWO) in southern Arizona. VERITAS started four-telescope operations in 2007 and collects about 1100 hours of good-weather data per year. The VERITAS collaboration has published over 100 journal articles since 2008 reporting on gamma-ray observations of a large variety of objects: Galactic sources like supernova remnants, pulsar wind nebulae, and binary systems; extragalactic sources like star forming galaxies, dwarf-spheroidal galaxies, and highly-variable active galactic nuclei. This note presents VTSCat: the catalog of high-level data products from all VERITAS publications. △ Less

Submitted 13 January, 2023; v1 submitted 11 January, 2023; originally announced January 2023.

Comments: Version with corrected author list

Journal ref: Research Notes of the AAS, Volume 7, Number 1, 2023

Electroconvective flow in presence of polyethylene glycol oligomer additives

Authors: Arpita Sharma, Ankush Mukherjee, Alexander Warren, Shuo Jin, Gaojin Li, Donald L. Koch, Lynden A. Archer

Abstract: Metal electrodeposition in batteries is fundamentally unstable and affected by different instabilities depending on operating conditions and chemical composition. Particularly at high charging rates, a hydrodynamic instability called electroconvection sets in that aggravates the situation by creating non-uniform ion flux and preferential deposition at the electrode. Here, we experimentally investi… ▽ More Metal electrodeposition in batteries is fundamentally unstable and affected by different instabilities depending on operating conditions and chemical composition. Particularly at high charging rates, a hydrodynamic instability called electroconvection sets in that aggravates the situation by creating non-uniform ion flux and preferential deposition at the electrode. Here, we experimentally investigate how oligomer additives interact with the hydrodynamic instability at a cation selective interface. From electrochemical measurements and direct visualization experiments, we find that electroconvection is delayed and suppressed at all voltage in the presence of oligomers. Our results also reveal that it is important to consider the role of polymers at the interface, in addition to their bulk effects, to understand the stabilization effect and its mechanism. △ Less

Submitted 20 November, 2022; originally announced November 2022.

Suppression of electroconvection due to van der Waals attraction of polymer additives towards the metal electrode

Authors: Ankush Mukherjee, Lynden A. Archer, Donald L. Koch

Abstract: Electroconvection in rechargeable batteries enhances the growth of dendrites at the electrode surface. The addition of low molecular weight polymers to the electrolyte in batteries results in the formation of a thin layer of higher polymer concentration near the electrode. This is due to van der Waals forces of attraction between the metal electrode and the polymers dissolved in the electrolyte. T… ▽ More Electroconvection in rechargeable batteries enhances the growth of dendrites at the electrode surface. The addition of low molecular weight polymers to the electrolyte in batteries results in the formation of a thin layer of higher polymer concentration near the electrode. This is due to van der Waals forces of attraction between the metal electrode and the polymers dissolved in the electrolyte. The van der Waals forces act as a restoring body force on the electrolyte and oppose the growth of perturbations. Using linear stability analysis, we show that this force opposes electroconvective flow. This increases the critical voltage required for the onset of electroconvection. △ Less

Submitted 27 October, 2022; originally announced October 2022.

Understanding The Reversible Electrodeposition of Al in Low-Cost Room Temperature Molten Salts

Authors: Regina Garcia-Mendez, Jingxu Zheng, David C. Bock, Cherno Jaye, Daniel A. Fischer, Amy C. Marschilok, Kenneth J. Takeuchi, Esther S. Takeuchi, Lynden A. Archer

Abstract: Aluminum is the most earth-abundant metal, is trivalent, is inert in ambient humid air, and has a density approximately four-times that of lithium at room temperature. These attributes make it an attractive material for cost-effective, long-duration storage of electrical energy in batteries. Scientific discoveries in the past decade have established that secondary Al batteries can be created by pa… ▽ More Aluminum is the most earth-abundant metal, is trivalent, is inert in ambient humid air, and has a density approximately four-times that of lithium at room temperature. These attributes make it an attractive material for cost-effective, long-duration storage of electrical energy in batteries. Scientific discoveries in the past decade have established that secondary Al batteries can be created by paring an Al anode with a graphite or transition metal oxide cathode, in imidazolium-based, room-temperature ionic-liquid-Aluminum chloride electrolytes. Here we report findings from a systematic study that sheds light on the structural requirements, physicochemical, and transport properties of the ionic liquid electrolytes responsible for the high reversibility of Al battery anodes. We find that the most important interfacial and transport properties of these electrolytes can be achieved in other electrolytes, including ammonium-based molten salts that are available at costs as much as twenty-times lower than the ionic liquid-Aluminum chloride melt. High Al reversibility in ammonium- and imidazolium-based electrolytes is specifically shown to require a critical ratio of the solvation species, where Lewis acidity and beneficial interfacial reactions continuously etch the alumina resistive interfacial layer and form beneficial solid electrolyte interphase at the anode. We report further that successful development of new electrolyte families that support high Al anode reversibility also provides a good platform for detailed studies of the working mechanisms of an intercalation graphite cathode using X-ray absorption spectroscopy. Our findings therefore open new opportunities for developing simple, cost-effective, room-temperature Al batteries that enable long-duration electrical energy storage. △ Less

Submitted 21 October, 2022; originally announced October 2022.

Comments: 44 pages, 26 figures, 1 table

Report number: 101452

Journal ref: Cell Reports Physical Science 4,101452 June 21 2023

Stochastic transitions: Paths over higher energy barriers can dominate in the early stages

Authors: S. P. Fitzgerald, A. Bailey Hass, G. Díaz Leines, A. J. Archer

Abstract: The time evolution of many physical, chemical, and biological systems can be modelled by stochastic transitions between the minima of the potential energy surface describing the system of interest. We show that in cases where there are two (or more) possible pathways that the system can take, the time available for the transition to occur is crucially important. The well-known results of reaction… ▽ More The time evolution of many physical, chemical, and biological systems can be modelled by stochastic transitions between the minima of the potential energy surface describing the system of interest. We show that in cases where there are two (or more) possible pathways that the system can take, the time available for the transition to occur is crucially important. The well-known results of reaction rate theory for determining the rates of the transitions apply in the long-time limit. However, at short times, the system can instead choose to pass over higher energy barriers with much higher probability, as long as the distance to travel in phase space is shorter. We construct two simple models to illustrate this general phenomenon. We also present an extension of the gMAM algorithm of Vanden-Eijnden and Heymann [J. Chem. Phys. {\bf 128}, 061103 (2008)] to determine the most likely path at both short and long times. △ Less

Submitted 20 October, 2022; originally announced October 2022.

Comments: 7 pages, 5 Figures

Thermodynamics, formation dynamics and structural correlations in the bulk amorphous phase of the phase-field crystal model

Authors: Shaho Abdalla, Andrew J. Archer, László Gránásy, Gyula I. Tóth

Abstract: We investigate bulk thermodynamic and microscopic structural properties of amorphous solids in the framework of the phase-field crystal (PFC) model. These are metastable states with a non-uniform density distribution having no long-range order. From extensive numerical simulations we determine the distribution of free energy density values in varying size amorphous systems and also the point-to-se… ▽ More We investigate bulk thermodynamic and microscopic structural properties of amorphous solids in the framework of the phase-field crystal (PFC) model. These are metastable states with a non-uniform density distribution having no long-range order. From extensive numerical simulations we determine the distribution of free energy density values in varying size amorphous systems and also the point-to-set correlation length, which is the radius of the largest volume of amorphous one can take while still having the particle arrangements within the volume determined by the particle ordering at the surface of the chosen volume. We find that in the thermodynamic limit, the free energy density of the amorphous tends to a value that has a slight dependence on the initial state from which it was formed -- i.e.\ it has a formation history dependence. The amorphous phase is observed to form on both sides of the liquid linear-stability limit, showing that the liquid to amorphous transition is first order, with an associated finite free energy barrier when the liquid is metastable. In our simulations this is demonstrated when noise in the initial density distribution is used to induce nucleation events from the metastable liquid. Depending on the strength of the initial noise, we observe a variety of nucleation pathways, in agreement with previous results for the PFC model, and which show that amorphous precursor mediated multi-step crystal nucleation can occur in colloidal systems. △ Less

Submitted 12 September, 2022; originally announced September 2022.

Comments: 13 pages, 16 figures

Rectangle--triangle soft-matter quasicrystals with hexagonal symmetry

Authors: Andrew J. Archer, Tomonari Dotera, Alastair M. Rucklidge

Abstract: Aperiodic (quasicrystalline) tilings, such as Penrose's tiling, can be built up from e.g. kites and darts, squares and equilateral triangles, rhombi or shield shaped tiles and can have a variety of different symmetries. However, almost all quasicrystals occurring in soft-matter are of the dodecagonal type. Here, we investigate a class of aperiodic tilings with hexagonal symmetry that are based on… ▽ More Aperiodic (quasicrystalline) tilings, such as Penrose's tiling, can be built up from e.g. kites and darts, squares and equilateral triangles, rhombi or shield shaped tiles and can have a variety of different symmetries. However, almost all quasicrystals occurring in soft-matter are of the dodecagonal type. Here, we investigate a class of aperiodic tilings with hexagonal symmetry that are based on rectangles and two types of equilateral triangles. We show how to design soft-matter systems of particles interacting via pair potentials containing two length-scales that form aperiodic stable states with two different examples of rectangle--triangle tilings. One of these is the bronze-mean tiling, while the other is a generalization. Our work points to how more general (beyond dodecagonal) quasicrystals can be designed in soft-matter. △ Less

Submitted 3 August, 2022; originally announced August 2022.

Comments: 15 pages, 13 figures. Submitted to Physical Review E. The data associated with this paper are openly available from the University of Leeds Data Repository at https://doi.org/10.5518/1188

MSC Class: 35B36 52C23

Journal ref: Phys. Rev. E 106, 044602 (2022)

Smectic Layering: Landau theory for a complex-tensor order parameter

Authors: Jack Paget, Una Alberti, Marco G. Mazza, Andrew J. Archer, Tyler N. Shendruk

Abstract: Composed of microscopic layers that stack along one direction while maintaining fluid-like positional disorder within layers, smectics are excellent systems for exploring topology, defects and geometric memory in complex confining geometries. However, the coexistence of crystalline-like characteristics in one direction and fluid-like disorder within layers makes lamellar liquid crystals notoriousl… ▽ More Composed of microscopic layers that stack along one direction while maintaining fluid-like positional disorder within layers, smectics are excellent systems for exploring topology, defects and geometric memory in complex confining geometries. However, the coexistence of crystalline-like characteristics in one direction and fluid-like disorder within layers makes lamellar liquid crystals notoriously difficult to model - especially in the presence of defects and large distortions. Nematic properties of smectics can be comprehensively described by the Q-tensor but to capture the features of the smectic layering alone, we develop a phenomenological Landau theory for a complex-tensor order parameter E, which is capable of describing the local degree of lamellar ordering, layer displacement, and orientation of the layers. This theory can account for both parallel and perpendicular elastic contributions. In addition to resolving the potential ambiguities inherent to complex scalar order parameter models, this model reduces to previous employed models of simple smectics, and opens new possibilities for numerical studies on smectics possessing many defects, within complex geometries and under extreme confinement. △ Less

Submitted 9 July, 2022; originally announced July 2022.

Comments: 23 pages, 3 figures

The Sabatier principle for Battery Anodes: Chemical Kinetics and Reversible Electrodeposition at Heterointerfaces

Authors: Jingxu Zheng, Yue Deng, Wenzao Li, Jiefu Yin, Patrick J. West, Tian Tang, Xiao Tong, David C. Bock, Shuo Jin, Qing Zhao, Regina Garcia-Mendez, Kenneth J. Takeuchi, Esther S. Takeuchi, Amy C. Marschilok, Lynden A. Archer

Abstract: How surface chemistry influences reactions occurring thereupon has been a long-standing question of broad scientific and technological interest for centuries. Recently, it has re-emerged as a critical question in a subdiscipline of chemistry - electrochemistry at heterointerphases, where the answers have implications for both how, and in what forms, humanity stores the rising quantities of renewab… ▽ More How surface chemistry influences reactions occurring thereupon has been a long-standing question of broad scientific and technological interest for centuries. Recently, it has re-emerged as a critical question in a subdiscipline of chemistry - electrochemistry at heterointerphases, where the answers have implications for both how, and in what forms, humanity stores the rising quantities of renewable electric power generated from solar and wind installations world-wide. Here we consider the relation between the surface chemistry at such interphases and the reversibility of electrochemical transformations at a rechargeable battery electrode. Conventional wisdom holds that stronger chemical interaction between the metal deposits and electrode promotes reversibility. We report instead that a moderate strength of chemical interaction between the deposit and the substrate, neither too weak nor too strong, enables highest reversibility and stability of the plating/stripping redox processes at a battery anode. Analogous to the empirical Sabatier principle for chemical heterogeneous catalysis, our finding arises from the confluence of competing processes - one driven by electrochemistry and the other by chemical alloying. Based on experimental evaluation of metal plating/stripping systems in battery anodes of contemporary interest, we show that such knowledge provides a powerful tool for designing key materials in highly reversible electrochemical energy storage technologies based on earth-abundant, low-cost metals. △ Less

Submitted 25 September, 2022; v1 submitted 14 March, 2022; originally announced March 2022.

Comments: accepted at Science Advances, in press

Complex-tensor theory of simple smectics

Authors: Jack Paget, Marco G. Mazza, Andew J. Archer, Tyler N. Shendruk

Abstract: Smectic materials represent a unique state between fluids and solids, characterized by orientational and partial positional order, making them notoriously difficult to model, particularly in confining geometries. We propose a complex order parameter tensor to describe the local degree of lamellar ordering, layer displacement and orientation. The theory accounts for both dislocations and disclinati… ▽ More Smectic materials represent a unique state between fluids and solids, characterized by orientational and partial positional order, making them notoriously difficult to model, particularly in confining geometries. We propose a complex order parameter tensor to describe the local degree of lamellar ordering, layer displacement and orientation. The theory accounts for both dislocations and disclinations, as well as arrested configurations and colloid-induced local ordering. It considerably simplifies numerics, facilitating studies on the dynamics of topologically complex lamellar systems. △ Less

Submitted 22 January, 2022; originally announced January 2022.

Comments: 6 pages, 3 figures

VERITAS contributions to the 37th International Cosmic Ray Conference

Authors: C. B. Adams, A. Archer, W. Benbow, A. Brill, J. H. Buckley, M. Capasso, J. L. Christiansen, A. J. Chromey, M. Errando, A. Falcone, K. A. Farrell, Q. Feng, G. M. Foote, L. Fortson, A. Furniss, A. Gent, G. H. Gillanders, C. Giuri, O. Gueta, D. Hanna, O. Hervet, J. Holder, B. Hona, T. B. Humensky, W. Jin , et al. (36 additional authors not shown)

Abstract: Compilation of papers presented by the VERITAS Collaboration at the 37th International Cosmic Ray Conference (ICRC), held July 12 through July 23, 2021 (online) in Berlin, Germany. Compilation of papers presented by the VERITAS Collaboration at the 37th International Cosmic Ray Conference (ICRC), held July 12 through July 23, 2021 (online) in Berlin, Germany. △ Less

Submitted 10 September, 2021; originally announced September 2021.

Comments: html page. 2021 ICRC, Berlin, Germany

Dynamical density functional theory for the drying and stratification of binary colloidal dispersions

Authors: Boshen He, Ignacio Martin-Fabiani, Roland Roth, Gyula I. Tóth, Andrew J. Archer

Abstract: We develop a dynamical density functional theory based model for the drying of colloidal films on planar surfaces. We consider mixtures of two different sizes of hard-sphere colloids. Depending on the solvent evaporation rate and the initial concentrations of the two species, we observe varying degrees of stratification in the final dried films. Our model predicts the various structures described… ▽ More We develop a dynamical density functional theory based model for the drying of colloidal films on planar surfaces. We consider mixtures of two different sizes of hard-sphere colloids. Depending on the solvent evaporation rate and the initial concentrations of the two species, we observe varying degrees of stratification in the final dried films. Our model predicts the various structures described in the literature previously from experiments and computer simulations, in particular the small-on-top stratified films. Our model also includes the influence of adsorption of particles to the interfaces. △ Less

Submitted 27 January, 2021; originally announced January 2021.

Comments: 13 pages, 12 figures, published in Langmuir

Snaking without subcriticality: grain boundaries as non-topological defects

Authors: Priya Subramanian, Andrew J. Archer, Edgar Knobloch, Alastair M. Rucklidge

Abstract: Non-topological defects such as grain boundaries abound in pattern forming systems, arising from local variations of pattern properties such as amplitude, wavelength, orientation, etc. We introduce the idea of treating such non-topological defects as spatially localised structures that are embedded in a background pattern, instead of treating them in an amplitude-phase decomposition. Using the two… ▽ More Non-topological defects such as grain boundaries abound in pattern forming systems, arising from local variations of pattern properties such as amplitude, wavelength, orientation, etc. We introduce the idea of treating such non-topological defects as spatially localised structures that are embedded in a background pattern, instead of treating them in an amplitude-phase decomposition. Using the two-dimensional quadratic-cubic Swift--Hohenberg equation as an example we obtain fully nonlinear equilibria that contain grain boundaries which are closed curves containing multiple penta-hepta defects separating regions of hexagons with different orientations. These states arise from local orientation mismatch between two stable hexagon patterns, one of which forms the localised grain and the other its background, and do not require a subcritical bifurcation connecting them. Multiple robust isolas that span a wide range of parameters are obtained even in the absence of a unique Maxwell point, underlining the importance of retaining pinning when analysing patterns with defects, an effect omitted from the amplitude-phase description. △ Less

Submitted 17 November, 2020; originally announced November 2020.

Comments: 16 pages, 12 figures and 2 movies in mp4 format

Journal ref: IMA Journal of Applied Mathematics (2021) 86, 1164-1180

Density Distribution in Soft Matter Crystals and Quasicrystals

Authors: Priya Subramanian, Daniel J. Ratliff, Alastair M. Rucklidge, Andrew J. Archer

Abstract: The density distribution in solids is often represented as a sum of Gaussian peaks (or similar functions) centred on lattice sites or via a Fourier sum. Here, we argue that representing instead the $\mathit{logarithm}$ of the density distribution via a Fourier sum is better. We show that truncating such a representation after only a few terms can be highly accurate for soft matter crystals. For qu… ▽ More The density distribution in solids is often represented as a sum of Gaussian peaks (or similar functions) centred on lattice sites or via a Fourier sum. Here, we argue that representing instead the $\mathit{logarithm}$ of the density distribution via a Fourier sum is better. We show that truncating such a representation after only a few terms can be highly accurate for soft matter crystals. For quasicrystals, this sum does not truncate so easily, nonetheless, representing the density profile in this way is still of great use, enabling us to calculate the phase diagram for a 3-dimensional quasicrystal forming system using an accurate non-local density functional theory. △ Less

Submitted 27 May, 2021; v1 submitted 20 October, 2020; originally announced October 2020.

Comments: 6 pages, 4 figures, 1 supplementary information PDF, 1 MATLAB program, 1 MAT datafile

Journal ref: Phys. Rev. Lett. 126, 218003 (2021)

Localized states in passive and active phase-field-crystal models

Authors: Max Philipp Holl, Andrew J. Archer, Svetlana V. Gurevich, Edgar Knobloch, Lukas Ophaus, Uwe Thiele

Abstract: The passive conserved Swift-Hohenberg equation (or phase-field-crystal [PFC] model) corresponds to a gradient dynamics for a single order parameter field related to density. It provides a simple microscopic description of the thermodynamic transition between liquid and crystalline states. In addition to spatially extended periodic structures, the model describes a large variety of steady spatially… ▽ More The passive conserved Swift-Hohenberg equation (or phase-field-crystal [PFC] model) corresponds to a gradient dynamics for a single order parameter field related to density. It provides a simple microscopic description of the thermodynamic transition between liquid and crystalline states. In addition to spatially extended periodic structures, the model describes a large variety of steady spatially localized structures. In appropriate bifurcation diagrams the corresponding solution branches exhibit characteristic slanted homoclinic snaking. In an active PFC model, encoding for instance the active motion of self-propelled colloidal particles, the gradient dynamics structure is broken by a coupling between density and an additional polarization field. Then, resting and traveling localized states are found with transitions characterized by parity-breaking drift bifurcations. Here, we first briefly review the snaking behavior of localized states in passive and active PFC models before discussing the bifurcation behavior of localized states in systems of (i) two coupled passive PFC equations described by common gradient dynamics, (ii) two coupled passive PFC where the coupling breaks the gradient dynamics structure, and (iii) a passive PFC coupled to an active PFC. △ Less

Submitted 6 October, 2020; originally announced October 2020.

Comments: submitted to the IMA Journal of Applied Mathematics' Special Issue on Homoclinic Snaking at 21, in memory of Patrick Woods

Efficient calculation of phase coexistence and phase diagrams: application to a binary phase-field crystal model

Authors: Max Philipp Holl, Andrew J. Archer, Uwe Thiele

Abstract: We show that one can employ well-established numerical continuation methods to efficiently calculate the phase diagram for thermodynamic systems. In particular, this involves the determination of lines of phase coexistence related to first order phase transitions and the continuation of triple points. To illustrate the method we apply it to a binary Phase-Field-Crystal model for the crystallisatio… ▽ More We show that one can employ well-established numerical continuation methods to efficiently calculate the phase diagram for thermodynamic systems. In particular, this involves the determination of lines of phase coexistence related to first order phase transitions and the continuation of triple points. To illustrate the method we apply it to a binary Phase-Field-Crystal model for the crystallisation of a mixture of two types of particles. The resulting phase diagram is determined for one- and two-dimensional domains. In the former case it is compared to the diagram obtained from a one-mode approximation. The various observed liquid and crystalline phases and their stable and metastable coexistence are discussed as well as the temperature-dependence of the phase diagrams. This includes the (dis)appearance of critical points and triple points. We also relate bifurcation diagrams for finite-size systems to the thermodynamics of phase transitions in the infinite-size limit. △ Less

Submitted 7 September, 2020; originally announced September 2020.

Journal ref: J. Phys.: Condens. Matter 33, 115401 (2021)

Length-scales in sheared soft matter depend sensitively on molecular interactions

Authors: A. Scacchi, M. G. Mazza, A. J. Archer

Abstract: The structure and degree of order in soft matter and other materials is intimately connected to the nature of the interactions between the particles. One important research goal is to find suitable control mechanisms, to enhance or suppress different structures. Using dynamical density functional theory, we investigate the interplay between external shear and the characteristic length-scales in th… ▽ More The structure and degree of order in soft matter and other materials is intimately connected to the nature of the interactions between the particles. One important research goal is to find suitable control mechanisms, to enhance or suppress different structures. Using dynamical density functional theory, we investigate the interplay between external shear and the characteristic length-scales in the interparticle correlations of a model system. We show that shear can controllably change the characteristic length-scale from one to another quite distinct value. Moreover, with specific small changes in the form of the particle interactions, the applied shear can either selectively enhance or suppress the different characteristic wavelengths of the system, thus showing how to tune these. Our results suggest that the nonlinear response to flow can be harnessed to design novel actively responsive materials. △ Less

Submitted 21 August, 2020; originally announced August 2020.

Comments: 8 pages, 2 figures

Journal ref: Phys. Rev. Research 2, 032064 (2020)

Quasicrystal formation in binary soft matter mixtures

Authors: A. Scacchi, W. R. C. Somerville, D. M. A. Buzza, A. J. Archer

Abstract: Using a strategy that may be applied in theory or in experiments, we identify the regime in which a model binary soft matter mixture forms quasicrystals. The system is described using classical density functional theory combined with integral equation theory. Quasicrystal formation requires particle ordering with two characteristic lengthscales in certain particular ratios. How the lengthscales ar… ▽ More Using a strategy that may be applied in theory or in experiments, we identify the regime in which a model binary soft matter mixture forms quasicrystals. The system is described using classical density functional theory combined with integral equation theory. Quasicrystal formation requires particle ordering with two characteristic lengthscales in certain particular ratios. How the lengthscales are related to the form of the pair interactions is reasonably well understood for one component systems, but less is known for mixtures. In our model mixture of big and small colloids confined to an interface, the two lengthscales stem from the range of the interactions between pairs of big particles and from the cross big-small interactions, respectively. The small-small lengthscale is not significant. Our strategy for finding quasicrystals involves tuning locations of maxima in the dispersion relation, or equivalently in the liquid state partial static structure factors. △ Less

Submitted 24 July, 2020; originally announced July 2020.

Comments: 8 pages, 5 figures

Journal ref: Phys. Rev. Research 2, 032043(R) (2020)

Suppression of Electroconvective and Morphological Instabilities by an Imposed Cross Flow of the Electrolyte

Authors: Gaojin Li, Alex Townsend, Lynden A. Archer, Donald L. Koch

Abstract: Electroconvection and its coupling with a morphological instability are important in many applications, including electrodialysis, batteries and fuel cells. In this work, we study the effects of a two-dimensional channel flow on the electroconvective and morphological instabilities using two approaches. In the bulk analysis, we consider the instability of the electroneutral bulk region driven by a… ▽ More Electroconvection and its coupling with a morphological instability are important in many applications, including electrodialysis, batteries and fuel cells. In this work, we study the effects of a two-dimensional channel flow on the electroconvective and morphological instabilities using two approaches. In the bulk analysis, we consider the instability of the electroneutral bulk region driven by a second kind electroosmosis slip velocity boundary condition and derive the asymptotic solutions for small and large wavenumbers. In the full analysis, we consider the entire region of the liquid electrolyte and use the ultraspherical spectral method to numerically solve the eigenvalue problems. Both studies show that the imposed flow significantly affects the electroconvective instability. The imposed flow generates a shielding effect by deforming the perturbed ion concentration field and hinders the ion transfer from low- to high- concentration regions which causes the instability. It fully suppresses the electroconvective instability at small wavenumbers and reduces the growth rate of the perturbations at large wavenumbers. The direct effect of the flow on the morphological instability is minor, while the suppression of the electroconvective instability may change the wavenumber of the most unstable mode of the coupled instabilities. For the electroconvective instability, the bulk analysis is qualitatively different from the full analysis at high wavenumbers. For the morphological instability, good agreement is found between the two studies at both small and large wavenumbers. △ Less

Submitted 9 July, 2020; originally announced July 2020.

Journal ref: Phys. Rev. Fluids 6, 033701 (2021)

VERITAS Discovery of VHE Emission from the Radio Galaxy 3C 264: A Multi-Wavelength Study

Authors: A. Archer, W. Benbow, R. Bird, A. Brill, M. Buchovecky, J. H. Buckley, M. T. Carini, J. L. Christiansen, A. J. Chromey, M. K. Daniel, M. Errando, A. Falcone, Q. Feng, P. Fortin, L. Fortson, A. Furniss, A. Gent, M. Georganopoulos, G. H. Gillanders, C. Giuri, O. Gueta, D. Hanna, T. Hassan, O. Hervet, J. Holder , et al. (45 additional authors not shown)

Abstract: The radio source 3C 264, hosted by the giant elliptical galaxy NGC 3862, was observed with VERITAS between February 2017 and May 2019. These deep observations resulted in the discovery of very-high-energy (VHE; E $>100$ GeV) $γ$-ray emission from this active galaxy. An analysis of $\sim$57 hours of quality-selected live time yields a detection at the position of the source, corresponding to a stat… ▽ More The radio source 3C 264, hosted by the giant elliptical galaxy NGC 3862, was observed with VERITAS between February 2017 and May 2019. These deep observations resulted in the discovery of very-high-energy (VHE; E $>100$ GeV) $γ$-ray emission from this active galaxy. An analysis of $\sim$57 hours of quality-selected live time yields a detection at the position of the source, corresponding to a statistical significance of 7.8 standard deviations above background. The observed VHE flux is variable on monthly time scales, with an elevated flux seen in 2018 observations. The VHE emission during this elevated state is well-characterized by a power-law spectrum with a photon index $Γ= 2.20 \pm 0.27$ and flux F($>315$ GeV) = ($7.6\pm 1.2_{\mathrm stat} \pm 2.3_{\mathrm syst})\times 10^{-13}$ cm$^{-2}$ s$^{-1}$, or approximately 0.7% of the Crab Nebula flux above the same threshold. 3C 264 ($z = 0.0217$) is the most distant radio galaxy detected at VHE, and the elevated state is thought to be similar to that of the famously outbursting jet in M 87. Consequently, extensive contemporaneous multi-wavelength data were acquired in 2018 at the time of the VHE high state. An analysis of these data, including VLBA, VLA, HST, Chandra and Swift observations in addition to the VERITAS data, is presented, along with a discussion of the resulting spectral energy distribution. △ Less

Submitted 6 May, 2020; originally announced May 2020.

Comments: 19 pages, 11 figures, Accepted for publication in Astrophysical Journal

Evidence for proton acceleration up to TeV energies based on VERITAS and Fermi-LAT observations of the Cas A SNR

Authors: A. U. Abeysekara, A. Archer, W. Benbow, R. Bird, R. Brose, M. Buchovecky, J. H. Buckley, A. J. Chromey, W. Cui, M. K. Daniel, S. Das, V. V. Dwarkadas, A. Falcone, Q. Feng, J. P. Finley, L. Fortson, A. Gent, G. H. Gillanders, C. Giuri, O. Gueta, D. Hanna, T. Hassan, O. Hervet, J. Holder, G. Hughes , et al. (38 additional authors not shown)

Abstract: We present a study of $γ$-ray emission from the core-collapse supernova remnant Cas~A in the energy range from 0.1GeV to 10TeV. We used 65 hours of VERITAS data to cover 200 GeV - 10 TeV, and 10.8 years of \textit{Fermi}-LAT data to cover 0.1-500 GeV. The spectral analysis of \textit{Fermi}-LAT data shows a significant spectral curvature around $1.3 \pm 0.4_{stat}$ GeV that is consistent with the… ▽ More We present a study of $γ$-ray emission from the core-collapse supernova remnant Cas~A in the energy range from 0.1GeV to 10TeV. We used 65 hours of VERITAS data to cover 200 GeV - 10 TeV, and 10.8 years of \textit{Fermi}-LAT data to cover 0.1-500 GeV. The spectral analysis of \textit{Fermi}-LAT data shows a significant spectral curvature around $1.3 \pm 0.4_{stat}$ GeV that is consistent with the expected spectrum from pion decay. Above this energy, the joint spectrum from \textit{Fermi}-LAT and VERITAS deviates significantly from a simple power-law, and is best described by a power-law with spectral index of $2.17\pm 0.02_{stat}$ with a cut-off energy of $2.3 \pm 0.5_{stat}$ TeV. These results, along with radio, X-ray and $γ$-ray data, are interpreted in the context of leptonic and hadronic models. Assuming a one-zone model, we exclude a purely leptonic scenario and conclude that proton acceleration up to at least 6 TeV is required to explain the observed $γ$-ray spectrum. From modeling of the entire multi-wavelength spectrum, a minimum magnetic field inside the remnant of $B_{\mathrm{min}}\approx150\,\mathrm{μG}$ is deduced. △ Less

Submitted 30 March, 2020; originally announced March 2020.

Comments: 33 pages, 9 Figures, 6 Tables

Pattern formation in two-dimensional hard-core/soft-shell systems with variable soft shell profiles

Authors: Walter R. C. Somerville, Adam D. Law, Marcel Rey, Nicolas Vogel, Andrew J. Archer, D. Martin A. Buzza

Abstract: Hard-core/soft shell (HCSS) particles have been shown to self-assemble into a remarkably rich variety of structures under compression due to the simple interplay between the hard-core and soft-shoulder length scales in their interactions. Most studies in this area model the soft shell interaction as a square shoulder potential. Although appealing from a theoretical point of view, the potential is… ▽ More Hard-core/soft shell (HCSS) particles have been shown to self-assemble into a remarkably rich variety of structures under compression due to the simple interplay between the hard-core and soft-shoulder length scales in their interactions. Most studies in this area model the soft shell interaction as a square shoulder potential. Although appealing from a theoretical point of view, the potential is physically unrealistic because there is no repulsive force in the soft shell regime, unlike in experimental HCSS systems. To make the model more realistic, here we consider HCSS particles with a range soft shell potential profiles beyond the standard square shoulder form and study the model using both minimum energy calculations and Monte Carlo simulations. We find that by tuning density and the soft shell profile, HCSS particles in the thin shell regime (i.e., shell to core ratio $r_1/r_0 \leq \sqrt{3}$) can form a large range of structures, including hexagons, chains, squares, rhomboids and two distinct zig-zag structures. Furthermore, by tuning the density and $r_1/r_0$, we find that HCSS particles with experimentally realistic linear ramp soft shoulder repulsions can form honeycombs and quasicrystals with 10-fold and 12-fold symmetry. Our study therefore suggests the exciting possibility of fabricating these exotic 2D structures experimentally through colloidal self-assembly. △ Less

Submitted 10 March, 2020; originally announced March 2020.

Comments: 12 pages, 7 figures

Journal ref: Soft Matter 16, 3564 (2020)

How ice grows from premelting films and water droplets

Authors: David N. Sibley, Pablo Llombart, Eva G. Noya, Andrew J. Archer, Luis G. MacDowell

Abstract: Close to the triple point, the surface of ice is covered by a thin liquid layer (so-called quasi-liquid layer) which crucially impacts growth and melting rates. Experimental probes cannot observe the growth processes below this layer, and classical models of growth by vapor deposition do not account for the formation of premelting films. Here, we develop a mesoscopic model of liquid-film mediated… ▽ More Close to the triple point, the surface of ice is covered by a thin liquid layer (so-called quasi-liquid layer) which crucially impacts growth and melting rates. Experimental probes cannot observe the growth processes below this layer, and classical models of growth by vapor deposition do not account for the formation of premelting films. Here, we develop a mesoscopic model of liquid-film mediated ice growth, and identify the various resulting growth regimes. At low saturation, freezing proceeds by terrace spreading, but the motion of the buried solid is conveyed through the liquid to the outer liquid-vapor interface. At higher saturations water droplets condense, a large crater forms below, and freezing proceeds undetectably beneath the droplet. Our approach is a general framework that naturally models freezing close to three phase coexistence and provides a first principle theory of ice growth and melting which may prove useful in the geosciences. △ Less

Submitted 21 January, 2021; v1 submitted 9 March, 2020; originally announced March 2020.

Comments: 32 pages, 10 figures

Journal ref: Nature Comm. 12, 239 (2021)

Rounded Layering Transitions on the Surface of Ice

Authors: Pablo Llombart, Eva G. Noya, David N. Sibley, Andrew J. Archer, Luis G. MacDowell

Abstract: Understanding the wetting properties of premelting films requires knowledge of the film's equation of state, which is not usually available. Here we calculate the disjoining pressure curve of premelting films, and perform a detailed thermodynamic characterization of premelting behavior on ice. Analysis of the density profiles reveals the signature of weak layering phenomena, from one to two and fr… ▽ More Understanding the wetting properties of premelting films requires knowledge of the film's equation of state, which is not usually available. Here we calculate the disjoining pressure curve of premelting films, and perform a detailed thermodynamic characterization of premelting behavior on ice. Analysis of the density profiles reveals the signature of weak layering phenomena, from one to two and from two to three water molecular layers. However, disjoining pressure curves, which closely follow expectations from a renormalized mean field liquid state theory, show that there are no layering phase transitions in the thermodynamic sense along the sublimation line. Instead, we find that transitions at mean field level are rounded due to capillary wave fluctuations. We see signatures that true first order layering transitions could arise at low temperatures, for pressures between the metastable line of water/vapor coexistence and the sublimation line. The extrapolation of the disjoining pressure curve above water vapor saturation displays a true first order phase transition from a thin to a thick film consistent with experimental observations. △ Less

Submitted 5 March, 2020; originally announced March 2020.

Comments: 6 pages + 3 figures + supplementary file

Journal ref: Phys. Rev. Lett., 124, 065702 (2020)

A decade of multi-wavelength observations of the TeV blazar 1ES 1215+303: Extreme shift of the synchrotron peak frequency and long-term optical-gamma-ray flux increase

Authors: Janeth Valverde, Deirdre Horan, Denis Bernard, Stephen Fegan, A. U. Abeysekara, A. Archer, W. Benbow, R. Bird, A. Brill, R. Brose, M. Buchovecky, J. H. Buckley, J. L. Christiansen, W. Cui, A. Falcone, Q. Feng, J. P. Finley, L. Fortson, A. Furniss, A. Gent, G. H. Gillanders, C. Giuri, O. Gueta, D. Hanna, T. Hassan , et al. (64 additional authors not shown)

Abstract: Blazars are known for their variability on a wide range of timescales at all wavelengths. Most studies of TeV gamma-ray blazars focus on short timescales, especially during flares. With a decade of observations from the Fermi-LAT and VERITAS, we present an extensive study of the long-term multi-wavelength radio-to-gamma-ray flux-density variability, with the addition of a couple of short-time radi… ▽ More Blazars are known for their variability on a wide range of timescales at all wavelengths. Most studies of TeV gamma-ray blazars focus on short timescales, especially during flares. With a decade of observations from the Fermi-LAT and VERITAS, we present an extensive study of the long-term multi-wavelength radio-to-gamma-ray flux-density variability, with the addition of a couple of short-time radio-structure and optical polarization observations of the blazar 1ES 1215+303 (z=0.130), with a focus on its gamma-ray emission from 100 MeV to 30 TeV. Multiple strong GeV gamma-ray flares, a long-term increase in the gamma-ray and optical flux baseline and a linear correlation between these two bands are observed over the ten-year period. Typical HBL behaviors are identified in the radio morphology and broadband spectrum of the source. Three stationary features in the innermost jet are resolved by VLBA at 43.1, 22.2, and 15.3 GHz. We employ a two-component synchrotron self-Compton model to describe different flux states of the source, including the epoch during which an extreme shift in energy of the synchrotron peak frequency from infrared to soft X-rays is observed. △ Less

Submitted 12 February, 2020; v1 submitted 10 February, 2020; originally announced February 2020.

Comments: Accepted for publication in the Astrophysical Journal

Journal ref: ApJ 891 (2020) 170

Adaptive stochastic continuation with a modified lifting procedure applied to complex systems

Authors: Clemens Willers, Uwe Thiele, Andrew J. Archer, David J. B. Lloyd, Oliver Kamps

Abstract: Many complex systems occurring in the natural or social sciences or economics are frequently described on a microscopic level, e.g., by lattice- or agent-based models. To analyze the states of such systems and their bifurcation structure on the level of macroscopic observables, one has to rely on equation-free methods like stochastic continuation. Here, we investigate how to improve stochastic con… ▽ More Many complex systems occurring in the natural or social sciences or economics are frequently described on a microscopic level, e.g., by lattice- or agent-based models. To analyze the states of such systems and their bifurcation structure on the level of macroscopic observables, one has to rely on equation-free methods like stochastic continuation. Here, we investigate how to improve stochastic continuation techniques by adaptively choosing the parameters of the algorithm. This allows one to obtain bifurcation diagrams quite accurately, especially near bifurcation points. We introduce lifting techniques which generate microscopic states with a naturally grown structure, which can be crucial for a reliable evaluation of macroscopic quantities. We show how to calculate fixed points of fluctuating functions by employing suitable linear fits. This procedure offers a simple measure of the statistical error. We demonstrate these improvements by applying the approach in analyses of (i) the Ising model in two dimensions, (ii) an active Ising model, and (iii) a stochastic Swift-Hohenberg model. We conclude by discussing the abilities and remaining problems of the technique. △ Less

Submitted 6 October, 2020; v1 submitted 5 February, 2020; originally announced February 2020.

Journal ref: Phys. Rev. E 102, 032210 (2020)

Probing the Properties of the Pulsar Wind in the Gamma-Ray Binary HESS J0632+057 with NuSTAR and VERITAS Observations

Authors: A. Archer, W. Benbow, R. Bird, A. Brill, R. Brose, M. Buchovecky, J. L. Christiansen, A. J. Chromey, W. Cui, A. Falcone, Q. Feng, J. P. Finley, L. Fortson, A. Furniss, A. Gent, G. H. Gillanders, C. Giuri, O. Gueta, D. Hanna, T. Hassan, O. Hervet, J. Holder, G. Hughes, T. B. Humensky, P. Kaaret , et al. (38 additional authors not shown)

Abstract: HESS J0632+057 is a gamma-ray binary composed of a compact object orbiting a Be star with a period of about $315$ days. Extensive X-ray and TeV gamma-ray observations have revealed a peculiar light curve containing two peaks, separated by a dip. We present the results of simultaneous observations in hard X-rays with NuSTAR and in TeV gamma-rays with VERITAS, performed in November and December 2017… ▽ More HESS J0632+057 is a gamma-ray binary composed of a compact object orbiting a Be star with a period of about $315$ days. Extensive X-ray and TeV gamma-ray observations have revealed a peculiar light curve containing two peaks, separated by a dip. We present the results of simultaneous observations in hard X-rays with NuSTAR and in TeV gamma-rays with VERITAS, performed in November and December 2017. These observations correspond to the orbital phases $φ\approx0.22$ and $0.3$, where the fluxes are rising towards the first light-curve peak. A significant variation of the spectral index from 1.77$\pm$0.05 to 1.56$\pm$0.05 is observed in the X-ray data. The multi-wavelength spectral energy distributions (SED) derived from the observations are interpreted in terms of a leptonic model, in which the compact object is assumed to be a pulsar and non-thermal radiation is emitted by high-energy electrons accelerated at the shock formed by the collision between the stellar and pulsar wind. The results of the SED fitting show that our data can be consistently described within this scenario, and allow us to estimate the magnetization of the pulsar wind at the location of the shock formation. The constraints on the pulsar-wind magnetization provided by our results are shown to be consistent with those obtained from other systems. △ Less

Submitted 12 December, 2019; v1 submitted 21 November, 2019; originally announced November 2019.

Comments: Accepted for publication in The Astrophysical Journal

Broadband control of water wave energy amplification in chirped arrays

Authors: Aidan J. Archer, Hugh A. Wolgamot, Jana Orszaghova, Luke G. Bennetts, Malte A. Peter, Richard V. Craster

Abstract: Water waves in natural environments are typically broadband, nonlinear and dynamic phenomena. Taking concepts developed for slow light in optics, we address the challenge of designing arrays to control the spatial distribution of wave energy, and amplify target frequencies at specified locations. Experiments on incident waves interacting with a chirped array of eight vertical cylinders demonstrate… ▽ More Water waves in natural environments are typically broadband, nonlinear and dynamic phenomena. Taking concepts developed for slow light in optics, we address the challenge of designing arrays to control the spatial distribution of wave energy, and amplify target frequencies at specified locations. Experiments on incident waves interacting with a chirped array of eight vertical cylinders demonstrate significant amplifications as predicted numerically, and provide motivation for application to energy harvesting. △ Less

Submitted 29 October, 2019; originally announced October 2019.

Comments: 5 pages, 5 figures

Journal ref: Phys. Rev. Fluids 5, 062801 (2020)

Measurement of the extragalactic background light spectral energy distribution with VERITAS

Authors: VERITAS collaboration, A. U. Abeysekara, A. Archer, W. Benbow, R. Bird, A. Brill, R. Brose, M. Buchovecky, J. L. Christiansen, W. Cui, M. K. Daniel, A. Falcone, Q. Feng, M. Fernandez-Alonso, J. P. Finley, L. Fortson, A. Furniss, A. Gent, C. Giuri, O. Gueta, D. Hanna, T. Hassan, O. Hervet, J. Holder, G. Hughes , et al. (37 additional authors not shown)

Abstract: The extragalactic background light (EBL), a diffuse photon field in the optical and infrared range, is a record of radiative processes over the Universe's history. Spectral measurements of blazars at very high energies ($>$100 GeV) enable the reconstruction of the spectral energy distribution (SED) of the EBL, as the blazar spectra are modified by redshift- and energy-dependent interactions of the… ▽ More The extragalactic background light (EBL), a diffuse photon field in the optical and infrared range, is a record of radiative processes over the Universe's history. Spectral measurements of blazars at very high energies ($>$100 GeV) enable the reconstruction of the spectral energy distribution (SED) of the EBL, as the blazar spectra are modified by redshift- and energy-dependent interactions of the gamma-ray photons with the EBL. The spectra of 14 VERITAS-detected blazars are included in a new measurement of the EBL SED that is independent of EBL SED models. The resulting SED covers an EBL wavelength range of 0.56--56 $μ$m, and is in good agreement with lower limits obtained by assuming that the EBL is entirely due to radiation from cataloged galaxies. △ Less

Submitted 1 October, 2019; originally announced October 2019.

Comments: Accepted for publication in The Astrophysical Journal

VERITAS contributions to the 36th International Cosmic Ray Conference

Authors: A. U. Abeysekara, A. Archer, W. Benbow, R. Bird, A. Brill, R. Brose, J. H. Buckley, J. L. Christiansen, A. J. Chromey, M. K. Daniel, S. Das, A. Falcone, Q. Feng, J. P. Finley, L. Fortson, A. Furniss, A. Gent, G. H. Gillanders, C. Giuri, O. Gueta, D. Hanna, T. Hassan, O. Hervet, J. Holder, G. Hughes , et al. (43 additional authors not shown)

Abstract: Compilation of papers presented by the VERITAS Collaboration at the 36th International Cosmic Ray Conference (ICRC), held July 24 through August 1, 2019 in Madison, Wisconsin. Compilation of papers presented by the VERITAS Collaboration at the 36th International Cosmic Ray Conference (ICRC), held July 24 through August 1, 2019 in Madison, Wisconsin. △ Less

Submitted 23 September, 2019; v1 submitted 17 September, 2019; originally announced September 2019.

Comments: html page. 2019 ICRC, Madison, Wisconsin

Small and Practical BERT Models for Sequence Labeling

Authors: Henry Tsai, Jason Riesa, Melvin Johnson, Naveen Arivazhagan, Xin Li, Amelia Archer

Abstract: We propose a practical scheme to train a single multilingual sequence labeling model that yields state of the art results and is small and fast enough to run on a single CPU. Starting from a public multilingual BERT checkpoint, our final model is 6x smaller and 27x faster, and has higher accuracy than a state-of-the-art multilingual baseline. We show that our model especially outperforms on low-re… ▽ More We propose a practical scheme to train a single multilingual sequence labeling model that yields state of the art results and is small and fast enough to run on a single CPU. Starting from a public multilingual BERT checkpoint, our final model is 6x smaller and 27x faster, and has higher accuracy than a state-of-the-art multilingual baseline. We show that our model especially outperforms on low-resource languages, and works on codemixed input text without being explicitly trained on codemixed examples. We showcase the effectiveness of our method by reporting on part-of-speech tagging and morphological prediction on 70 treebanks and 48 languages. △ Less

Submitted 30 August, 2019; originally announced September 2019.

Comments: 11 pages including appendices; accepted to appear at EMNLP-IJCNLP 2019

First order phase transitions and the thermodynamic limit

Authors: Uwe Thiele, Tobias Frohoff-Hülsmann, Sebastian Engelnkemper, Edgar Knobloch, Andrew J. Archer

Abstract: We consider simple mean field continuum models for first order liquid-liquid demixing and solid-liquid phase transitions and show how the Maxwell construction at phase coexistence emerges on going from finite-size closed systems to the thermodynamic limit. The theories considered are the Cahn-Hilliard model of phase separation, which is also a model for the liquid-gas transition, and the phase fie… ▽ More We consider simple mean field continuum models for first order liquid-liquid demixing and solid-liquid phase transitions and show how the Maxwell construction at phase coexistence emerges on going from finite-size closed systems to the thermodynamic limit. The theories considered are the Cahn-Hilliard model of phase separation, which is also a model for the liquid-gas transition, and the phase field crystal model of the solid-liquid transition. Our results show that states comprising the Maxwell line depend strongly on the mean density with spatially localized structures playing a key role in the approach to the thermodynamic limit. △ Less

Submitted 29 January, 2020; v1 submitted 29 August, 2019; originally announced August 2019.

Journal ref: Uwe Thiele et al 2019 New J. Phys. 21 123021

Deriving phase field crystal theory from dynamical density functional theory: consequences of the approximations

Authors: Andrew J. Archer, Daniel J. Ratliff, Alastair M. Rucklidge, Priya Subramanian

Abstract: Phase field crystal (PFC) theory, extensively used for modelling the structure of solids, can be derived from dynamical density functional theory (DDFT) via a sequence of approximations. Standard derivations neglect a term of form $\nabla\cdot[n\nabla L n]$, where $n$ is the scaled density profile and $L$ is a linear operator. We show that this term makes a significant contribution to the stabilit… ▽ More Phase field crystal (PFC) theory, extensively used for modelling the structure of solids, can be derived from dynamical density functional theory (DDFT) via a sequence of approximations. Standard derivations neglect a term of form $\nabla\cdot[n\nabla L n]$, where $n$ is the scaled density profile and $L$ is a linear operator. We show that this term makes a significant contribution to the stability of the crystal, and dropping this term from the theory forces another approximation, that of replacing the logarithmic term from the ideal gas contribution to the free energy with its truncated Taylor expansion, to yield a polynomial in $n$. However, the consequences of doing this are the presence of an additional spinodal in the phase diagram, so the liquid is predicted first to freeze and then to melt again as the density is increased; and other periodic structures are erroneously predicted to be thermodynamic equilibria. A second approximation is to replace $L$ by a gradient expansion. This leads to the possibility of solutions failing to exist above a certain value of the average density. We illustrate these conclusions with a simple model two-dimensional fluid. The consequences of the PFC approximations are that the phase diagram is both qualitatively incorrect, in that it has a stripe phase, and quantitatively incorrect (by orders of magnitude) regarding the properties of the crystal. Thus, although PFC models are successful as phenomenological models of crystallisation, we find it impossible to derive the PFC model as an accurate approximation to DDFT, without introducing spurious artefacts. However, making a simple one-mode approximation for the logarithm of the density distribution is surprisingly accurate, which gives a tantalising hint that accurate PFC-type theories may instead be derived for the field $\log(ρ(x))$, rather than for the density profile itself. △ Less

Submitted 7 August, 2019; originally announced August 2019.

Comments: 28 pages, 9 figures, to appear in Physical Review E

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

Which wavenumbers determine the thermodynamic stability of soft matter quasicrystals?

Authors: D. J. Ratliff, A. J. Archer, P. Subramanian, A. M. Rucklidge

Abstract: For soft matter to form quasicrystals an important ingredient is to have two characteristic lengthscales in the interparticle interactions. To be more precise, for stable quasicrystals, periodic modulations of the local density distribution with two particular wavenumbers should be favored, and the ratio of these wavenumbers should be close to certain special values. So, for simple models, the ans… ▽ More For soft matter to form quasicrystals an important ingredient is to have two characteristic lengthscales in the interparticle interactions. To be more precise, for stable quasicrystals, periodic modulations of the local density distribution with two particular wavenumbers should be favored, and the ratio of these wavenumbers should be close to certain special values. So, for simple models, the answer to the title question is that only these two ingredients are needed. However, for more realistic models, where in principle all wavenumbers can be involved, other wavenumbers are also important, specifically those of the second and higher reciprocal lattice vectors. We identify features in the particle pair interaction potentials which can suppress or encourage density modes with wavenumbers associated with one of the regular crystalline orderings that compete with quasicrystals, enabling either the enhancement or suppression of quasicrystals in a generic class of systems. △ Less

Submitted 12 July, 2019; originally announced July 2019.

Comments: 6 pages, 3 figures

MSC Class: 82D60 (Primary) 35B36; 52C23 (Secondary)

Journal ref: Phys. Rev. Lett. 123, 148004 (2019)

Kinetic Monte Carlo and hydrodynamic modelling of droplet dynamics on surfaces, including evaporation and condensation

Authors: Mounirah Areshi, Dmitri Tseluiko, Andrew J. Archer

Abstract: We present a lattice-gas (generalised Ising) model for liquid droplets on solid surfaces. The time evolution in the model involves two processes: (i) Single-particle moves which are determined by a kinetic Monte Carlo algorithm. These incorporate into the model particle diffusion over the surface and within the droplets and also evaporation and condensation, i.e. the exchange of particles between… ▽ More We present a lattice-gas (generalised Ising) model for liquid droplets on solid surfaces. The time evolution in the model involves two processes: (i) Single-particle moves which are determined by a kinetic Monte Carlo algorithm. These incorporate into the model particle diffusion over the surface and within the droplets and also evaporation and condensation, i.e. the exchange of particles between droplets and the surrounding vapour. (ii) Larger-scale collective moves, modelling advective hydrodynamic fluid motion, determined by considering the dynamics predicted by a thin-film equation. The model enables us to relate how macroscopic quantities such as the contact angle and the surface tension depend on the microscopic interaction parameters between the particles and with the solid surface. We present results for droplets joining, spreading, sliding under gravity, dewetting, the effects of evaporation, the interplay of diffusive and advective dynamics, and how all this behaviour depends on the temperature and other parameters. △ Less

Submitted 2 October, 2019; v1 submitted 19 June, 2019; originally announced June 2019.

Comments: 19 pages, 15 figures

Journal ref: Phys. Rev. Fluids 4, 104006 (2019)

A Search for Pulsed Very High-Energy Gamma Rays from Thirteen Young Pulsars in Archival VERITAS Data

Authors: A. Archer, W. Benbow, R. Bird, R. Brose, M. Buchovecky, J. H. Buckley, A. J. Chromey, W. Cui, A. Falcone, Q. Feng, J. P. Finley, L. Fortson, A. Furniss, A. Gent, O. Gueta, D. Hanna, T. Hassan, O. Hervet, J. Holder, G. Hughes, T. B. Humensky, C. A. Johnson, P. Kaaret, P. Kar, N. Kelley-Hoskins , et al. (36 additional authors not shown)

Abstract: We conduct a search for periodic emission in the very high-energy gamma-ray band (VHE; E > 100 GeV) from a total of 13 pulsars in an archival VERITAS data set with a total exposure of over 450 hours. The set of pulsars includes many of the brightest young gamma-ray pulsars visible in the Northern Hemisphere. The data analysis resulted in non-detections of pulsed VHE gamma rays from each pulsar. Up… ▽ More We conduct a search for periodic emission in the very high-energy gamma-ray band (VHE; E > 100 GeV) from a total of 13 pulsars in an archival VERITAS data set with a total exposure of over 450 hours. The set of pulsars includes many of the brightest young gamma-ray pulsars visible in the Northern Hemisphere. The data analysis resulted in non-detections of pulsed VHE gamma rays from each pulsar. Upper limits on a potential VHE gamma-ray flux are derived at the 95% confidence level above three energy thresholds using two methods. These are the first such searches for pulsed VHE emission from each of the pulsars, and the obtained limits constrain a possible flux component manifesting at VHEs as is seen for the Crab pulsar. △ Less

Submitted 19 April, 2019; originally announced April 2019.

Comments: 32 pages, 20 figures, accepted for publication in ApJ

Binding potentials for vapour nanobubbles on surfaces using density functional theory

Authors: Hanyu Yin, David N. Sibley, Andrew J. Archer

Abstract: We calculate density profiles of a simple model fluid in contact with a planar surface using density functional theory (DFT), in particular for the case where there is a vapour layer intruding between the wall and the bulk liquid. We apply the method of Hughes et al. [J. Chem. Phys. 142, 074702 (2015)] to calculate the density profiles for varying (specified) amounts of the vapour adsorbed at the… ▽ More We calculate density profiles of a simple model fluid in contact with a planar surface using density functional theory (DFT), in particular for the case where there is a vapour layer intruding between the wall and the bulk liquid. We apply the method of Hughes et al. [J. Chem. Phys. 142, 074702 (2015)] to calculate the density profiles for varying (specified) amounts of the vapour adsorbed at the wall. This is equivalent to varying the thickness $h$ of the vapour at the surface. From the resulting sequence of density profiles we calculate the thermodynamic grand potential as $h$ is varied and thereby determine the binding potential as a function of $h$. The binding potential obtained via this coarse-graining approach allows us to determine the disjoining pressure in the film and also to predict the shape of vapour nano-bubbles on the surface. Our microscopic DFT based approach captures information from length scales much smaller than some commonly used models in continuum mechanics. △ Less

Submitted 13 April, 2019; originally announced April 2019.

Comments: 15 pages, 15 figures