Showing 1–50 of 184 results for author: Elliott, J

Dimension reduction in quantum sampling of stochastic processes

Authors: Chengran Yang, Marta Florido-Llin`as, Mile Gu, Thomas J. Elliott

Abstract: Quantum technologies offer a promising route to the efficient sampling and analysis of stochastic processes, with potential applications across the sciences. Such quantum advantages rely on the preparation of a quantum sample state of the stochastic process, which requires a memory system to propagate correlations between the past and future of the process. Here, we introduce a method of lossy qua… ▽ More Quantum technologies offer a promising route to the efficient sampling and analysis of stochastic processes, with potential applications across the sciences. Such quantum advantages rely on the preparation of a quantum sample state of the stochastic process, which requires a memory system to propagate correlations between the past and future of the process. Here, we introduce a method of lossy quantum dimension reduction that allows this memory to be compressed, not just beyond classical limits, but also beyond current state-of-the-art quantum stochastic sampling approaches. We investigate the trade-off between the saving in memory resources from this compression, and the distortion it introduces. We show that our approach can be highly effective in low distortion compression of both Markovian and strongly non-Markovian processes alike. We further discuss the application of our results to quantum stochastic modelling more broadly. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Libfork: portable continuation-stealing with stackless coroutines

Authors: Conor John Williams, James Elliott

Abstract: Fully-strict fork-join parallelism is a powerful model for shared-memory programming due to its optimal time scaling and strong bounds on memory scaling. The latter is rarely achieved due to the difficulty of implementing continuation stealing in traditional High Performance Computing (HPC) languages -- where it is often impossible without modifying the compiler or resorting to non-portable techni… ▽ More Fully-strict fork-join parallelism is a powerful model for shared-memory programming due to its optimal time scaling and strong bounds on memory scaling. The latter is rarely achieved due to the difficulty of implementing continuation stealing in traditional High Performance Computing (HPC) languages -- where it is often impossible without modifying the compiler or resorting to non-portable techniques. We demonstrate how stackless coroutines (a new feature in C++20) can enable fully-portable continuation stealing and present libfork a lock-free fine-grained parallelism library, combining coroutines with user-space, geometric segmented-stacks. We show our approach is able to achieve optimal time/memory scaling, both theoretically and empirically, across a variety of benchmarks. Compared to openMP (libomp), libfork is on average 7.2x faster and consumes 10x less memory. Similarly, compared to Intel's TBB, libfork is on average 2.7x faster and consumes 6.2x less memory. Additionally, we introduce non-uniform memory access (NUMA) optimizations for schedulers that demonstrate performance matching busy-waiting schedulers. △ Less

Submitted 28 February, 2024; originally announced February 2024.

Embedding memory-efficient stochastic simulators as quantum trajectories

Authors: Thomas J. Elliott, Mile Gu

Abstract: By exploiting the complexity intrinsic to quantum dynamics, quantum technologies promise a whole host of computational advantages. One such advantage lies in the field of stochastic modelling, where it has been shown that quantum stochastic simulators can operate with a lower memory overhead than their best classical counterparts. This advantage is particularly pronounced for continuous-time stoch… ▽ More By exploiting the complexity intrinsic to quantum dynamics, quantum technologies promise a whole host of computational advantages. One such advantage lies in the field of stochastic modelling, where it has been shown that quantum stochastic simulators can operate with a lower memory overhead than their best classical counterparts. This advantage is particularly pronounced for continuous-time stochastic processes; however, the corresponding quantum stochastic simulators heretofore prescribed operate only on a quasi-continuous-time basis, and suffer an ever-increasing circuit complexity with increasing temporal resolution. Here, by establishing a correspondence with quantum trajectories -- a method for modelling open quantum systems -- we show how truly continuous-time quantum stochastic simulators can be embedded in such open quantum systems, bridging this gap and obviating previous constraints. We further show how such an embedding can be made for discrete-time stochastic processes, which manifest as jump-only trajectories, and discuss how viewing the correspondence in the reverse direction provides new means of studying structural complexity in quantum systems themselves. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: 16 pages, 4 figures

Journal ref: Phys. Rev. A 109, 022434 (2024)

Challenges of Securing Massively Multiplayer Online Games

Authors: Kolten Sinclair, Steven Womack, Jacob Elliott, Benjamin Stafford, Sundar Krishnan

Abstract: When it comes to security in the modern world, things have improved a lot since the early 2000s. Hypertext Transfer Protocol Secure (HTTPS) and Transport Layer Security (TLS) have made the transfer of our data across the internet much safer than years prior, and the advent of VPNs and private browsing have only compounded that. However, the gaming industry has been notoriously behind the curve whe… ▽ More When it comes to security in the modern world, things have improved a lot since the early 2000s. Hypertext Transfer Protocol Secure (HTTPS) and Transport Layer Security (TLS) have made the transfer of our data across the internet much safer than years prior, and the advent of VPNs and private browsing have only compounded that. However, the gaming industry has been notoriously behind the curve when it comes to security, most notably with Massively Multiplayer Online (MMO) games, which due to the intrinsic nature of their architecture, have an astounding amount of ground to cover. In this paper, the authors discuss the challenges that MMO developers face when trying to design a secure game, as well as some more modern approaches to security that will help improve the industry moving forward. The authors also highlight a few real-life examples of exploits and breaches that have happened and look at how they were mitigated. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Possible Applications of Dissolution Dynamic Nuclear Polarization in Conjunction with Zero- to Ultralow-Field Nuclear Magnetic Resonance

Authors: Danila A. Barskiy, John W. Blanchard, Dmitry Budker, Quentin Stern, James Eills, Stuart J. Elliott, Roman Picazo-Frutos, Antoine Garcon, Sami Jannin, Igor V. Koptyug

Abstract: The combination of a powerful and broadly applicable nuclear hyperpolarization technique with emerging (near-)zero-field modalities offer novel opportunities in a broad range of nuclear magnetic resonance spectroscopy and imaging applications, including biomedical diagnostics, monitoring catalytic reactions within metal reactors and many others. These are discussed along with a roadmap for future… ▽ More The combination of a powerful and broadly applicable nuclear hyperpolarization technique with emerging (near-)zero-field modalities offer novel opportunities in a broad range of nuclear magnetic resonance spectroscopy and imaging applications, including biomedical diagnostics, monitoring catalytic reactions within metal reactors and many others. These are discussed along with a roadmap for future developments. △ Less

Submitted 13 July, 2023; originally announced July 2023.

Comments: 12 pages, 5 figures

Robust internal model control approach for position control of systems with sandwiched backlash

Authors: Yoav Vered, Stephen J. Elliott

Abstract: This paper describes the design of a robust controller for position control in systems with sandwiched backlash. The backlash, which is nonsmooth and nonlinear, is inevitable in the operation of many systems, but it can have destructive effects on the stability and performance of feedback systems. In this work, a robust controller is designed using a modified linear internal model control framewor… ▽ More This paper describes the design of a robust controller for position control in systems with sandwiched backlash. The backlash, which is nonsmooth and nonlinear, is inevitable in the operation of many systems, but it can have destructive effects on the stability and performance of feedback systems. In this work, a robust controller is designed using a modified linear internal model control framework. Different controller architectures are considered and compared based on an experimental case study. The experimental testbased is composed of a three-platform structure driven by a stepper motor. The backlash is introduced into the system in a non-destructive and controllable manner by closing an internal nonlinear feedback loop around the stepper motor. The robustness of the designed controller to a large amount of backlash is verified experimentally, and while the stability is maintained, some residual vibrations are observed. The effects on the residual vibration levels of including nonlinear elements in the controller and changing the controller's settling time are also examined experimentally. The robustness to changes and mismodelling of the linear system, with and without the backlash, is described, as is the tracking of a smooth sinusoidal command signal with a growing amount of backlash. Based on the case study, it is concluded that combining the linear internal model control design method with a small dead zone results in a highly robust controller both with respect to the backlash and to changes in the linear system, which ensure stability and good performance. The required robustness is achieved by tuning the controller's settling time and the dead zone width parameters. △ Less

Submitted 12 July, 2023; originally announced July 2023.

Comments: Preprint

Causal classification of spatiotemporal quantum correlations

Authors: Minjeong Song, Varun Narasimhachar, Bartosz Regula, Thomas J. Elliott, Mile Gu

Abstract: From correlations in measurement outcomes alone, can two otherwise isolated parties establish whether such correlations are atemporal? That is, can they rule out that they have been given the same system at two different times? Classical statistics says no, yet quantum theory disagrees. Here, we introduce the necessary and sufficient conditions by which such quantum correlations can be identified… ▽ More From correlations in measurement outcomes alone, can two otherwise isolated parties establish whether such correlations are atemporal? That is, can they rule out that they have been given the same system at two different times? Classical statistics says no, yet quantum theory disagrees. Here, we introduce the necessary and sufficient conditions by which such quantum correlations can be identified as atemporal. We demonstrate the asymmetry of atemporality under time reversal, and reveal it to be a measure of spatial quantum correlation distinct from entanglement. Our results indicate that certain quantum correlations possess an intrinsic arrow of time, and enable classification of general quantum correlations across space-time based on their (in)compatibility with various underlying causal structures. △ Less

Submitted 15 June, 2023; originally announced June 2023.

Active Noise Control in The New Century: The Role and Prospect of Signal Processing

Authors: Dongyuan Shi, Bhan Lam, Woon-Seng Gan, Jordan Cheer, Stephen J. Elliott

Abstract: Since Paul Leug's 1933 patent application for a system for the active control of sound, the field of active noise control (ANC) has not flourished until the advent of digital signal processors forty years ago. Early theoretical advancements in digital signal processing and processors laid the groundwork for the phenomenal growth of the field, particularly over the past quarter-century. The widespr… ▽ More Since Paul Leug's 1933 patent application for a system for the active control of sound, the field of active noise control (ANC) has not flourished until the advent of digital signal processors forty years ago. Early theoretical advancements in digital signal processing and processors laid the groundwork for the phenomenal growth of the field, particularly over the past quarter-century. The widespread commercial success of ANC in aircraft cabins, automobile cabins, and headsets demonstrates the immeasurable public health and economic benefits of ANC. This article continues where Elliott and Nelson's 1993 Signal Processing Magazine article and Elliott's 1997 50th anniversary commentary on ANC left off, tracing the technical developments and applications in ANC spurred by the seminal texts of Nelson and Elliott (1991), Kuo and Morgan (1996), Hansen and Snyder (1996), and Elliott (2001) since the turn of the century. This article focuses on technical developments pertaining to real-world implementations, such as improving algorithmic convergence, reducing system latency, and extending control to non-stationary and/or broadband noise, as well as the commercial transition challenges from analog to digital ANC systems. Finally, open issues and the future of ANC in the era of artificial intelligence are discussed. △ Less

Submitted 6 July, 2023; v1 submitted 2 June, 2023; originally announced June 2023.

Comments: Submitted to inter.noise 2023, Chiba, Japan

Adversarial Machine Learning and Cybersecurity: Risks, Challenges, and Legal Implications

Authors: Micah Musser, Andrew Lohn, James X. Dempsey, Jonathan Spring, Ram Shankar Siva Kumar, Brenda Leong, Christina Liaghati, Cindy Martinez, Crystal D. Grant, Daniel Rohrer, Heather Frase, Jonathan Elliott, John Bansemer, Mikel Rodriguez, Mitt Regan, Rumman Chowdhury, Stefan Hermanek

Abstract: In July 2022, the Center for Security and Emerging Technology (CSET) at Georgetown University and the Program on Geopolitics, Technology, and Governance at the Stanford Cyber Policy Center convened a workshop of experts to examine the relationship between vulnerabilities in artificial intelligence systems and more traditional types of software vulnerabilities. Topics discussed included the extent… ▽ More In July 2022, the Center for Security and Emerging Technology (CSET) at Georgetown University and the Program on Geopolitics, Technology, and Governance at the Stanford Cyber Policy Center convened a workshop of experts to examine the relationship between vulnerabilities in artificial intelligence systems and more traditional types of software vulnerabilities. Topics discussed included the extent to which AI vulnerabilities can be handled under standard cybersecurity processes, the barriers currently preventing the accurate sharing of information about AI vulnerabilities, legal issues associated with adversarial attacks on AI systems, and potential areas where government support could improve AI vulnerability management and mitigation. This report is meant to accomplish two things. First, it provides a high-level discussion of AI vulnerabilities, including the ways in which they are disanalogous to other types of vulnerabilities, and the current state of affairs regarding information sharing and legal oversight of AI vulnerabilities. Second, it attempts to articulate broad recommendations as endorsed by the majority of participants at the workshop. △ Less

Submitted 23 May, 2023; originally announced May 2023.

Inverse projection of axisymmetric orientation distributions

Authors: Philipp A. Kloza, James A. Elliott

Abstract: We show that the projection of an axisymmetric three-dimensional orientation distribution to two dimensions can be cast into an Abel transform. Based on this correspondence, we derive an exact integral inverse, which allows for the quantification of three-dimensional uniaxial alignment of rodlike units from two-dimensional sliced images, thus providing an alternative to X-ray or tomographic analys… ▽ More We show that the projection of an axisymmetric three-dimensional orientation distribution to two dimensions can be cast into an Abel transform. Based on this correspondence, we derive an exact integral inverse, which allows for the quantification of three-dimensional uniaxial alignment of rodlike units from two-dimensional sliced images, thus providing an alternative to X-ray or tomographic analysis. A matrix representation of the projection and its inverse is derived, providing a direct relationship between two- and three-dimensional order parameters for both polar and non-polar systems. △ Less

Submitted 10 May, 2023; originally announced May 2023.

Comments: 7 pages, 3 figures, 1 table, submitted to Europhysics Letters

The James Webb Space Telescope Mission

Authors: Jonathan P. Gardner, John C. Mather, Randy Abbott, James S. Abell, Mark Abernathy, Faith E. Abney, John G. Abraham, Roberto Abraham, Yasin M. Abul-Huda, Scott Acton, Cynthia K. Adams, Evan Adams, David S. Adler, Maarten Adriaensen, Jonathan Albert Aguilar, Mansoor Ahmed, Nasif S. Ahmed, Tanjira Ahmed, Rüdeger Albat, Loïc Albert, Stacey Alberts, David Aldridge, Mary Marsha Allen, Shaune S. Allen, Martin Altenburg , et al. (983 additional authors not shown)

Abstract: Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astrono… ▽ More Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astronomers will celebrate their accomplishments for the life of the mission, potentially as long as 20 years, and beyond. This report and the scientific discoveries that follow are extended thank-you notes to the 20,000 team members. The telescope is working perfectly, with much better image quality than expected. In this and accompanying papers, we give a brief history, describe the observatory, outline its objectives and current observing program, and discuss the inventions and people who made it possible. We cite detailed reports on the design and the measured performance on orbit. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: Accepted by PASP for the special issue on The James Webb Space Telescope Overview, 29 pages, 4 figures

ShadowNav: Crater-Based Localization for Nighttime and Permanently Shadowed Region Lunar Navigation

Authors: Abhishek Cauligi, R. Michael Swan, Masahiro Ono, Shreyansh Daftry, John Elliott, Larry Matthies, Deegan Atha

Abstract: There has been an increase in interest in missions that drive significantly longer distances per day than what has currently been performed. Further, some of these proposed missions require autonomous driving and absolute localization in darkness. For example, the Endurance A mission proposes to drive 1200km of its total traverse at night. The lack of natural light available during such missions l… ▽ More There has been an increase in interest in missions that drive significantly longer distances per day than what has currently been performed. Further, some of these proposed missions require autonomous driving and absolute localization in darkness. For example, the Endurance A mission proposes to drive 1200km of its total traverse at night. The lack of natural light available during such missions limits what can be used as visual landmarks and the range at which landmarks can be observed. In order for planetary rovers to traverse long ranges, onboard absolute localization is critical to the ability of the rover to maintain its planned trajectory and avoid known hazardous regions. Currently, to accomplish absolute localization, a ground in the loop (GITL) operation is performed wherein a human operator matches local maps or images from onboard with orbital images and maps. This GITL operation limits the distance that can be driven in a day to a few hundred meters, which is the distance that the rover can maintain acceptable localization error via relative methods. Previous work has shown that using craters as landmarks is a promising approach for performing absolute localization on the moon during the day. In this work we present a method of absolute localization that utilizes craters as landmarks and matches detected crater edges on the surface with known craters in orbital maps. We focus on a localization method based on a perception system which has an external illuminator and a stereo camera. We evaluate (1) both monocular and stereo based surface crater edge detection techniques, (2) methods of scoring the crater edge matches for optimal localization, and (3) localization performance on simulated Lunar surface imagery at night. We demonstrate that this technique shows promise for maintaining absolute localization error of less than 10m required for most planetary rover missions. △ Less

Submitted 11 January, 2023; originally announced January 2023.

Comments: IEEE Aerospace Conference 2023

Constant Chemical Potential-Quantum Mechanical-Molecular Dynamics simulations of the Graphene-electrolyte double layer

Authors: Nicodemo Di Pasquale, Aaron R. Finney, Joshua Elliott, Paola Carbone, Matteo Salvalaglio

Abstract: We present the coupling of two frameworks -- the pseudo-open boundary simulation method known as constant potential Molecular Dynamics simulations (C$μ$MD), combined with QMMD calculations -- to describe the properties of graphene electrodes in contact with electrolytes. The resulting C$μ$QMMD model was then applied to three ionic solutions (LiCl, NaCl and KCl in water) at bulk solution concentrat… ▽ More We present the coupling of two frameworks -- the pseudo-open boundary simulation method known as constant potential Molecular Dynamics simulations (C$μ$MD), combined with QMMD calculations -- to describe the properties of graphene electrodes in contact with electrolytes. The resulting C$μ$QMMD model was then applied to three ionic solutions (LiCl, NaCl and KCl in water) at bulk solution concentrations ranging from 0.5 M up to 6 M in contact with a charged graphene electrode. The new approach we are describing here provides a simulation protocol to control the concentration of the electrolyte solutions while including the effects of a fully polarizable electrode surface. Thanks to this coupling, we are able to accurately model both the electrode and solution side of the double layer and provide a thorough analysis of the properties of electrolytes at charged interfaces, such as the screening ability of the electrolyte and the electrostatic potential profile. We also report the calculation of the integral electrochemical double layer capacitance in the whole range of concentrations analysed for each ionic species, while the QM simulations provide access to the differential and integral quantum capacitance. We highlight how subtle features, such as the adsorption of potassium at the interface or the tendency of the ions to form clusters, emerge from our simulations, contribute to explaining the ability of graphene to store charge and suggest implications for desalination. △ Less

Submitted 9 December, 2022; v1 submitted 7 December, 2022; originally announced December 2022.

Comments: 28 pages, 10 figures

Implementing quantum dimensionality reduction for non-Markovian stochastic simulation

Authors: Kang-Da Wu, Chengran Yang, Ren-Dong He, Mile Gu, Guo-Yong Xiang, Chuan-Feng Li, Guang-Can Guo, Thomas J. Elliott

Abstract: Complex systems are embedded in our everyday experience. Stochastic modelling enables us to understand and predict the behaviour of such systems, cementing its utility across the quantitative sciences. Accurate models of highly non-Markovian processes -- where the future behaviour depends on events that happened far in the past -- must track copious amounts of information about past observations,… ▽ More Complex systems are embedded in our everyday experience. Stochastic modelling enables us to understand and predict the behaviour of such systems, cementing its utility across the quantitative sciences. Accurate models of highly non-Markovian processes -- where the future behaviour depends on events that happened far in the past -- must track copious amounts of information about past observations, requiring high-dimensional memories. Quantum technologies can ameliorate this cost, allowing models of the same processes with lower memory dimension than corresponding classical models. Here we implement such memory-efficient quantum models for a family of non-Markovian processes using a photonic setup. We show that with a single qubit of memory our implemented quantum models can attain higher precision than possible with any classical model of the same memory dimension. This heralds a key step towards applying quantum technologies in complex systems modelling. △ Less

Submitted 18 October, 2023; v1 submitted 26 August, 2022; originally announced August 2022.

Comments: 11+2 pages, 5 figures

Journal ref: Nature Communications 14, 2624 (2023)

The Science Performance of JWST as Characterized in Commissioning

Authors: Jane Rigby, Marshall Perrin, Michael McElwain, Randy Kimble, Scott Friedman, Matt Lallo, René Doyon, Lee Feinberg, Pierre Ferruit, Alistair Glasse, Marcia Rieke, George Rieke, Gillian Wright, Chris Willott, Knicole Colon, Stefanie Milam, Susan Neff, Christopher Stark, Jeff Valenti, Jim Abell, Faith Abney, Yasin Abul-Huda, D. Scott Acton, Evan Adams, David Adler , et al. (601 additional authors not shown)

Abstract: This paper characterizes the actual science performance of the James Webb Space Telescope (JWST), as determined from the six month commissioning period. We summarize the performance of the spacecraft, telescope, science instruments, and ground system, with an emphasis on differences from pre-launch expectations. Commissioning has made clear that JWST is fully capable of achieving the discoveries f… ▽ More This paper characterizes the actual science performance of the James Webb Space Telescope (JWST), as determined from the six month commissioning period. We summarize the performance of the spacecraft, telescope, science instruments, and ground system, with an emphasis on differences from pre-launch expectations. Commissioning has made clear that JWST is fully capable of achieving the discoveries for which it was built. Moreover, almost across the board, the science performance of JWST is better than expected; in most cases, JWST will go deeper faster than expected. The telescope and instrument suite have demonstrated the sensitivity, stability, image quality, and spectral range that are necessary to transform our understanding of the cosmos through observations spanning from near-earth asteroids to the most distant galaxies. △ Less

Submitted 10 April, 2023; v1 submitted 12 July, 2022; originally announced July 2022.

Comments: 5th version as accepted to PASP; 31 pages, 18 figures; https://iopscience.iop.org/article/10.1088/1538-3873/acb293

Journal ref: PASP 135 048001 (2023)

Bit complexity for computing one point in each connected component of a smooth real algebraic set

Authors: Jesse Elliott, Mark Giesbrecht, Eric Schost

Abstract: We analyze the bit complexity of an algorithm for the computation of at least one point in each connected component of a smooth real algebraic set. This work is a continuation of our analysis of the hypersurface case (On the bit complexity of finding points in connected components of a smooth real hypersurface, ISSAC'20). In this paper, we extend the analysis to more general cases. Let… ▽ More We analyze the bit complexity of an algorithm for the computation of at least one point in each connected component of a smooth real algebraic set. This work is a continuation of our analysis of the hypersurface case (On the bit complexity of finding points in connected components of a smooth real hypersurface, ISSAC'20). In this paper, we extend the analysis to more general cases. Let $F=(f_1,..., f_p)$ in $\mathbb{Z}[X_1, ... , X_n]^p$ be a sequence of polynomials with $V = V(F) \subset \mathbb{C}^n$ a smooth and equidimensional variety and $\langle F \rangle \subset \mathbb{C}[X_1, ..., X_n]$ a radical ideal. To compute at least one point in each connected component of $V \cap \mathbb{R}^n$, our starting point is an algorithm by Safey El Din and Schost (Polar varieties and computation of one point in each connected component of a smooth real algebraic set, ISSAC'03). This algorithm uses random changes of variables that are proven to generically ensure certain desirable geometric properties. The cost of the algorithm was given in an algebraic complexity model; here, we analyze the bit complexity and the error probability, and we provide a quantitative analysis of the genericity statements. In particular, we are led to use Lagrange systems to describe polar varieties, as they make it simpler to rely on techniques such as weak transversality and an effective Nullstellensatz. △ Less

Submitted 9 July, 2022; originally announced July 2022.

Comments: Journal of Symbolic Computation, to appear

Wettability of graphite under 2D confinement

Authors: Zixuan Wei, Mara Chiricotto, Joshua D. Elliott, Fausto Martelli, Paola Carbone

Abstract: The thermodynamics of solid/liquid interfaces under nanoconfinement has tremendous implications for liquid transport properties. Here using molecular dynamics, we investigate graphite nanoslits and study how the water/graphite interfacial tension changes with the degree of confinement. We found that, for nanochannel heights between 0.7nm and 2.6nm, graphite becomes more hydrophobic than in bulk, a… ▽ More The thermodynamics of solid/liquid interfaces under nanoconfinement has tremendous implications for liquid transport properties. Here using molecular dynamics, we investigate graphite nanoslits and study how the water/graphite interfacial tension changes with the degree of confinement. We found that, for nanochannel heights between 0.7nm and 2.6nm, graphite becomes more hydrophobic than in bulk, and that the value of the surface tension oscillates before eventually converging towards a constant value for larger slits. The value of the surface tension is correlated with the slip length of the fluid and explained in terms of the effective and interfacial density, hydration pressure and friction coefficient. The study clearly indicates that there is a critical channel height of 0.9nm (achievable experimentally1) at which the surface tension reaches its highest value, but the water diffusion across the channel is at its minimum. The structural analysis shows that for this pore size a transition between a 2D and 3D hydrogen bond network is accompanied by an abrupt increase in conformational entropy. Our results show that the wettability of solid surfaces can change under nanoconfinement and the data can be used to interpret the experimental permeability data. △ Less

Submitted 16 May, 2022; originally announced May 2022.

Hidden Quantum Memory: Is Memory There When Somebody Looks?

Authors: Philip Taranto, Thomas J. Elliott, Simon Milz

Abstract: In classical physics, memoryless dynamics and Markovian statistics are one and the same. This is not true for quantum dynamics, first and foremost because quantum measurements are invasive. Going beyond measurement invasiveness, here we derive a novel distinction between classical and quantum processes, namely the possibility of hidden quantum memory. While Markovian statistics of classical proces… ▽ More In classical physics, memoryless dynamics and Markovian statistics are one and the same. This is not true for quantum dynamics, first and foremost because quantum measurements are invasive. Going beyond measurement invasiveness, here we derive a novel distinction between classical and quantum processes, namely the possibility of hidden quantum memory. While Markovian statistics of classical processes can always be reproduced by a memoryless dynamical model, our main result shows that this is not true in quantum mechanics: We first provide an example of quantum non-Markovianity whose manifestation depends on whether or not a previous measurement is performed -- an impossible phenomenon for memoryless dynamics; we then strengthen this result by demonstrating statistics that are Markovian independent of how they are probed, but are nonetheless still incompatible with memoryless quantum dynamics. Thus, we establish the existence of Markovian statistics gathered by probing a quantum process that nevertheless fundamentally require memory for their creation. △ Less

Submitted 24 April, 2023; v1 submitted 18 April, 2022; originally announced April 2022.

Comments: 7 + 9 pages, 5 figures

Journal ref: Quantum 7, 991 (2023)

Simultaneously Enhanced Tenacity, Rupture Work, and Thermal Conductivity of Carbon Nanotubes Fibers by Increasing the Effective Tube Contribution

Authors: Xiao Zhang, Michael De Volder, Wenbin Zhou, Liron Issman, Xiaojun Wei, Adarsh Kaniyoor, Jeronimo Terrones Portas, Fiona Smail, Zibo Wang, Yanchun Wang, Huaping Liu, Weiya Zhou, James Elliott, Sishen Xie, Adam Boies

Abstract: Although individual carbon nanotubes (CNTs) are superior as constituents to polymer chains, the mechanical and thermal properties of CNT fibers (CNTFs) remain inferior to commercial synthetic fibers due to the lack of synthesis methods to embed CNTs effectively in superstructures. The application of conventional techniques for mechanical enhancement resulted in a mild improvement of target propert… ▽ More Although individual carbon nanotubes (CNTs) are superior as constituents to polymer chains, the mechanical and thermal properties of CNT fibers (CNTFs) remain inferior to commercial synthetic fibers due to the lack of synthesis methods to embed CNTs effectively in superstructures. The application of conventional techniques for mechanical enhancement resulted in a mild improvement of target properties while achieving parity at best on others. In this work, a Double-Drawing technique is developed to deform continuously grown CNTFs and rearrange the constituent CNTs in both mesoscale and nanoscale morphology. Consequently, the mechanical and thermal properties of the resulting CNTFs can be jointly improved, and simultaneously reach their highest performances with specific strength (tenacity) $\rm\sim3.30\,N\,tex^{-1}$, work of rupture $\rm\sim70\,J\,g^{-1}$, and thermal conductivity $\rm\sim354\,W\,m^{-1}\,K^{-1}$, despite starting from commercial low-crystallinity materials ($I{\rm_G}:I{\rm_D}\sim5$). The processed CNTFs are more versatile than comparable carbon fiber, Zylon, Dyneema, and Kevlar. Furthermore, based on evidence of load transfer efficiency on individual CNTs measured with In-Situ Stretching Raman, we find the main contributors to property enhancements are (1) the increased proportion of load-bearing CNT bundles and (2) the extension of effective length of tubes attached on these bundles. △ Less

Submitted 15 April, 2022; v1 submitted 9 April, 2022; originally announced April 2022.

Comments: 27 pages, 8 figures, research article (main text+supporting information)

Report number: eabq3515

Journal ref: Science Advances 2022

Do specific ion effects influence the physical chemistry of aqueous graphene-based supercapacitors? Perspectives from multiscale QMMD simulations

Authors: Joshua D. Elliott, Mara Chiricotto, Alessandro Troisi, Paola Carbone

Abstract: Whether or not specific ion effects determine the charge storage properties of aqueous graphene and graphite-based supercapacitors remains a highly debated topic. In this work we present a multiscale quantum mechanics classical molecular dynamics investigation of aqueous mono- and divalent salt electrolytes in contact with fully polarizable charged graphene sheets. By computing both the electroche… ▽ More Whether or not specific ion effects determine the charge storage properties of aqueous graphene and graphite-based supercapacitors remains a highly debated topic. In this work we present a multiscale quantum mechanics classical molecular dynamics investigation of aqueous mono- and divalent salt electrolytes in contact with fully polarizable charged graphene sheets. By computing both the electrochemical double layer and quantum capacitance we observe a constant electrode specific capacitance with cationic radii and charge. Counterintuitively, we determine that a switch in the cation adsorption mechanism from inner to outer Helmholtz layers leads to negligible changes to the EDL capacitance, this appears to be due to the robust electronic structure of the graphene electrodes. However, the ability of ions (such as K+) with a relatively low hydration free energy to penetrate the inner Helmholtz plane and adsorb directly on the electrode surface is found to slow their diffusion parallel to the interface. Ions in the outer Helmholtz layer are found to have higher diffusivity at the surface due to their position in ion channels between water layers. Our results show that surface effects such as the surface polarization and the partial dehydration and local structuring of ions on the surface underpin the behaviour of cations at the interface and add a vital new perspective on trends in ion mobilities seen under confinement. △ Less

Submitted 9 November, 2022; v1 submitted 4 March, 2022; originally announced March 2022.

Comments: 23 pages, 5 figures 4 tables

Regime Switching Entropic Risk Measures on Crude Oil Pricing

Authors: Babacar Seck, Robert J. Elliott

Abstract: This paper introduces a new type of risk measures, namely regime switching entropic risk measures, and study their applicability through simulations. The state of the economy is incorporated into the entropic risk formulation by using a Markov chain. Closed formulae of the risk measure are obtained for futures on crude oil derivatives. The applicability of these new types of risk measures is based… ▽ More This paper introduces a new type of risk measures, namely regime switching entropic risk measures, and study their applicability through simulations. The state of the economy is incorporated into the entropic risk formulation by using a Markov chain. Closed formulae of the risk measure are obtained for futures on crude oil derivatives. The applicability of these new types of risk measures is based on the study of the risk aversion parameter and the convenience yield. The numerical results show a term structure and a mean-reverting behavior of the convenience yield. △ Less

Submitted 24 December, 2021; originally announced December 2021.

Comments: 20 pages, 6 figures

MSC Class: 91B24; 91B30; 91B70

Backward stochastic differential equations with regime-switching and sublinear expectations

Authors: Engel John C. Dela Vega, Robert J. Elliott

Abstract: This paper introduces a backward stochastic differential equation driven by both Brownian motion and a Markov chain (BSDEBM). Regime-switching is also incorporated through its driver. The existence and uniqueness of the solution of the BSDEBM are proved. A comparison theorem is also derived. Filtration consistent sublinear expectations are defined and characterized as solutions to the BSDEBM. The… ▽ More This paper introduces a backward stochastic differential equation driven by both Brownian motion and a Markov chain (BSDEBM). Regime-switching is also incorporated through its driver. The existence and uniqueness of the solution of the BSDEBM are proved. A comparison theorem is also derived. Filtration consistent sublinear expectations are defined and characterized as solutions to the BSDEBM. The bid and ask prices are then represented using sublinear expectations. △ Less

Submitted 4 December, 2021; originally announced December 2021.

Comments: 19 pages

MSC Class: 60H10 (Primary) 60J28 (Secondary)

A stochastic control approach to bid-ask price modelling

Authors: Engel John C. Dela Vega, Robert J. Elliott

Abstract: This paper develops a model for the bid and ask prices of a European type asset by formulating a stochastic control problem. The state process is governed by a modified geometric Brownian motion whose drift and diffusion coefficients depend on a Markov chain. A Girsanov theorem for Markov chains is implemented for the change of coefficients, including the diffusion coefficient which cannot be chan… ▽ More This paper develops a model for the bid and ask prices of a European type asset by formulating a stochastic control problem. The state process is governed by a modified geometric Brownian motion whose drift and diffusion coefficients depend on a Markov chain. A Girsanov theorem for Markov chains is implemented for the change of coefficients, including the diffusion coefficient which cannot be changed by the usual Girsanov theorem for Brownian motion. The price of a European type asset is then determined using an Esscher transform and a system of partial differential equations. A dynamic programming principle and a maximum/minimum principle associated with the stochastic control problem are then derived to model bid and ask prices. These prices are not quotes of traders or market makers but represent estimates in our model on which reasonable quantities could be traded. △ Less

Submitted 4 December, 2021; originally announced December 2021.

Comments: 22 pages

MSC Class: 93E20 (Primary) 91G20; 60J27 (Secondary)

Quantum adaptive agents with efficient long-term memories

Authors: Thomas J. Elliott, Mile Gu, Andrew J. P. Garner, Jayne Thompson

Abstract: Central to the success of adaptive systems is their ability to interpret signals from their environment and respond accordingly -- they act as agents interacting with their surroundings. Such agents typically perform better when able to execute increasingly complex strategies. This comes with a cost: the more information the agent must recall from its past experiences, the more memory it will need… ▽ More Central to the success of adaptive systems is their ability to interpret signals from their environment and respond accordingly -- they act as agents interacting with their surroundings. Such agents typically perform better when able to execute increasingly complex strategies. This comes with a cost: the more information the agent must recall from its past experiences, the more memory it will need. Here we investigate the power of agents capable of quantum information processing. We uncover the most general form a quantum agent need adopt to maximise memory compression advantages, and provide a systematic means of encoding their memory states. We show these encodings can exhibit extremely favourable scaling advantages relative to memory-minimal classical agents, particularly when information must be retained about events increasingly far into the past. △ Less

Submitted 11 January, 2022; v1 submitted 24 August, 2021; originally announced August 2021.

Comments: 16 pages, 4 figures

Journal ref: Physical Review X 12, 011007 (2022)

The RSNA-ASNR-MICCAI BraTS 2021 Benchmark on Brain Tumor Segmentation and Radiogenomic Classification

Authors: Ujjwal Baid, Satyam Ghodasara, Suyash Mohan, Michel Bilello, Evan Calabrese, Errol Colak, Keyvan Farahani, Jayashree Kalpathy-Cramer, Felipe C. Kitamura, Sarthak Pati, Luciano M. Prevedello, Jeffrey D. Rudie, Chiharu Sako, Russell T. Shinohara, Timothy Bergquist, Rong Chai, James Eddy, Julia Elliott, Walter Reade, Thomas Schaffter, Thomas Yu, Jiaxin Zheng, Ahmed W. Moawad, Luiz Otavio Coelho, Olivia McDonnell , et al. (78 additional authors not shown)

Abstract: The BraTS 2021 challenge celebrates its 10th anniversary and is jointly organized by the Radiological Society of North America (RSNA), the American Society of Neuroradiology (ASNR), and the Medical Image Computing and Computer Assisted Interventions (MICCAI) society. Since its inception, BraTS has been focusing on being a common benchmarking venue for brain glioma segmentation algorithms, with wel… ▽ More The BraTS 2021 challenge celebrates its 10th anniversary and is jointly organized by the Radiological Society of North America (RSNA), the American Society of Neuroradiology (ASNR), and the Medical Image Computing and Computer Assisted Interventions (MICCAI) society. Since its inception, BraTS has been focusing on being a common benchmarking venue for brain glioma segmentation algorithms, with well-curated multi-institutional multi-parametric magnetic resonance imaging (mpMRI) data. Gliomas are the most common primary malignancies of the central nervous system, with varying degrees of aggressiveness and prognosis. The RSNA-ASNR-MICCAI BraTS 2021 challenge targets the evaluation of computational algorithms assessing the same tumor compartmentalization, as well as the underlying tumor's molecular characterization, in pre-operative baseline mpMRI data from 2,040 patients. Specifically, the two tasks that BraTS 2021 focuses on are: a) the segmentation of the histologically distinct brain tumor sub-regions, and b) the classification of the tumor's O[6]-methylguanine-DNA methyltransferase (MGMT) promoter methylation status. The performance evaluation of all participating algorithms in BraTS 2021 will be conducted through the Sage Bionetworks Synapse platform (Task 1) and Kaggle (Task 2), concluding in distributing to the top ranked participants monetary awards of $60,000 collectively. △ Less

Submitted 12 September, 2021; v1 submitted 5 July, 2021; originally announced July 2021.

Comments: 19 pages, 2 figures, 1 table

Quantum coarse-graining for extreme dimension reduction in modelling stochastic temporal dynamics

Authors: Thomas J. Elliott

Abstract: Stochastic modelling of complex systems plays an essential, yet often computationally intensive role across the quantitative sciences. Recent advances in quantum information processing have elucidated the potential for quantum simulators to exhibit memory advantages for such tasks. Heretofore, the focus has been on lossless memory compression, wherein the advantage is typically in terms of lesseni… ▽ More Stochastic modelling of complex systems plays an essential, yet often computationally intensive role across the quantitative sciences. Recent advances in quantum information processing have elucidated the potential for quantum simulators to exhibit memory advantages for such tasks. Heretofore, the focus has been on lossless memory compression, wherein the advantage is typically in terms of lessening the amount of information tracked by the model, while -- arguably more practical -- reductions in memory dimension are not always possible. Here we address the case of lossy compression for quantum stochastic modelling of continuous-time processes, introducing a method for coarse-graining in quantum state space that drastically reduces the requisite memory dimension for modelling temporal dynamics whilst retaining near-exact statistics. In contrast to classical coarse-graining, this compression is not based on sacrificing temporal resolution, and brings memory-efficient, high-fidelity stochastic modelling within reach of present quantum technologies. △ Less

Submitted 21 June, 2021; v1 submitted 14 May, 2021; originally announced May 2021.

Comments: 14 pages, 9 figures

Journal ref: PRX Quantum 2, 020342 (2021)

Memory compression and thermal efficiency of quantum implementations of non-deterministic hidden Markov models

Authors: Thomas J. Elliott

Abstract: Stochastic modelling is an essential component of the quantitative sciences, with hidden Markov models (HMMs) often playing a central role. Concurrently, the rise of quantum technologies promises a host of advantages in computational problems, typically in terms of the scaling of requisite resources such as time and memory. HMMs are no exception to this, with recent results highlighting quantum im… ▽ More Stochastic modelling is an essential component of the quantitative sciences, with hidden Markov models (HMMs) often playing a central role. Concurrently, the rise of quantum technologies promises a host of advantages in computational problems, typically in terms of the scaling of requisite resources such as time and memory. HMMs are no exception to this, with recent results highlighting quantum implementations of deterministic HMMs exhibiting superior memory and thermal efficiency relative to their classical counterparts. In many contexts however, non-deterministic HMMs are viable alternatives; compared to them the advantages of current quantum implementations do not always hold. Here, we provide a systematic prescription for constructing quantum implementations of non-deterministic HMMs that re-establish the quantum advantages against this broader class. Crucially, we show that whenever the classical implementation suffers from thermal dissipation due to its need to process information in a time-local manner, our quantum implementations will both mitigate some of this dissipation, and achieve an advantage in memory compression. △ Less

Submitted 21 June, 2021; v1 submitted 13 May, 2021; originally announced May 2021.

Comments: 10 pages, 1 figure, 1 table

Journal ref: Physical Review A, 103, 052615 (2021)

Dynamically polarisable force-fields for surface simulations via multi-output classification Neural Networks

Authors: Nicodemo Di Pasquale, Joshua D. Elliott, Panagiotis Hadjidoukas, Paola Carbone

Abstract: We present a general procedure to introduce electronic polarization into classical Molecular Dynamics (MD) force-fields using a Neural Network (NN) model. We apply this framework to the simulation of a solid-liquid interface where the polarization of the surface is essential to correctly capture the main features of the system. By introducing a multi-input, multi-output NN and treating the surface… ▽ More We present a general procedure to introduce electronic polarization into classical Molecular Dynamics (MD) force-fields using a Neural Network (NN) model. We apply this framework to the simulation of a solid-liquid interface where the polarization of the surface is essential to correctly capture the main features of the system. By introducing a multi-input, multi-output NN and treating the surface polarization as a discrete classification problem, for which NNs are known to excel, we are able to obtain very good accuracy in terms of quality of predictions. Through the definition of a custom loss function we are able to impose a physically motivated constraint within the NN itself making this model extremely versatile, especially in the modelling of different surface charge states. The NN is validated considering the redistribution of electronic charge density within a graphene based electrode in contact with aqueous electrolyte solution, a system highly relevant to the development of next generation low-cost supercapacitors. We compare the performances of our NN/MD model against Quantum Mechanics/Molecular dynamics simulations where we obtain a most satisfactorily agreement. △ Less

Submitted 30 March, 2021; originally announced March 2021.

Comments: 6 figures

Quantum-inspired identification of complex cellular automata

Authors: Matthew Ho, Andri Pradana, Thomas J. Elliott, Lock Yue Chew, Mile Gu

Abstract: Elementary cellular automata (ECA) present iconic examples of complex systems. Though described only by one-dimensional strings of binary cells evolving according to nearest-neighbour update rules, certain ECA rules manifest complex dynamics capable of universal computation. Yet, the classification of precisely which rules exhibit complex behaviour remains a significant challenge. Here we approach… ▽ More Elementary cellular automata (ECA) present iconic examples of complex systems. Though described only by one-dimensional strings of binary cells evolving according to nearest-neighbour update rules, certain ECA rules manifest complex dynamics capable of universal computation. Yet, the classification of precisely which rules exhibit complex behaviour remains a significant challenge. Here we approach this question using tools from quantum stochastic modelling, where quantum statistical memory -- the memory required to model a stochastic process using a class of quantum machines -- can be used to quantify the structure of a stochastic process. By viewing ECA rules as transformations of stochastic patterns, we ask: Does an ECA generate structure as quantified by the quantum statistical memory, and if so, how quickly? We illustrate how the growth of this measure over time correctly distinguishes simple ECA from complex counterparts. Moreover, it provides a more refined means for quantitatively identifying complex ECAs -- providing a spectrum on which we can rank the complexity of ECA by the rate in which they generate structure. △ Less

Submitted 20 March, 2024; v1 submitted 25 March, 2021; originally announced March 2021.

Comments: 22 pages, 9 figures

Journal ref: Eur. Phys. J. Plus 138 (6) 540 (2023)

Efficient exploration and calibration of a semi-analytical model of galaxy formation with deep learning

Authors: Edward J. Elliott, Carlton M. Baugh, Cedric G. Lacey

Abstract: We implement a sample-efficient method for rapid and accurate emulation of semi-analytical galaxy formation models over a wide range of model outputs. We use ensembled deep learning algorithms to produce a fast emulator of an updated version of the GALFORM model from a small number of training examples. We use the emulator to explore the model's parameter space, and apply sensitivity analysis tech… ▽ More We implement a sample-efficient method for rapid and accurate emulation of semi-analytical galaxy formation models over a wide range of model outputs. We use ensembled deep learning algorithms to produce a fast emulator of an updated version of the GALFORM model from a small number of training examples. We use the emulator to explore the model's parameter space, and apply sensitivity analysis techniques to better understand the relative importance of the model parameters. We uncover key tensions between observational datasets by applying a heuristic weighting scheme in a Markov chain Monte Carlo framework and exploring the effects of requiring improved fits to certain datasets relative to others. Furthermore, we demonstrate that this method can be used to successfully calibrate the model parameters to a comprehensive list of observational constraints. In doing so, we re-discover previous GALFORM fits in an automatic and transparent way, and discover an improved fit by applying a heavier weighting to the fit to the metallicities of early-type galaxies. The deep learning emulator requires a fraction of the model evaluations needed in similar emulation approaches, achieving an out-of-sample mean absolute error at the knee of the K-band luminosity function of 0.06 dex with less than 1000 model evaluations. We demonstrate that this is an extremely efficient, inexpensive and transparent way to explore multi-dimensional parameter spaces, and can be applied more widely beyond semi-analytical galaxy formation models. △ Less

Submitted 1 March, 2021; originally announced March 2021.

Comments: 21 pages, 15 figures, submitted to MNRAS

Deterministic single-atom source of quasi-superradiant $N$-photon pulses

Authors: Caspar Groiseau, Alexander E. J. Elliott, Stuart J. Masson, Scott Parkins

Abstract: We propose a single-atom, cavity quantum electrodynamics system, compatible with recently demonstrated, fiber-integrated micro- and nano-cavity setups, for the on-demand production of optical number-state, $0N$-state, and binomial-code-state pulses. The scheme makes use of Raman transitions within an entire atomic ground-state hyperfine level and operates with laser and cavity fields detuned from… ▽ More We propose a single-atom, cavity quantum electrodynamics system, compatible with recently demonstrated, fiber-integrated micro- and nano-cavity setups, for the on-demand production of optical number-state, $0N$-state, and binomial-code-state pulses. The scheme makes use of Raman transitions within an entire atomic ground-state hyperfine level and operates with laser and cavity fields detuned from the atomic transition by much more than the excited-state hyperfine splitting. This enables reduction of the dynamics to that of a simple, cavity-damped Tavis-Cummings model with the collective spin determined by the total angular momentum of the ground hyperfine level. △ Less

Submitted 30 November, 2020; originally announced December 2020.

Comments: 14 pages, 10 figures

Journal ref: Phys. Rev. Lett. 127, 033602 (2021)

Asymptotic expansions of weighted prime power counting functions

Authors: Jesse Elliott

Abstract: We prove several asymptotic continued fraction expansions of $π(x)$, $Π(x)$, $\operatorname{li}(x)$, $\operatorname{Ri}(x)$, and related functions, where $π(x)$ is the prime counting function, $Π(x) = \sum_{k = 1}^\infty \frac{1}{k}π(\sqrt[k]{x})$ is the Riemann prime counting function, and $\operatorname{Ri}(x) = \sum_{k=1}^\infty \frac{ μ(k)}{k} \operatorname{li}(\sqrt[k]{x})$ is Riemann's appro… ▽ More We prove several asymptotic continued fraction expansions of $π(x)$, $Π(x)$, $\operatorname{li}(x)$, $\operatorname{Ri}(x)$, and related functions, where $π(x)$ is the prime counting function, $Π(x) = \sum_{k = 1}^\infty \frac{1}{k}π(\sqrt[k]{x})$ is the Riemann prime counting function, and $\operatorname{Ri}(x) = \sum_{k=1}^\infty \frac{ μ(k)}{k} \operatorname{li}(\sqrt[k]{x})$ is Riemann's approximation to the prime counting function. We also determine asymptotic continued fraction expansions of the function $\sum_{p \leq x} p^s$ for all $s \in \mathbb{C}$ with $\operatorname{Re}(s) > -1$, and of the functions $\sum_{a^x < p \leq a^{x+1}} \frac{1}{p}$ and $\log \prod_{a^x < p \leq a^{x+1}} (1 -1/p)^{-1}$ for all real numbers $a > 1$. We also determine the first few terms of an asymptotic continued fraction expansion of the function $π(ax)-π(bx)$ for $a > b > 0$. As a corollary of these results, we determine the best rational approximations of the "linearized" verions of these various functions. △ Less

Submitted 4 October, 2020; originally announced October 2020.

Comments: 26 pages, 4 figures

MSC Class: 11N05; 30B70; 44A15

Group actions, power mean orbit size, and musical scales

Authors: Jesse Elliott

Abstract: We provide an application of the theory of group actions to the study of musical scales. For any group $G$, finite $G$-set $S$, and real number $t$, we define the {\it $t$-power diameter} $\operatorname{diam}_t(G,S)$ to be the size of any maximal orbit of $S$ divided by the $t$-power mean orbit size of the elements of $S$. The symmetric group $S_{11}$ acts on the set of all tonic scales, where a {… ▽ More We provide an application of the theory of group actions to the study of musical scales. For any group $G$, finite $G$-set $S$, and real number $t$, we define the {\it $t$-power diameter} $\operatorname{diam}_t(G,S)$ to be the size of any maximal orbit of $S$ divided by the $t$-power mean orbit size of the elements of $S$. The symmetric group $S_{11}$ acts on the set of all tonic scales, where a {\it tonic scale} is a subset of $\mathbb{Z}_{12}$ containing $0$. We show that, for all $t \in [-1,1]$, among all the subgroups $G$ of $S_{11}$, the $t$-power diameter of the $G$-set of all heptatonic scales is largest for the subgroup $Γ$, and its conjugate subgroups, generated by $\{(1 \ 2),(3 \ 4),(5 \ 6),(8 \ 9),(10 \ 11)\}$. The unique maximal $Γ$-orbit consists of the 32 thāts of Hindustani classical music popularized by Bhatkhande. This analysis provides a reason why these 32 scales, among all 462 heptatonic scales, are of mathematical interest. We also apply our analysis, to a lesser degree, to hexatonic and pentatonic scales. △ Less

Submitted 4 October, 2020; originally announced October 2020.

Comments: 20 pages

MSC Class: 05E18; 26E60

Flexibility of brain regions during working memory curtails cognitive consequences to lack of sleep

Authors: Nina Lauharatanahirun, Kanika Bansal, Steven M. Thurman, Jean M. Vettel, Barry Giesbrecht, Scott Grafton, James C. Elliott, Erin Flynn-Evans, Emily Falk, Javier O. Garcia

Abstract: Previous research has shown a clear relationship between sleep and memory, examining the impact of sleep deprivation on key cognitive processes over very short durations or in special populations. Here, we show, in a longitudinal 16 week study, that naturalistic, unfettered sleep modulations in healthy adults have significant impacts on the brain. Using a dynamic networks approach combined with hi… ▽ More Previous research has shown a clear relationship between sleep and memory, examining the impact of sleep deprivation on key cognitive processes over very short durations or in special populations. Here, we show, in a longitudinal 16 week study, that naturalistic, unfettered sleep modulations in healthy adults have significant impacts on the brain. Using a dynamic networks approach combined with hierarchical statistical modelling, we show that the flexibility of particular brain regions that span a large network including regions in occipital, temporal, and frontal cortex increased when participants performed a working memory task following low sleep episodes. Critically, performance itself did not change as a function of sleep, implying adaptability in brain networks to compensate for having a poor night's sleep by recruiting the necessary resources to complete the task. We further explore whether this compensatory effect is driven by a (i) increase in the recruitment of network resources over time and/or (ii) an expansion of the network itself. Our results add to the literature linking sleep and memory, provide an analytical framework in which to investigate compensatory modulations in the brain, and highlight the brain's resilience to day-to-day fluctuations of external pressures to performance. △ Less

Submitted 15 September, 2020; originally announced September 2020.

A QM/MD coupling method to model the ion-induced polarization of graphene

Authors: Joshua D Elliott, Alessandro Troisi, Paola Carbone

Abstract: We report a new Quantum Mechanical/Molecular Dynamics (QM/MD) simulation loop to model the coupling between the electron and atom dynamics in solid/liquid interfacial systems. The method can describe simultaneously both the quantum mechanical surface polarizability emerging from the proximity to the electrolyte, and the electrolyte structure and dynamics. In the current set up Density Functional T… ▽ More We report a new Quantum Mechanical/Molecular Dynamics (QM/MD) simulation loop to model the coupling between the electron and atom dynamics in solid/liquid interfacial systems. The method can describe simultaneously both the quantum mechanical surface polarizability emerging from the proximity to the electrolyte, and the electrolyte structure and dynamics. In the current set up Density Functional Tight Binding calculations for the electronic structure calculations of the surface are coupled with classical molecular dynamics to simulate the electrolyte solution. The reduced computational cost of the QM part makes the coupling with a classical simulation engine computationally feasible and allows simulation of large systems for hundreds of nanoseconds. We tested the method by simulating both a non-charged graphene flake and a non-charged and charged infinite graphene sheet immersed in an NaCl electrolyte solution. We found that, when no bias is applied, ions preferentially remained in solution and only cations are mildly attracted to the surface of the graphene. This preferential adsorption of cations vs anions seems to persist also when the surface is moderately charged and rules out any substantial ions/surface charge transfer. △ Less

Submitted 19 June, 2020; v1 submitted 26 March, 2020; originally announced March 2020.

Comments: 22 pages, 8 figures

Journal ref: Journal of Chemical Theory and Computation 2020 16 8 5253

Harmonic numbers and the prime counting function

Authors: Jesse Elliott

Abstract: We provide approximations to the prime counting function by various discretized versions of the logarithmic integral function, expressed solely in terms of the harmonic numbers. We demonstrate with explicit error bounds that these approximations are at least as good as the logarithmic integral approximation. As a corollary, we provide some reformulations of the Riemann hypothesis in terms of the p… ▽ More We provide approximations to the prime counting function by various discretized versions of the logarithmic integral function, expressed solely in terms of the harmonic numbers. We demonstrate with explicit error bounds that these approximations are at least as good as the logarithmic integral approximation. As a corollary, we provide some reformulations of the Riemann hypothesis in terms of the prime counting function and the harmonic numbers. △ Less

Submitted 3 January, 2021; v1 submitted 6 February, 2020; originally announced February 2020.

Comments: 31 pages

MSC Class: 11N05; 11M26

Koopmans Meets Bethe-Salpeter: Excitonic Optical Spectra without GW

Authors: Joshua Elliott, Nicola Colonna, Margherita Marsili, Nicola Marzari, Paolo Umari

Abstract: The Bethe-Salpeter Equation (BSE) can be applied to compute from first-principles optical spectra that include the effects of screened electron-hole interactions. As input, BSE calculations require single-particle states, quasiparticle energy levels and the screened Coulomb interaction, which are typically obtained with many-body perturbation theory, whose cost limits the scope of possible applica… ▽ More The Bethe-Salpeter Equation (BSE) can be applied to compute from first-principles optical spectra that include the effects of screened electron-hole interactions. As input, BSE calculations require single-particle states, quasiparticle energy levels and the screened Coulomb interaction, which are typically obtained with many-body perturbation theory, whose cost limits the scope of possible applications. This work tries to address this practical limitation, instead deriving spectral energies from Koopmans-compliant functionals and introducing a new methodology for handling the screened Coulomb interaction. The explicit calculation of the $W$ matrix is bypassed via a direct minimization scheme applied on top of a maximally localised Wannier function basis. We validate and benchmark this approach by computing the low-lying excited states of the molecules in Thiel's set, and the optical absorption spectrum of a $\text{C}_{60}$ fullerene. The results show the same trends as quantum chemical methods and are in excellent agreement with previous simulations carried out at the TD-DFT or $G_{0}W_{0}$-\text{BSE} level. Conveniently, the new framework reduces the parameter space controlling the accuracy of the calculation, thereby simplifying the simulation of charge-neutral excitations, offering the potential to expand the applicability of first-principles spectroscopies to larger systems of applied interest. △ Less

Submitted 31 December, 2019; originally announced December 2019.

Journal ref: J. Chem. Theory Comput. 2019, 15, 6, 3710-3720

Robust inference of memory structure for efficient quantum modelling of stochastic processes

Authors: Matthew Ho, Mile Gu, Thomas J. Elliott

Abstract: A growing body of work has established the modelling of stochastic processes as a promising area of application for quantum techologies; it has been shown that quantum models are able to replicate the future statistics of a stochastic process whilst retaining less information about the past than any classical model must -- even for a purely classical process. Such memory-efficient models open a po… ▽ More A growing body of work has established the modelling of stochastic processes as a promising area of application for quantum techologies; it has been shown that quantum models are able to replicate the future statistics of a stochastic process whilst retaining less information about the past than any classical model must -- even for a purely classical process. Such memory-efficient models open a potential future route to study complex systems in greater detail than ever before, and suggest profound consequences for our notions of structure in their dynamics. Yet, to date methods for constructing these quantum models are based on having a prior knowledge of the optimal classical model. Here, we introduce a protocol for blind inference of the memory structure of quantum models -- tailored to take advantage of quantum features -- direct from time-series data, in the process highlighting the robustness of their structure to noise. This in turn provides a way to construct memory-efficient quantum models of stochastic processes whilst circumventing certain drawbacks that manifest solely as a result of classical information processing in classical inference protocols. △ Less

Submitted 20 February, 2020; v1 submitted 7 November, 2019; originally announced November 2019.

Comments: 10 pages, 6 figures

Journal ref: Phys. Rev. A 101, 032327 (2020)

Measures of distinguishability between stochastic processes

Authors: Chengran Yang, Felix C. Binder, Mile Gu, Thomas J. Elliott

Abstract: Quantifying how distinguishable two stochastic processes are lies at the heart of many fields, such as machine learning and quantitative finance. While several measures have been proposed for this task, none have universal applicability and ease of use. In this Letter, we suggest a set of requirements for a well-behaved measure of process distinguishability. Moreover, we propose a family of measur… ▽ More Quantifying how distinguishable two stochastic processes are lies at the heart of many fields, such as machine learning and quantitative finance. While several measures have been proposed for this task, none have universal applicability and ease of use. In this Letter, we suggest a set of requirements for a well-behaved measure of process distinguishability. Moreover, we propose a family of measures, called divergence rates, that satisfy all of these requirements. Focussing on a particular member of this family -- the co-emission divergence rate -- we show that it can be computed efficiently, behaves qualitatively similar to other commonly-used measures in their regimes of applicability, and remains well-behaved in scenarios where other measures break down. △ Less

Submitted 18 September, 2019; originally announced September 2019.

Journal ref: Phys. Rev. E 101, 062137 (2020)

Extreme dimensionality reduction with quantum modelling

Authors: Thomas J. Elliott, Chengran Yang, Felix C. Binder, Andrew J. P. Garner, Jayne Thompson, Mile Gu

Abstract: Effective and efficient forecasting relies on identification of the relevant information contained in past observations -- the predictive features -- and isolating it from the rest. When the future of a process bears a strong dependence on its behaviour far into the past, there are many such features to store, necessitating complex models with extensive memories. Here, we highlight a family of sto… ▽ More Effective and efficient forecasting relies on identification of the relevant information contained in past observations -- the predictive features -- and isolating it from the rest. When the future of a process bears a strong dependence on its behaviour far into the past, there are many such features to store, necessitating complex models with extensive memories. Here, we highlight a family of stochastic processes whose minimal classical models must devote unboundedly many bits to tracking the past. For this family, we identify quantum models of equal accuracy that can store all relevant information within a single two-dimensional quantum system (qubit). This represents the ultimate limit of quantum compression and highlights an immense practical advantage of quantum technologies for the forecasting and simulation of complex systems. △ Less

Submitted 23 December, 2020; v1 submitted 6 September, 2019; originally announced September 2019.

Comments: 6+3 pages, 3+1 figures

Journal ref: Phys. Rev. Lett. 125, 260501 (2020)

Semantic interoperability and characterization of data provenance in computational molecular engineering

Authors: M. T. Horsch, C. Niethammer, G. Boccardo, P. Carbone, S. Chiacchiera, M. Chiricotto, J. D. Elliott, V. Lobaskin, P. Neumann, P. Schiffels, M. A. Seaton, I. T. Todorov, J. Vrabec, W. L. Cavalcanti

Abstract: By introducing a common representational system for metadata that describe the employed simulation workflows, diverse sources of data and platforms in computational molecular engineering, such as workflow management systems, can become interoperable at the semantic level. To achieve semantic interoperability, the present work introduces two ontologies that provide a formal specification of the ent… ▽ More By introducing a common representational system for metadata that describe the employed simulation workflows, diverse sources of data and platforms in computational molecular engineering, such as workflow management systems, can become interoperable at the semantic level. To achieve semantic interoperability, the present work introduces two ontologies that provide a formal specification of the entities occurring in a simulation workflow and the relations between them: The software ontology VISO is developed to represent software packages and their features, and OSMO, an ontology for simulation, modelling, and optimization, is introduced on the basis of MODA, a previously developed semi-intuitive graph notation for workflows in materials modelling. As a proof of concept, OSMO is employed to describe a use case of the TaLPas workflow management system, a scheduler and workflow optimizer for particle-based simulations. △ Less

Submitted 15 November, 2019; v1 submitted 29 July, 2019; originally announced August 2019.

Definitions of local density in density-dependent potentials for mixtures

Authors: Peter Vanya, James A. Elliott

Abstract: Density-dependent potentials are frequently used in materials simulations due to their approximate description of many-body effects at minimal computational cost. However, in order to apply such models to multi-component systems, an appropriate definition of total local particle density is required. Here, we discuss two definitions of local density in the context of many-body dissipative particle… ▽ More Density-dependent potentials are frequently used in materials simulations due to their approximate description of many-body effects at minimal computational cost. However, in order to apply such models to multi-component systems, an appropriate definition of total local particle density is required. Here, we discuss two definitions of local density in the context of many-body dissipative particle dynamics. We show that only a potential which combines local densities from all particle types in its argument gives physically meaningful results for all composition ratios. Drawing on the ideas from metal potentials, we redefine local density such that it can accommodate different inter-type interactions despite the constraint to keep the main interaction parameter constant, known as Warren's no-go theorem, and generalise the many-body potential to heterogeneous systems. We then show via simulation how liquid-liquid and liquid-solid coexistence can arise just by tuning the interaction parameters. △ Less

Submitted 10 September, 2020; v1 submitted 31 July, 2019; originally announced July 2019.

Comments: 5 pages, 5 figures, published version

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

Semimartingale dynamics and estimation for a semi-Markov chain

Authors: Robert J. Elliott

Abstract: We consider a finite state discrete time process X. Without loss of generality the finite state space can be identified with the set of unit vectors {e1, e2, . . . , eN} with ei = (0, . . . , 0, 1, 0, . . . , 0)0 2 RN. For a Markov chain the times the process stays in any state are geometrically distributed. This condition is relaxed for a semi-Markov chain. We first derive the semimartingale dyna… ▽ More We consider a finite state discrete time process X. Without loss of generality the finite state space can be identified with the set of unit vectors {e1, e2, . . . , eN} with ei = (0, . . . , 0, 1, 0, . . . , 0)0 2 RN. For a Markov chain the times the process stays in any state are geometrically distributed. This condition is relaxed for a semi-Markov chain. We first derive the semimartingale dynamics for a semi-Markov chain. We then consider the situation where the chain is observed in noise. We suggest how to estimate the occupation times in the states and derive filters and smoothers for quantities associated with the chain. △ Less

Submitted 30 April, 2019; originally announced May 2019.

Surveying structural complexity in quantum many-body systems

Authors: Whei Yeap Suen, Thomas J. Elliott, Jayne Thompson, Andrew J. P. Garner, John R. Mahoney, Vlatko Vedral, Mile Gu

Abstract: Quantum many-body systems exhibit a rich and diverse range of exotic behaviours, owing to their underlying non-classical structure. These systems present a deep structure beyond those that can be captured by measures of correlation and entanglement alone. Using tools from complexity science, we characterise such structure. We investigate the structural complexities that can be found within the pat… ▽ More Quantum many-body systems exhibit a rich and diverse range of exotic behaviours, owing to their underlying non-classical structure. These systems present a deep structure beyond those that can be captured by measures of correlation and entanglement alone. Using tools from complexity science, we characterise such structure. We investigate the structural complexities that can be found within the patterns that manifest from the observational data of these systems. In particular, using two prototypical quantum many-body systems as test cases - the one-dimensional quantum Ising and Bose-Hubbard models - we explore how different information-theoretic measures of complexity are able to identify different features of such patterns. This work furthers the understanding of fully-quantum notions of structure and complexity in quantum systems and dynamics. △ Less

Submitted 18 March, 2022; v1 submitted 23 December, 2018; originally announced December 2018.

Comments: 9 pages, 5 figures

Managing uncertainty in data-derived densities to accelerate density functional theory

Authors: Andrew T. Fowler, Chris J. Pickard, James A. Elliott

Abstract: Faithful representations of atomic environments and general models for regression can be harnessed to learn electron densities that are close to the ground state. One of the applications of data-derived electron densities is to orbital-free density functional theory. However, extrapolations of densities learned from a training set to dissimilar structures could result in inaccurate results, which… ▽ More Faithful representations of atomic environments and general models for regression can be harnessed to learn electron densities that are close to the ground state. One of the applications of data-derived electron densities is to orbital-free density functional theory. However, extrapolations of densities learned from a training set to dissimilar structures could result in inaccurate results, which would limit the applicability of the method. Here, we show that a non-Bayesian approach can produce estimates of uncertainty which can successfully distinguish accurate from inaccurate predictions of electron density. We apply our approach to density functional theory where we initialise calculations with data-derived densities only when we are confident about their quality. This results in a guaranteed acceleration to self-consistency for configurations that are similar to those seen during training and could be useful for sampling based methods, where previous ground state densities cannot be used to initialise subsequent calculations. △ Less

Submitted 28 February, 2019; v1 submitted 5 December, 2018; originally announced December 2018.

Optimal stochastic modelling with unitary quantum dynamics

Authors: Qing Liu, Thomas. J. Elliott, Felix. C. Binder, Carlo Di Franco, Mile Gu

Abstract: Identifying and extracting the past information relevant to the future behaviour of stochastic processes is a central task in the quantitative sciences. Quantum models offer a promising approach to this, allowing for accurate simulation of future trajectories whilst using less past information than any classical counterpart. Here we introduce a class of phase-enhanced quantum models, representing… ▽ More Identifying and extracting the past information relevant to the future behaviour of stochastic processes is a central task in the quantitative sciences. Quantum models offer a promising approach to this, allowing for accurate simulation of future trajectories whilst using less past information than any classical counterpart. Here we introduce a class of phase-enhanced quantum models, representing the most general means of causal simulation with a unitary quantum circuit. We show that the resulting constructions can display advantages over previous state-of-art methods - both in the amount of information they need to store about the past, and in the minimal memory dimension they require to store this information. Moreover, we find that these two features are generally competing factors in optimisation - leading to an ambiguity in what constitutes the optimal model - a phenomenon that does not manifest classically. Our results thus simultaneously offer new quantum advantages for stochastic simulation, and illustrate further qualitative differences in behaviour between classical and quantum notions of complexity. △ Less

Submitted 23 October, 2018; originally announced October 2018.

Comments: 9 pages, 5 figures

Journal ref: Phys. Rev. A 99, 062110 (2019)

Asymptotic expansions of the prime counting function

Authors: Jesse Elliott

Abstract: We provide several asymptotic expansions of the prime counting function $π(x)$ and related functions. We define an {\it asymptotic continued fraction expansion} of a complex-valued function of a real or complex variable to be a possibly divergent continued fraction whose approximants provide an asymptotic expansion of the given function. We show that, for each positive integer $n$, two well-known… ▽ More We provide several asymptotic expansions of the prime counting function $π(x)$ and related functions. We define an {\it asymptotic continued fraction expansion} of a complex-valued function of a real or complex variable to be a possibly divergent continued fraction whose approximants provide an asymptotic expansion of the given function. We show that, for each positive integer $n$, two well-known continued fraction expansions of the exponential integral function $E_n(z)$ correspondingly yield two asymptotic continued fraction expansions of $π(x)/x$. We prove this by first establishing some general results about asymptotic continued fraction expansions. We show, for instance, that the "best"' rational function approximations of a function possessing an asymptotic Jacobi continued fraction expansion are precisely the approximants of the continued fraction, and as a corollary we determine all of the best rational function approximations of the function $π(e^x)/e^x$. Finally, we generalize our results on $π(x)$ to any arithmetic semigroup satisfying Axiom A, and thus to any number field. △ Less

Submitted 17 August, 2021; v1 submitted 18 September, 2018; originally announced September 2018.

Comments: 35 pages, 2 figures

MSC Class: 11N05; 11N45; 11N80; 11R42; 30B70; 44A15

Relativistic supernova 2009bb exploded close to an atomic gas cloud

Authors: Michał J. Michałowski, G. Gentile, T. Kruhler, H. Kuncarayakti, P. Kamphuis, J. Hjorth, S. Berta, V. D'Elia, J. Elliott, L. Galbany, J. Greiner, L. K. Hunt, M. P. Koprowski, E. Le Floc'h, A. Nicuesa Guelbenzu, E. Palazzi, J. Rasmussen, A. Rossi, S. Savaglio, A. de Ugarte Postigo, P. van der Werf, S. D. Vergani

Abstract: The potential similarity of the powering mechanisms of relativistic SNe and GRBs allowed us to make a prediction that relativistic SNe are born in environments similar to those of GRBs, that is, ones which are rich in atomic gas. Here we embark on testing this hypothesis by analysing the properties of the host galaxy NGC 3278 of the relativistic SN 2009bb. This is the first time the atomic gas pro… ▽ More The potential similarity of the powering mechanisms of relativistic SNe and GRBs allowed us to make a prediction that relativistic SNe are born in environments similar to those of GRBs, that is, ones which are rich in atomic gas. Here we embark on testing this hypothesis by analysing the properties of the host galaxy NGC 3278 of the relativistic SN 2009bb. This is the first time the atomic gas properties of a relativistic SN host are provided and the first time resolved 21 cm-hydrogen-line (HI) information is provided for a host of an SN of any type in the context of the SN position. We obtained radio observations with ATCA covering the HI line, and optical integral field unit spectroscopy observations with MUSE. The atomic gas distribution of NGC 3278 is not centred on the optical galaxy centre, but instead around a third of atomic gas resides in the region close to the SN position. This galaxy has a few times lower atomic and molecular gas masses than predicted from its SFR. SN 2009bb exploded close to the region with the highest SFR density and the lowest age (~5.5 Myr, corresponding to the initial mass of the progenitor star ~36 Mo). As for GRB hosts, the gas properties of NGC 3278 are consistent with a recent inflow of gas from the intergalactic medium, which explains the concentration of atomic gas close to the SN position and the enhanced SFR. Super-solar metallicity at the position of the SN (unlike for most GRBs) may mean that relativistic explosions signal a recent inflow of gas (and subsequent star formation), and their type (GRBs or SNe) is determined either i) by the metallicity of the inflowing gas, so that metal-poor gas results in a GRB explosion and metal-rich gas results in a relativistic SN explosion without an accompanying GRB, or ii) by the efficiency of gas mixing, or iii) by the type of the galaxy. △ Less

Submitted 2 August, 2018; originally announced August 2018.

Comments: Accepted by A&A. 18 pages, 9 figures, 7 tables. The ATCA and VLT/MUSE data will be available through CDS linked in the published version

Journal ref: A&A 618, A104 (2018)

Invariance of experimental observables with respect to coarse-graining in standard and many-body dissipative particle dynamics

Authors: Peter Vanya, Jonathan Sharman, James A. Elliott

Abstract: Dissipative particle dynamics (DPD) is a well-established mesoscale simulation method. However, there have been long-standing ambiguities regarding the dependence of its (purely repulsive) force field parameter on temperature as well as the variation of the resulting experimental observables, such as diffusivity or surface tension, with coarse-graining (CG) degree. Here, we revisit the role of the… ▽ More Dissipative particle dynamics (DPD) is a well-established mesoscale simulation method. However, there have been long-standing ambiguities regarding the dependence of its (purely repulsive) force field parameter on temperature as well as the variation of the resulting experimental observables, such as diffusivity or surface tension, with coarse-graining (CG) degree. Here, we revisit the role of the CG degree and rederive the temperature dependence in standard DPD simulations. Consequently, we derive a scaling of the input variables that renders the system properties invariant with respect to CG degree, and illustrate the versatility of the method by computing the surface tensions of binary solvent mixtures. We then extend this procedure to many-body dissipative particle dynamics (MDPD) and, by computing surface tensions of the same mixtures at a range of CG degrees, demonstrate that this newer method, which has not been widely applied so far, is also capable of simulating complex fluids of practical interest. △ Less

Submitted 21 February, 2019; v1 submitted 11 May, 2018; originally announced May 2018.

Comments: 12 pages, 7 figures, published version

Journal ref: J. Chem. Phys. 150, 064101 (2019)

Liquid phase parametrisation and solidification in many-body dissipative particle dynamics

Authors: Peter Vanya, Phillip Crout, Jonathan Sharman, James A. Elliott

Abstract: Many-body dissipative particle dynamics (MDPD) is a mesoscale method capable of reproducing liquid-vapour coexistence in a single simulation. Despite having been introduced more than a decade ago, this method remains broadly unexplored and, as a result, relatively unused for modelling of industrially important soft matter systems. In this work, we systematically investigate the structure and prope… ▽ More Many-body dissipative particle dynamics (MDPD) is a mesoscale method capable of reproducing liquid-vapour coexistence in a single simulation. Despite having been introduced more than a decade ago, this method remains broadly unexplored and, as a result, relatively unused for modelling of industrially important soft matter systems. In this work, we systematically investigate the structure and properties of an MDPD fluid. We show that, besides the liquid phase, the MDPD potential can also yield a gas phase and a thermodynamically stable solid phase with a bcc lattice, but lacking a proper stress-strain relation. For the liquid phase, we determine the dependence of density and surface tension on the interaction parameters, and devise a top-down parametrisation protocol for real liquids. △ Less

Submitted 21 February, 2019; v1 submitted 11 May, 2018; originally announced May 2018.

Comments: 8 pages, 6 figures, updated tables, published version

Journal ref: Phys. Rev. E 98, 033310 (2018)