-
On picking operations in e-commerce warehouses: Insights from the complete-information counterpart
Authors:
Catherine Lorenz,
Alena Otto,
Michel Gendreau
Abstract:
Major players in e-commerce process dynamically incoming orders in real-time and already use advanced anticipation techniques, like AI, to predict characteristics of future orders. However, at the warehousing level, there are still no unambiguous recommendations on integrating anticipation with intelligent online optimization algorithms, nor an unbiased benchmark to assess the improvement potentia…
▽ More
Major players in e-commerce process dynamically incoming orders in real-time and already use advanced anticipation techniques, like AI, to predict characteristics of future orders. However, at the warehousing level, there are still no unambiguous recommendations on integrating anticipation with intelligent online optimization algorithms, nor an unbiased benchmark to assess the improvement potential of advanced anticipation over myopic techniques, as optimal online solutions are usually unavailable. In this paper, we compute and analyze Complete-Information Optimal policy Solutions (CIOSs), of an exact perfect anticipation algorithm with full knowledge of future customer orders' arrival times and items, for picking operations in picker-to-parts warehouses. We provide analytical properties and leverage CIOSs to uncover decision patterns that enhance simpler algorithms, improving both makespan (costs) and order turnover (delivery speed). Using a metric similar to the gap to optimality, we quantify the gains from optimization elements and the improvements remaining for advanced anticipation mechanisms. To compute CIOSs, we design the first exact algorithm for the Order Batching, Sequencing, and Routing Problem with Release Times (OBSRP-R), based on dynamic programming. Our analysis advises on the largely overlooked intervention - the dynamic adjustment of started batches, and the strategic relocation of an idle picker towards future picking locations. The former affected over 60% of CIOS orders and triggered consistent improvements across warehousing policies. The latter occurred before 39%-62% CIOS orders and could decrease a myopic policy's gap to optimal turnover by 4.3% on average (up to 14%). Notably, we challenge the debated concept of strategic waiting, revealing why the latter resembles an "all-in" gamble and harms both makespan and order turnover when intervention is allowed.
△ Less
Submitted 18 October, 2024;
originally announced October 2024.
-
Picking Operations in Warehouses with Dynamically Arriving Orders: How Good is Reoptimization?
Authors:
Catherine Lorenz,
Alena Otto,
Michel Gendreau
Abstract:
E-commerce operations are essentially online, with customer orders arriving dynamically. However, very little is known about the performance of online policies for warehousing with respect to optimality, particularly for order picking and batching operations, which constitute a substantial portion of the total operating costs in warehouses. We aim to close this gap for one of the most prominent dy…
▽ More
E-commerce operations are essentially online, with customer orders arriving dynamically. However, very little is known about the performance of online policies for warehousing with respect to optimality, particularly for order picking and batching operations, which constitute a substantial portion of the total operating costs in warehouses. We aim to close this gap for one of the most prominent dynamic algorithms, namely reoptimization (Reopt), which reoptimizes the current solution each time when a new order arrives. We examine Reopt in the Online Order Batching, Sequencing, and Routing Problem (OOBSRP), in both cases when the picker uses either a manual pushcart or a robotic cart. Moreover, we examine the noninterventionist Reopt in the case of a manual pushcart, wherein picking instructions are provided exclusively at the depot. We establish analytical performance bounds employing worst-case and probabilistic analysis. We demonstrate that, under generic stochastic assumptions, Reopt is almost surely asymptotically optimal and, notably, we validate its near-optimal performance in computational experiments across a broad range of warehouse settings. These results underscore Reopt's relevance as a method for online warehousing applications.
△ Less
Submitted 19 September, 2024;
originally announced September 2024.
-
Digital Ecosystem for FAIR Time Series Data Management in Environmental System Science
Authors:
J. Bumberger,
M. Abbrent,
N. Brinckmann,
J. Hemmen,
R. Kunkel,
C. Lorenz,
P. Lünenschloß,
B. Palm,
T. Schnicke,
C. Schulz,
H. van der Schaaf,
D. Schäfer
Abstract:
Addressing the challenges posed by climate change, biodiversity loss, and environmental pollution requires comprehensive monitoring and effective data management strategies that are applicable across various scales in environmental system science. This paper introduces a versatile and transferable digital ecosystem for managing time series data, designed to adhere to the FAIR principles (Findable,…
▽ More
Addressing the challenges posed by climate change, biodiversity loss, and environmental pollution requires comprehensive monitoring and effective data management strategies that are applicable across various scales in environmental system science. This paper introduces a versatile and transferable digital ecosystem for managing time series data, designed to adhere to the FAIR principles (Findable, Accessible, Interoperable, and Reusable). The system is highly adaptable, cloud-ready, and suitable for deployment in a wide range of settings, from small-scale projects to large-scale monitoring initiatives. The ecosystem comprises three core components: the Sensor Management System (SMS) for detailed metadata registration and management; time$.$IO, a platform for efficient time series data storage, transfer, and real-time visualization; and the System for Automated Quality Control (SaQC), which ensures data integrity through real-time analysis and quality assurance. The modular architecture, combined with standardized protocols and interfaces, ensures that the ecosystem can be easily transferred and deployed across different environments and institutions. This approach enhances data accessibility for a broad spectrum of stakeholders, including researchers, policymakers, and the public, while fostering collaboration and advancing scientific research in environmental monitoring.
△ Less
Submitted 17 September, 2024; v1 submitted 5 September, 2024;
originally announced September 2024.
-
Insider trading in discrete time Kyle games
Authors:
Christoph Kühn,
Christopher Lorenz
Abstract:
We present a new discrete time version of Kyle's (1985) classic model of insider trading, formulated as a generalised extensive form game. The model has three kinds of traders: an insider, random noise traders, and a market maker. The insider aims to exploit her informational advantage and maximise expected profits while the market maker observes the total order flow and sets prices accordingly.…
▽ More
We present a new discrete time version of Kyle's (1985) classic model of insider trading, formulated as a generalised extensive form game. The model has three kinds of traders: an insider, random noise traders, and a market maker. The insider aims to exploit her informational advantage and maximise expected profits while the market maker observes the total order flow and sets prices accordingly.
First, we show how the multi-period model with finitely many pure strategies can be reduced to a (static) social system in the sense of Debreu (1952) and prove the existence of a sequential Kyle equilibrium, following Kreps and Wilson (1982). This works for any probability distribution with finite support of the noise trader's demand and the true value, and for any finite information flow of the insider. In contrast to Kyle (1985) with normal distributions, equilibria exist in general only in mixed strategies and not in pure strategies.
In the single-period model we establish bounds for the insider's strategy in equilibrium. Finally, we prove the existence of an equilibrium for the game with a continuum of actions, by considering an approximating sequence of games with finitely many actions. Because of the lack of compactness of the set of measurable price functions, standard infinite-dimensional fixed point theorems are not applicable.
△ Less
Submitted 1 July, 2024; v1 submitted 1 December, 2023;
originally announced December 2023.
-
Elastic Microphase Separation Produces Robust Bicontinuous Materials
Authors:
Carla Fernández-Rico,
Sanjay Schreiber,
Hamza Oudich,
Charlotta Lorenz,
Alba Sicher,
Tianqi Sai,
Stefanie Heyden,
Pietro Carrara,
Laura De Lorenzis,
Robert W. Style,
Eric R. Dufresne
Abstract:
Bicontinuous microstructures are essential to the function of diverse natural and synthetic systems. Their synthesis has been based on two approaches: arrested phase separation or self-assembly of block copolymers. The former is attractive for its chemical simplicity, the latter for its thermodynamic robustness. Here, we introduce Elastic MicroPhase Separation (EMPS) as an alternative approach to…
▽ More
Bicontinuous microstructures are essential to the function of diverse natural and synthetic systems. Their synthesis has been based on two approaches: arrested phase separation or self-assembly of block copolymers. The former is attractive for its chemical simplicity, the latter for its thermodynamic robustness. Here, we introduce Elastic MicroPhase Separation (EMPS) as an alternative approach to make bicontinuous microstructures. Conceptually, EMPS balances the molecular-scale forces that drive demixing with large-scale elasticity to encode a thermodynamic length scale. This process features a continuous phase transition, reversible without hysteresis. Practically, we trigger EMPS by simply super-saturating an elastomeric matrix with a liquid. This results in uniform bicontinuous materials with a well-defined microscopic length-scale tuned by the matrix stiffness. The versatility and robustness of EMPS is further demonstrated by fabricating bicontinuous materials with superior mechanical properties and controlled anisotropy and microstructural gradients.
△ Less
Submitted 22 April, 2023;
originally announced April 2023.
-
Pathology Synthesis of 3D-Consistent Cardiac MR Images using 2D VAEs and GANs
Authors:
Sina Amirrajab,
Cristian Lorenz,
Juergen Weese,
Josien Pluim,
Marcel Breeuwer
Abstract:
We propose a method for synthesizing cardiac magnetic resonance (MR) images with plausible heart pathologies and realistic appearances for the purpose of generating labeled data for the application of supervised deep-learning (DL) training. The image synthesis consists of label deformation and label-to-image translation tasks. The former is achieved via latent space interpolation in a VAE model, w…
▽ More
We propose a method for synthesizing cardiac magnetic resonance (MR) images with plausible heart pathologies and realistic appearances for the purpose of generating labeled data for the application of supervised deep-learning (DL) training. The image synthesis consists of label deformation and label-to-image translation tasks. The former is achieved via latent space interpolation in a VAE model, while the latter is accomplished via a label-conditional GAN model. We devise three approaches for label manipulation in the latent space of the trained VAE model; i) \textbf{intra-subject synthesis} aiming to interpolate the intermediate slices of a subject to increase the through-plane resolution, ii) \textbf{inter-subject synthesis} aiming to interpolate the geometry and appearance of intermediate images between two dissimilar subjects acquired with different scanner vendors, and iii) \textbf{pathology synthesis} aiming to synthesize a series of pseudo-pathological synthetic subjects with characteristics of a desired heart disease. Furthermore, we propose to model the relationship between 2D slices in the latent space of the VAE prior to reconstruction for generating 3D-consistent subjects from stacking up 2D slice-by-slice generations. We demonstrate that such an approach could provide a solution to diversify and enrich an available database of cardiac MR images and to pave the way for the development of generalizable DL-based image analysis algorithms. We quantitatively evaluate the quality of the synthesized data in an augmentation scenario to achieve generalization and robustness to multi-vendor and multi-disease data for image segmentation. Our code is available at https://github.com/sinaamirrajab/CardiacPathologySynthesis.
△ Less
Submitted 30 May, 2023; v1 submitted 9 September, 2022;
originally announced September 2022.
-
sim2real: Cardiac MR Image Simulation-to-Real Translation via Unsupervised GANs
Authors:
Sina Amirrajab,
Yasmina Al Khalil,
Cristian Lorenz,
Jurgen Weese,
Josien Pluim,
Marcel Breeuwer
Abstract:
There has been considerable interest in the MR physics-based simulation of a database of virtual cardiac MR images for the development of deep-learning analysis networks. However, the employment of such a database is limited or shows suboptimal performance due to the realism gap, missing textures, and the simplified appearance of simulated images. In this work we 1) provide image simulation on vir…
▽ More
There has been considerable interest in the MR physics-based simulation of a database of virtual cardiac MR images for the development of deep-learning analysis networks. However, the employment of such a database is limited or shows suboptimal performance due to the realism gap, missing textures, and the simplified appearance of simulated images. In this work we 1) provide image simulation on virtual XCAT subjects with varying anatomies, and 2) propose sim2real translation network to improve image realism. Our usability experiments suggest that sim2real data exhibits a good potential to augment training data and boost the performance of a segmentation algorithm.
△ Less
Submitted 9 August, 2022;
originally announced August 2022.
-
Very large-scale neighborhood search for drone routing with energy replenishment
Authors:
Catherine Lorenz,
Nicola Mimmo,
Alena Otto,
Daniele Vigo
Abstract:
The Drone Routing Problem with Energy replenishment (DRP-E) belongs to a general class of routing problems with intermediate stops and synchronization constraints. In DRP-E, the drone has to visit a set of nodes and routinely requires battery swaps from a (potentially) mobile replenishment station. Contrary to widespread restrictions in the drone routing literature, several destinations may be vis…
▽ More
The Drone Routing Problem with Energy replenishment (DRP-E) belongs to a general class of routing problems with intermediate stops and synchronization constraints. In DRP-E, the drone has to visit a set of nodes and routinely requires battery swaps from a (potentially) mobile replenishment station. Contrary to widespread restrictions in the drone routing literature, several destinations may be visited in between two consecutive battery swaps. In this paper, we propose a nontrivial very large-scale neighbourhood for DRP-E, which synergetically leverages two large-sized polynomially solvable DRP-E SubProblems (SP1 and SP2). The number of feasible solutions in the resulting neighborhood is a multiple of those in SP1 and SP2, and, thus, exponential in the input size of the problem, whereas the computational time to search it remains polynomial. The proposed polynomial two-stage dynamic programming algorithm VLSN to search this neighborhood can be flexibly adjusted to the desired trade-off between accuracy and computational time. For instance, the search procedure can be converted into an exact algorithm of competitive runtime for DRP-E. In computational tests, the developed solution methods outperform current state-of-the art heuristics for DRP-E by a significant margin. A case study based on a search for missing persons demonstrates that VLSN easily accommodates additional practice relevant features and outperforms the state-of-the-art solution in disaster relief by 20%.
△ Less
Submitted 28 July, 2022;
originally announced July 2022.
-
sympy2c: from symbolic expressions to fast C/C++ functions and ODE solvers in Python
Authors:
Uwe Schmitt,
Beatrice Moser,
Christiane S. Lorenz,
Alexandre Refregier
Abstract:
Computer algebra systems play an important role in science as they facilitate the development of new theoretical models. The resulting symbolic equations are often implemented in a compiled programming language in order to provide fast and portable codes for practical applications. We describe sympy2c, a new Python package designed to bridge the gap between the symbolic development and the numeric…
▽ More
Computer algebra systems play an important role in science as they facilitate the development of new theoretical models. The resulting symbolic equations are often implemented in a compiled programming language in order to provide fast and portable codes for practical applications. We describe sympy2c, a new Python package designed to bridge the gap between the symbolic development and the numerical implementation of a theoretical model. sympy2c translates symbolic equations implemented in the SymPy Python package to C/C++ code that is optimized using symbolic transformations. The resulting functions can be conveniently used as an extension module in Python. sympy2c is used within the PyCosmo Python package to solve the Einstein-Boltzmann equations, a large system of ODEs describing the evolution of linear perturbations in the Universe. After reviewing the functionalities and usage of sympy2c, we describe its implementation and optimization strategies. This includes, in particular, a novel approach to generate optimized ODE solvers making use of the sparsity of the symbolic Jacobian matrix. We demonstrate its performance using the Einstein-Boltzmann equations as a test case. sympy2c is widely applicable and may prove useful for various areas of computational physics. sympy2c is publicly available at https://cosmology.ethz.ch/research/software-lab/sympy2c.html
△ Less
Submitted 22 March, 2022;
originally announced March 2022.
-
Synergy between cosmological and laboratory searches in neutrino physics
Authors:
Martina Gerbino,
Evan Grohs,
Massimiliano Lattanzi,
Kevork N. Abazajian,
Nikita Blinov,
Thejs Brinckmann,
Mu-Chun Chen,
Zelimir Djurcic,
Peizhi Du,
Miguel Escudero,
Steffen Hagstotz,
Kevin J. Kelly,
Christiane S. Lorenz,
Marilena Loverde,
Pablo Martínez-Miravé,
Olga Mena,
Joel Meyers,
Walter Pettus,
Ninetta Saviano,
Anna M. Suliga,
Volodymyr Takhistov,
Mariam Tórtola,
José W. F. Valle,
Benjamin Wallisch
Abstract:
The intersection of the cosmic and neutrino frontiers is a rich field where much discovery space still remains. Neutrinos play a pivotal role in the hot big bang cosmology, influencing the dynamics of the universe over numerous decades in cosmological history. Recent studies have made tremendous progress in understanding some properties of cosmological neutrinos, primarily their energy density. Up…
▽ More
The intersection of the cosmic and neutrino frontiers is a rich field where much discovery space still remains. Neutrinos play a pivotal role in the hot big bang cosmology, influencing the dynamics of the universe over numerous decades in cosmological history. Recent studies have made tremendous progress in understanding some properties of cosmological neutrinos, primarily their energy density. Upcoming cosmological probes will measure the energy density of relativistic particles with higher precision, but could also start probing other properties of the neutrino spectra. When convolved with results from terrestrial experiments, cosmology can become even more acute at probing new physics related to neutrinos or even Beyond the Standard Model (BSM). Any discordance between laboratory and cosmological data sets may reveal new BSM physics and/or suggest alternative models of cosmology. We give examples of the intersection between terrestrial and cosmological probes in the neutrino sector, and briefly discuss the possibilities of what different laboratory experiments may see in conjunction with cosmological observatories.
△ Less
Submitted 26 December, 2023; v1 submitted 14 March, 2022;
originally announced March 2022.
-
Symbolic Implementation of Extensions of the $\texttt{PyCosmo}$ Boltzmann Solver
Authors:
Beatrice Moser,
Christiane S. Lorenz,
Uwe Schmitt,
Alexandre Refregier,
Janis Fluri,
Raphael Sgier,
Federica Tarsitano,
Lavinia Heisenberg
Abstract:
$\texttt{PyCosmo}$ is a Python-based framework for the fast computation of cosmological model predictions. One of its core features is the symbolic representation of the Einstein-Boltzmann system of equations. Efficient $\texttt{C/C++}$ code is generated from the $\texttt{SymPy}$ symbolic expressions making use of the $\texttt{sympy2c}…
▽ More
$\texttt{PyCosmo}$ is a Python-based framework for the fast computation of cosmological model predictions. One of its core features is the symbolic representation of the Einstein-Boltzmann system of equations. Efficient $\texttt{C/C++}$ code is generated from the $\texttt{SymPy}$ symbolic expressions making use of the $\texttt{sympy2c}$ package. This enables easy extensions of the equation system for the implementation of new cosmological models. We illustrate this with three extensions of the $\texttt{PyCosmo}$ Boltzmann solver to include a dark energy component with a constant equation of state, massive neutrinos and a radiation streaming approximation. We describe the $\texttt{PyCosmo}$ framework, highlighting new features, and the symbolic implementation of the new models. We compare the $\texttt{PyCosmo}$ predictions for the $Λ$CDM model extensions with $\texttt{CLASS}$, both in terms of accuracy and computational speed. We find a good agreement, to better than 0.1% when using high-precision settings and a comparable computational speed. Links to the Python Package Index (PyPI) page of the code release and to the PyCosmo Hub, an online platform where the package is installed, are available at: https://cosmology.ethz.ch/research/software-lab/PyCosmo.html.
△ Less
Submitted 17 June, 2022; v1 submitted 15 December, 2021;
originally announced December 2021.
-
EuCAPT White Paper: Opportunities and Challenges for Theoretical Astroparticle Physics in the Next Decade
Authors:
R. Alves Batista,
M. A. Amin,
G. Barenboim,
N. Bartolo,
D. Baumann,
A. Bauswein,
E. Bellini,
D. Benisty,
G. Bertone,
P. Blasi,
C. G. Böhmer,
Ž. Bošnjak,
T. Bringmann,
C. Burrage,
M. Bustamante,
J. Calderón Bustillo,
C. T. Byrnes,
F. Calore,
R. Catena,
D. G. Cerdeño,
S. S. Cerri,
M. Chianese,
K. Clough,
A. Cole,
P. Coloma
, et al. (112 additional authors not shown)
Abstract:
Astroparticle physics is undergoing a profound transformation, due to a series of extraordinary new results, such as the discovery of high-energy cosmic neutrinos with IceCube, the direct detection of gravitational waves with LIGO and Virgo, and many others. This white paper is the result of a collaborative effort that involved hundreds of theoretical astroparticle physicists and cosmologists, und…
▽ More
Astroparticle physics is undergoing a profound transformation, due to a series of extraordinary new results, such as the discovery of high-energy cosmic neutrinos with IceCube, the direct detection of gravitational waves with LIGO and Virgo, and many others. This white paper is the result of a collaborative effort that involved hundreds of theoretical astroparticle physicists and cosmologists, under the coordination of the European Consortium for Astroparticle Theory (EuCAPT). Addressed to the whole astroparticle physics community, it explores upcoming theoretical opportunities and challenges for our field of research, with particular emphasis on the possible synergies among different subfields, and the prospects for solving the most fundamental open questions with multi-messenger observations.
△ Less
Submitted 19 October, 2021;
originally announced October 2021.
-
Combined $13\times2$-point analysis of the Cosmic Microwave Background and Large-Scale Structure: implications for the $S_8$-tension and neutrino mass constraints
Authors:
Raphael Sgier,
Christiane Lorenz,
Alexandre Refregier,
Janis Fluri,
Dominik Zürcher,
Federica Tarsitano
Abstract:
We present cosmological constraints for the flat $Λ$CDM model, including the sum of neutrino masses, by performing a multi-probe analysis of a total of 13 tomographic auto- and cross-angular power spectra. This is achieved by combining, at map level, the latest primary CMB and CMB-lensing measurements from the Planck 2018 data release, as well as spectroscopic galaxy samples from BOSS DR12, and th…
▽ More
We present cosmological constraints for the flat $Λ$CDM model, including the sum of neutrino masses, by performing a multi-probe analysis of a total of 13 tomographic auto- and cross-angular power spectra. This is achieved by combining, at map level, the latest primary CMB and CMB-lensing measurements from the Planck 2018 data release, as well as spectroscopic galaxy samples from BOSS DR12, and the latest Kilo-Degree Survey (KiDS-1000) tomographic weak lensing shear data release. Our analysis includes auto- and cross-correlations as well as calibration parameters for all cosmological probes, thus providing a self-calibration of the combined data sets. We find a good fit (reduced $χ^2$=1.7) for the combined probes with calibration parameters only moderately different from their nominal value, thus giving a possible interpretation of the tension between the early- and late-Universe probes. The resulting value for the structure growth parameter is $S_8 = 0.754 \pm 0.016$ (68\% CL). We also obtain a $\sim$2.3$σ$ constraint on the neutrino mass sum of $\sum m_ν= 0.51^{+0.21}_{-0.24}$ eV (68\% CL), which is compatible with current particle physics limits. We perform several tests by fixing the neutrino mass sum to a low value, considering narrower priors on the multiplicative bias parameters for cosmic shear, and by fixing all calibration parameters to their expected values. These tests result in worse fits compared to our fiducial run, especially for the case when all calibration parameters are fixed. This latter test also yields a lower upper limit of the neutrino mass sum. We discuss how the interplay between the cosmological and calibration parameters impact the $S_8$-tension and the constraints on the neutrino mass sum. [abridged]
△ Less
Submitted 7 October, 2021;
originally announced October 2021.
-
Impact of climate change on West Nile virus distribution in South America
Authors:
Camila Lorenz,
Thiago Salomao de Azevedo,
Francisco Chiaravalloti-Neto
Abstract:
West Nile virus (WNV) is a vector-borne pathogen of global relevance and is currently the most widely distributed flavivirus of encephalitis worldwide. This virus infects birds, humans, horses, and other mammals, and its transmission cycle occurs in urban and rural areas. Climate conditions have direct and indirect impacts on vector abundance and virus dynamics within the mosquito. The significanc…
▽ More
West Nile virus (WNV) is a vector-borne pathogen of global relevance and is currently the most widely distributed flavivirus of encephalitis worldwide. This virus infects birds, humans, horses, and other mammals, and its transmission cycle occurs in urban and rural areas. Climate conditions have direct and indirect impacts on vector abundance and virus dynamics within the mosquito. The significance of environmental variables as drivers in WNV epidemiology is increasing under the current climate change scenario. In this study, we used a machine learning algorithm to model WNV distributions in South America. Our model evaluated eight environmental variables (type of biome, annual temperature, seasonality of temperature, daytime temperature variation, thermal amplitude, seasonality of precipitation, annual rainfall, and elevation) for their contribution to the occurrence of WNV since its introduction in South America (2004). Our results showed that environmental variables can directly alter the occurrence of WNV, with lower precipitation and higher temperatures associated with increased virus incidence. High-risk areas may be modified in the coming years, becoming more evident with high greenhouse gas emission levels. Countries such as Bolivia and Paraguay will be greatly affected, drastically changing their current WNV distribution. Several Brazilian areas will also increase the likelihood of presenting WNV, mainly in the Northeast and Midwest regions and the Pantanal biome. The Galapagos Islands will also probably increase their geographic range suitable for WNV occurrence. It is necessary to develop preventive policies to minimize potential WNV infection in humans and enhance active epidemiological surveillance in birds, humans, and other mammals before it becomes a more significant public health problem in South America.
△ Less
Submitted 1 April, 2021;
originally announced April 2021.
-
Effect of social isolation in dengue cases in the state of Sao Paulo, Brazil: an analysis during the COVID-19 pandemic
Authors:
Gleice Margarete de Souza Conceição,
Gerson Laurindo Barbosa,
Camila Lorenz,
Ana Carolina Dias Bocewicz,
Lidia Maria Reis Santana,
Cristiano Corrêa de Azevedo Marques,
Francisco Chiaravalloti-Neto
Abstract:
Background: Studies have shown that human mobility is an important factor in dengue epidemiology. Changes in mobility resulting from COVID-19 pandemic set up a real-life situation to test this hypothesis. Our objective was to evaluate the effect of reduced mobility due to this pandemic in the occurrence of dengue in the state of São Paulo, Brazil. Method: It is an ecological study of time series,…
▽ More
Background: Studies have shown that human mobility is an important factor in dengue epidemiology. Changes in mobility resulting from COVID-19 pandemic set up a real-life situation to test this hypothesis. Our objective was to evaluate the effect of reduced mobility due to this pandemic in the occurrence of dengue in the state of São Paulo, Brazil. Method: It is an ecological study of time series, developed between January and August 2020. We use the number of confirmed dengue cases and residential mobility, on a daily basis, from secondary information sources. Mobility was represented by the daily percentage variation of residential population isolation, obtained from the Google database. We modeled the relationship between dengue occurrence and social distancing by negative binomial regression, adjusted for seasonality. We represent the social distancing dichotomously (isolation versus no isolation) and consider lag for isolation from the dates of occurrence of dengue. Results: The risk of dengue decreased around 9.1% (95% CI: 14.2 to 3.7) in the presence of isolation, considering a delay of 20 days between the degree of isolation and the dengue first symptoms. Conclusions: We have shown that mobility can play an important role in the epidemiology of dengue and should be considered in surveillance and control activities
△ Less
Submitted 15 March, 2021;
originally announced March 2021.
-
Examining socioeconomic factors to understand the hospital case-fatality rates of COVID-19 in the city of Sao Paulo, Brazil
Authors:
Camila Lorenz,
Patricia Marques Moralejo Bermudi,
Marcelo Antunes Failla,
Breno Souza de Aguiar,
Tatiana Natasha Toporcov,
Francisco Chiaravalloti Neto,
Ligia Vizeu Barrozo
Abstract:
Understanding differences in hospital case-fatality rates (HCFRs) of coronavirus disease (COVID-19) may help evaluate its severity and the capacity of the healthcare system to reduce mortality. We examined the variability in HCFRs of COVID-19 in relation to spatial inequalities in socioeconomic factors across the city of Sao Paulo, Brazil. We found that HCFRs were higher for men and for individual…
▽ More
Understanding differences in hospital case-fatality rates (HCFRs) of coronavirus disease (COVID-19) may help evaluate its severity and the capacity of the healthcare system to reduce mortality. We examined the variability in HCFRs of COVID-19 in relation to spatial inequalities in socioeconomic factors across the city of Sao Paulo, Brazil. We found that HCFRs were higher for men and for individuals aged 60 years and older. Our models identified per capita income as a significant factor that is negatively associated with the HCFRs of COVID-19, even after adjusting by age, sex and presence of risk factors.
△ Less
Submitted 28 February, 2021;
originally announced March 2021.
-
Reconstruction of the neutrino mass as a function of redshift
Authors:
Christiane S. Lorenz,
Lena Funcke,
Matthias Löffler,
Erminia Calabrese
Abstract:
We reconstruct the neutrino mass as a function of redshift, z, from current cosmological data using both standard binned priors and linear spline priors with variable knots. Using cosmic microwave background temperature, polarization and lensing data, in combination with distance measurements from baryonic acoustic oscillations and supernovae, we find that the neutrino mass is consistent with…
▽ More
We reconstruct the neutrino mass as a function of redshift, z, from current cosmological data using both standard binned priors and linear spline priors with variable knots. Using cosmic microwave background temperature, polarization and lensing data, in combination with distance measurements from baryonic acoustic oscillations and supernovae, we find that the neutrino mass is consistent with $\sum m_ν(z)$ = const. We obtain a larger bound on the neutrino mass at low redshifts coinciding with the onset of dark energy domination, $\sum m_ν(z = 0)$ < 1.46 eV (95% CL). This result can be explained either by the well-known degeneracy between $\sum m_ν$ and $Ω_Λ$ at low redshifts, or by models in which neutrino masses are generated very late in the Universe. We finally convert our results into cosmological limits for models with non-relativistic neutrino decay and find $\sum m_ν$ < 0.21 eV (95% CL), which would be out of reach for the KATRIN experiment.
△ Less
Submitted 29 November, 2021; v1 submitted 26 February, 2021;
originally announced February 2021.
-
Spatiotemporal dynamic of COVID-19 mortality in the city of Sao Paulo, Brazil: shifting the high risk from the best to the worst socio-economic conditions
Authors:
Patricia Marques Moralejo Bermudi,
Camila Lorenz,
Breno Souza de Aguiar,
Marcelo Antunes Failla,
Ligia Vizeu Barrozo,
Francisco Chiaravalloti-Neto
Abstract:
Currently, Brazil has one of the fastest increasing COVID-19 epidemics in the world, that has caused at least 94 thousand confirmed deaths until now. The city of Sao Paulo is particularly vulnerable because it is the most populous in the country. Analyzing the spatiotemporal dynamics of COVID-19 is important to help the urgent need to integrate better actions to face the pandemic. Thus, this study…
▽ More
Currently, Brazil has one of the fastest increasing COVID-19 epidemics in the world, that has caused at least 94 thousand confirmed deaths until now. The city of Sao Paulo is particularly vulnerable because it is the most populous in the country. Analyzing the spatiotemporal dynamics of COVID-19 is important to help the urgent need to integrate better actions to face the pandemic. Thus, this study aimed to analyze the COVID-19 mortality, from March to July 2020, considering the spatio-time architectures, the socio-economic context of the population, and using a fine granular level, in the most populous city in Brazil. For this, we conducted an ecological study, using secondary public data from the mortality information system. We describe mortality rates for each epidemiological week and the entire period by sex and age. We modelled the deaths using spatiotemporal and spatial architectures and Poisson probability distributions in a latent Gaussian Bayesian model approach. We obtained the relative risks for temporal and spatiotemporal trends and socio-economic conditions. To reduce possible sub notification, we considered the confirmed and suspected deaths. Our findings showed an apparent stabilization of the temporal trend, at the end of the period, but that may change in the future. Mortality rate increased with increasing age and was higher in men. The risk of death was greater in areas with the worst social conditions throughout the study period. However, this was not a uniform pattern over time, since we identified a shift from the high risk in the areas with best socio-economic conditions to the worst ones. Our study contributed by emphasizing the importance of geographic screening in areas with a higher risk of death, and, currently, worse socio-economic contexts, as a crucial aspect to reducing disease mortality and health inequities, through integrated public health actions.
△ Less
Submitted 5 August, 2020;
originally announced August 2020.
-
XCAT-GAN for Synthesizing 3D Consistent Labeled Cardiac MR Images on Anatomically Variable XCAT Phantoms
Authors:
Sina Amirrajab,
Samaneh Abbasi-Sureshjani,
Yasmina Al Khalil,
Cristian Lorenz,
Juergen Weese,
Josien Pluim,
Marcel Breeuwer
Abstract:
Generative adversarial networks (GANs) have provided promising data enrichment solutions by synthesizing high-fidelity images. However, generating large sets of labeled images with new anatomical variations remains unexplored. We propose a novel method for synthesizing cardiac magnetic resonance (CMR) images on a population of virtual subjects with a large anatomical variation, introduced using th…
▽ More
Generative adversarial networks (GANs) have provided promising data enrichment solutions by synthesizing high-fidelity images. However, generating large sets of labeled images with new anatomical variations remains unexplored. We propose a novel method for synthesizing cardiac magnetic resonance (CMR) images on a population of virtual subjects with a large anatomical variation, introduced using the 4D eXtended Cardiac and Torso (XCAT) computerized human phantom. We investigate two conditional image synthesis approaches grounded on a semantically-consistent mask-guided image generation technique: 4-class and 8-class XCAT-GANs. The 4-class technique relies on only the annotations of the heart; while the 8-class technique employs a predicted multi-tissue label map of the heart-surrounding organs and provides better guidance for our conditional image synthesis. For both techniques, we train our conditional XCAT-GAN with real images paired with corresponding labels and subsequently at the inference time, we substitute the labels with the XCAT derived ones. Therefore, the trained network accurately transfers the tissue-specific textures to the new label maps. By creating 33 virtual subjects of synthetic CMR images at the end-diastolic and end-systolic phases, we evaluate the usefulness of such data in the downstream cardiac cavity segmentation task under different augmentation strategies. Results demonstrate that even with only 20% of real images (40 volumes) seen during training, segmentation performance is retained with the addition of synthetic CMR images. Moreover, the improvement in utilizing synthetic images for augmenting the real data is evident through the reduction of Hausdorff distance up to 28% and an increase in the Dice score up to 5%, indicating a higher similarity to the ground truth in all dimensions.
△ Less
Submitted 31 July, 2020; v1 submitted 27 July, 2020;
originally announced July 2020.
-
Integrative dynamic structural biology unveils conformers essential for the oligomerization of a large GTPase
Authors:
Thomas-Otavio Peulen,
Carola S. Hengstenberg,
Ralf Biehl,
Mykola Dimura,
Charlotte Lorenz,
Alessandro Valeri,
Semra Ince,
Tobias Vöpel,
Bela Faragó,
Holger Gohlke,
Johann P. Klare,
Andreas M. Stadler,
Claus A. M. Seidel,
Christian Herrmann
Abstract:
Guanylate binding proteins (GBPs) are soluble dynamin-like proteins with structured domains that undergo a conformational transition for GTP-controlled oligomerization to exert their function as part of the innate immune system of mammalian cells - attacking intra-cellular parasites by disrupting their membranes. The structural basis and mechanism of this process is unknown. Therefore, we apply ne…
▽ More
Guanylate binding proteins (GBPs) are soluble dynamin-like proteins with structured domains that undergo a conformational transition for GTP-controlled oligomerization to exert their function as part of the innate immune system of mammalian cells - attacking intra-cellular parasites by disrupting their membranes. The structural basis and mechanism of this process is unknown. Therefore, we apply neutron spin echo, X-ray scattering, fluorescence, and EPR spectroscopy as techniques for integrative dynamic structural biology to human GBP1 (hGBP1). We mapped hGBP1's essential dynamics from nanoseconds to milliseconds by motional spectra of sub-domains. We find a GTP-independent flexibility of the C-terminal effector domain in the $μ$s-regime and structurally characterize conformers being essential that hGBP1 can open like a pocketknife for oligomerization. This unveils the intrinsic flexibility, a GTP-triggered association of the GTPase-domains and assembly-dependent GTP-hydrolysis as functional design principles of hGBP1 that control its reversible oligomerization in polar assemblies and the subsequent formation of condensates.
△ Less
Submitted 8 April, 2020;
originally announced April 2020.
-
4D Semantic Cardiac Magnetic Resonance Image Synthesis on XCAT Anatomical Model
Authors:
Samaneh Abbasi-Sureshjani,
Sina Amirrajab,
Cristian Lorenz,
Juergen Weese,
Josien Pluim,
Marcel Breeuwer
Abstract:
We propose a hybrid controllable image generation method to synthesize anatomically meaningful 3D+t labeled Cardiac Magnetic Resonance (CMR) images. Our hybrid method takes the mechanistic 4D eXtended CArdiac Torso (XCAT) heart model as the anatomical ground truth and synthesizes CMR images via a data-driven Generative Adversarial Network (GAN). We employ the state-of-the-art SPatially Adaptive De…
▽ More
We propose a hybrid controllable image generation method to synthesize anatomically meaningful 3D+t labeled Cardiac Magnetic Resonance (CMR) images. Our hybrid method takes the mechanistic 4D eXtended CArdiac Torso (XCAT) heart model as the anatomical ground truth and synthesizes CMR images via a data-driven Generative Adversarial Network (GAN). We employ the state-of-the-art SPatially Adaptive De-normalization (SPADE) technique for conditional image synthesis to preserve the semantic spatial information of ground truth anatomy. Using the parameterized motion model of the XCAT heart, we generate labels for 25 time frames of the heart for one cardiac cycle at 18 locations for the short axis view. Subsequently, realistic images are generated from these labels, with modality-specific features that are learned from real CMR image data. We demonstrate that style transfer from another cardiac image can be accomplished by using a style encoder network. Due to the flexibility of XCAT in creating new heart models, this approach can result in a realistic virtual population to address different challenges the medical image analysis research community is facing such as expensive data collection. Our proposed method has a great potential to synthesize 4D controllable CMR images with annotations and adaptable styles to be used in various supervised multi-site, multi-vendor applications in medical image analysis.
△ Less
Submitted 20 May, 2020; v1 submitted 17 February, 2020;
originally announced February 2020.
-
Open Power System Data - Frictionless data for electricity system modelling
Authors:
Frauke Wiese,
Ingmar Schlecht,
Wolf-Dieter Bunke,
Clemens Gerbaulet,
Lion Hirth,
Martin Jahn,
Friedrich Kunz,
Casimir Lorenz,
Jonathan Mühlenpfordt,
Juliane Reimann,
Wolf-Peter Schill
Abstract:
The quality of electricity system modelling heavily depends on the input data used. Although a lot of data is publicly available, it is often dispersed, tedious to process and partly contains errors. We argue that a central provision of input data for modelling has the character of a public good: it reduces overall societal costs for quantitative energy research as redundant work is avoided, and i…
▽ More
The quality of electricity system modelling heavily depends on the input data used. Although a lot of data is publicly available, it is often dispersed, tedious to process and partly contains errors. We argue that a central provision of input data for modelling has the character of a public good: it reduces overall societal costs for quantitative energy research as redundant work is avoided, and it improves transparency and reproducibility in electricity system modelling. This paper describes the Open Power System Data platform that aims at realising the efficiency and quality gains of centralised data provision by collecting, checking, processing, aggregating, documenting and publishing data required by most modellers. We conclude that the platform can provide substantial benefits to energy system analysis by raising efficiency of data pre-processing, providing a method for making data pre-processing for energy system modelling traceable, flexible and reproducible and improving the quality of original data published by data providers.
△ Less
Submitted 3 January, 2019; v1 submitted 10 December, 2018;
originally announced December 2018.
-
Core Cosmology Library: Precision Cosmological Predictions for LSST
Authors:
Nora Elisa Chisari,
David Alonso,
Elisabeth Krause,
C. Danielle Leonard,
Philip Bull,
Jérémy Neveu,
Antonia Villarreal,
Sukhdeep Singh,
Thomas McClintock,
John Ellison,
Zilong Du,
Joe Zuntz,
Alexander Mead,
Shahab Joudaki,
Christiane S. Lorenz,
Tilman Troester,
Javier Sanchez,
Francois Lanusse,
Mustapha Ishak,
Renée Hlozek,
Jonathan Blazek,
Jean-Eric Campagne,
Husni Almoubayyed,
Tim Eifler,
Matthew Kirby
, et al. (5 additional authors not shown)
Abstract:
The Core Cosmology Library (CCL) provides routines to compute basic cosmological observables to a high degree of accuracy, which have been verified with an extensive suite of validation tests. Predictions are provided for many cosmological quantities, including distances, angular power spectra, correlation functions, halo bias and the halo mass function through state-of-the-art modeling prescripti…
▽ More
The Core Cosmology Library (CCL) provides routines to compute basic cosmological observables to a high degree of accuracy, which have been verified with an extensive suite of validation tests. Predictions are provided for many cosmological quantities, including distances, angular power spectra, correlation functions, halo bias and the halo mass function through state-of-the-art modeling prescriptions available in the literature. Fiducial specifications for the expected galaxy distributions for the Large Synoptic Survey Telescope (LSST) are also included, together with the capability of computing redshift distributions for a user-defined photometric redshift model. A rigorous validation procedure, based on comparisons between CCL and independent software packages, allows us to establish a well-defined numerical accuracy for each predicted quantity. As a result, predictions for correlation functions of galaxy clustering, galaxy-galaxy lensing and cosmic shear are demonstrated to be within a fraction of the expected statistical uncertainty of the observables for the models and in the range of scales of interest to LSST. CCL is an open source software package written in C, with a python interface and publicly available at https://github.com/LSSTDESC/CCL.
△ Less
Submitted 5 April, 2019; v1 submitted 14 December, 2018;
originally announced December 2018.
-
Time-varying neutrino mass from a supercooled phase transition: current cosmological constraints and impact on the $Ω_m$-$σ_8$ plane
Authors:
Christiane S. Lorenz,
Lena Funcke,
Erminia Calabrese,
Steen Hannestad
Abstract:
In this paper we investigate a time-varying neutrino mass model, motivated by the mild tension between cosmic microwave background (CMB) measurements of the matter fluctuations and those obtained from low-redshift data. We modify the minimal case of the model proposed by Dvali and Funcke (2016) that predicts late neutrino mass generation in a post-recombination cosmic phase transition, by assuming…
▽ More
In this paper we investigate a time-varying neutrino mass model, motivated by the mild tension between cosmic microwave background (CMB) measurements of the matter fluctuations and those obtained from low-redshift data. We modify the minimal case of the model proposed by Dvali and Funcke (2016) that predicts late neutrino mass generation in a post-recombination cosmic phase transition, by assuming that neutrino asymmetries allow for the presence of relic neutrinos in the late-time Universe. We show that, if the transition is supercooled, current cosmological data (including CMB temperature, polarization and lensing, baryon acoustic oscillations, and Type Ia supernovae) prefer the scale factor $a_s$ of the phase transition to be very large, peaking at $a_s\sim 1$, and therefore supporting a cosmological scenario in which neutrinos are almost massless until very recent times. We find that in this scenario the cosmological bound on the total sum of the neutrino masses today is significantly weakened compared to the standard case of constant-mass neutrinos, with $\sum m_ν<4.8$~eV at 95\% confidence, and in agreement with the model predictions. The main reason for this weaker bound is a large correlation arising between the dark energy and neutrino components in the presence of false vacuum energy that converts into the non-zero neutrino masses after the transition. This result provides new targets for the coming KATRIN and PTOLEMY experiments. We also show that the time-varying neutrino mass model considered here does not provide a clear explanation to the existing cosmological $Ω_m$-$σ_8$ discrepancies.
△ Less
Submitted 14 December, 2018; v1 submitted 5 November, 2018;
originally announced November 2018.
-
Deep Learning Based Rib Centerline Extraction and Labeling
Authors:
Matthias Lenga,
Tobias Klinder,
Christian Bürger,
Jens von Berg,
Astrid Franz,
Cristian Lorenz
Abstract:
Automated extraction and labeling of rib centerlines is a typically needed prerequisite for more advanced assisted reading tools that help the radiologist to efficiently inspect all 24 ribs in a CT volume. In this paper, we combine a deep learning-based rib detection with a dedicated centerline extraction algorithm applied to the detection result for the purpose of fast, robust and accurate rib ce…
▽ More
Automated extraction and labeling of rib centerlines is a typically needed prerequisite for more advanced assisted reading tools that help the radiologist to efficiently inspect all 24 ribs in a CT volume. In this paper, we combine a deep learning-based rib detection with a dedicated centerline extraction algorithm applied to the detection result for the purpose of fast, robust and accurate rib centerline extraction and labeling from CT volumes. More specifically, we first apply a fully convolutional neural network (FCNN) to generate a probability map for detecting the first rib pair, the twelfth rib pair, and the collection of all intermediate ribs. In a second stage, a newly designed centerline extraction algorithm is applied to this multi-label probability map. Finally, the distinct detection of first and twelfth rib separately, allows to derive individual rib labels by simple sorting and counting the detected centerlines. We applied our method to CT volumes from 116 patients which included a variety of different challenges and achieved a centerline accuracy of 0.787 mm with respect to manual centerline annotations.
This article is a preprint version of: Lenga M., Klinder T., Bürger C., von Berg J., Franz A., Lorenz C. (2019) Deep Learning Based Rib Centerline Extraction and Labeling. In: Vrtovec T., Yao J., Zheng G., Pozo J. (eds) Computational Methods and Clinical Applications in Musculoskeletal Imaging. MSKI 2018. Lecture Notes in Computer Science, vol 11404. Springer, Cham
△ Less
Submitted 14 January, 2019; v1 submitted 19 September, 2018;
originally announced September 2018.
-
Driving Spin-Boson Models From Equilibrium Using Exact Quantum Dynamics
Authors:
G. M. G. McCaul,
C. D. Lorenz,
L. Kantorovich
Abstract:
We present an application of the Extended Stochastic Liouville Equation (ESLE) Phys. Rev. B 95, 125124, which gives an exact solution for the reduced density matrix of an open system surrounded by a harmonic heat bath. This method considers the extended system (the open system and the bath) being thermally equilibrated prior to the action of a time dependent perturbation, as opposed to the usual a…
▽ More
We present an application of the Extended Stochastic Liouville Equation (ESLE) Phys. Rev. B 95, 125124, which gives an exact solution for the reduced density matrix of an open system surrounded by a harmonic heat bath. This method considers the extended system (the open system and the bath) being thermally equilibrated prior to the action of a time dependent perturbation, as opposed to the usual assumption that system and bath are initially partitioned. This is an exact technique capable of accounting for arbitrary parameter regimes of the model. Here we present our first numerical implementation of the method in the simplest case of a Caldeira-Leggett representation of the bath Hamiltonian, and apply it to a spin-Boson system driven from coupled equilibrium. We observe significant behaviours in both the transient dynamics and asymptotic states of the reduced density matrix not present in the usual approximation.
△ Less
Submitted 18 December, 2017;
originally announced December 2017.
-
The impact of relativistic effects on cosmological parameter estimation
Authors:
Christiane S. Lorenz,
David Alonso,
Pedro G. Ferreira
Abstract:
Future surveys will access large volumes of space and hence very long wavelength fluctuations of the matter density and gravitational field. It has been argued that the set of secondary effects that affect the galaxy distribution, relativistic in nature, will bring new, complementary cosmological constraints. We study this claim in detail by focusing on a subset of wide-area future surveys: Stage-…
▽ More
Future surveys will access large volumes of space and hence very long wavelength fluctuations of the matter density and gravitational field. It has been argued that the set of secondary effects that affect the galaxy distribution, relativistic in nature, will bring new, complementary cosmological constraints. We study this claim in detail by focusing on a subset of wide-area future surveys: Stage-4 cosmic microwave background experiments and photometric redshift surveys. In particular, we look at the magnification lensing contribution to galaxy clustering and general relativistic corrections to all observables. We quantify the amount of information encoded in these effects in terms of the tightening of the final cosmological constraints as well as the potential bias in inferred parameters associated with neglecting them. We do so for a wide range of cosmological parameters, covering neutrino masses, standard dark-energy parametrizations and scalar-tensor gravity theories. Our results show that, while the effect of lensing magnification to number counts does not contain a significant amount of information when galaxy clustering is combined with cosmic shear measurements, this contribution does play a significant role in biasing estimates on a host of parameter families if unaccounted for. Since the amplitude of the magnification term is controlled by the slope of the source number counts with apparent magnitude, $s(z)$, we also estimate the accuracy to which this quantity must be known to avoid systematic parameter biases, finding that future surveys will need to determine $s(z)$ to the $\sim$5-10\% level. On the contrary, large-scale general-relativistic corrections are irrelevant both in terms of information content and parameter bias for most cosmological parameters, but significant for the level of primordial non-Gaussianity.
△ Less
Submitted 14 February, 2018; v1 submitted 6 October, 2017;
originally announced October 2017.
-
Distinguishing between Neutrinos and time-varying Dark Energy through Cosmic Time
Authors:
Christiane S. Lorenz,
Erminia Calabrese,
David Alonso
Abstract:
We study the correlations between parameters characterizing neutrino physics and the evolution of dark energy. Using a fluid approach, we show that time-varying dark energy models exhibit degeneracies with the cosmic neutrino background over extended periods of the cosmic history, leading to a degraded estimation of the total mass and number of species of neutrinos. We investigate how to break deg…
▽ More
We study the correlations between parameters characterizing neutrino physics and the evolution of dark energy. Using a fluid approach, we show that time-varying dark energy models exhibit degeneracies with the cosmic neutrino background over extended periods of the cosmic history, leading to a degraded estimation of the total mass and number of species of neutrinos. We investigate how to break degeneracies and combine multiple probes across cosmic time to anchor the behaviour of the two components. We use Planck CMB data and BAO measurements from the BOSS, SDSS and 6dF surveys to present current limits on the model parameters, and then forecast the future reach from the CMB Stage-4 and DESI experiments. We show that a multi-probe analysis of current data provides only marginal improvement on the determination of the individual parameters and no reduction of the correlations. Future observations will better distinguish the neutrino mass and preserve the current sensitivity to the number of species even in case of a time-varying dark energy component.
△ Less
Submitted 25 August, 2017; v1 submitted 2 June, 2017;
originally announced June 2017.
-
Modelling a Bistable System Strongly Coupled to a Debye Bath: A Quasiclassical Approach Based on the Generalised Langevin Equation
Authors:
L. Stella,
H. Ness,
C. D. Lorenz,
L. Kantorovich
Abstract:
Bistable systems present two degenerate metastable configurations separated by an energy barrier. Thermal or quantum fluctuations can promote the transition between the configurations at a rate which depends on the dynamical properties of the local environment (i.e., a thermal bath). In the case of classical systems, strong system-bath interaction has been successfully modelled by the Generalised…
▽ More
Bistable systems present two degenerate metastable configurations separated by an energy barrier. Thermal or quantum fluctuations can promote the transition between the configurations at a rate which depends on the dynamical properties of the local environment (i.e., a thermal bath). In the case of classical systems, strong system-bath interaction has been successfully modelled by the Generalised Langevin Equation (GLE) formalism. Here we show that the efficient GLE algorithm introduced in Phys. Rev. B 89, 134303 (2014) can be extended to include some crucial aspects of the quantum fluctuations. In particular, the expected isotopic effect is observed along with the convergence of the quantum and classical transition rates in the strong coupling limit. Saturation of the transition rates at low temperature is also retrieved, in qualitative, yet not quantitative, agreement with the analytic predictions. The discrepancies in the tunnelling regime are due to an incorrect sampling close to the barrier top. The domain of applicability of the quasiclassical GLE is also discussed.
△ Less
Submitted 12 April, 2017;
originally announced April 2017.
-
High-precision mass measurements for the isobaric multiplet mass equation at A = 52
Authors:
D. A. Nesterenko,
A. Kankainen,
L. Canete,
M. Block,
D. Cox,
T. Eronen,
C. Fahlander,
U. Forsberg,
J. Gerl,
P. Golubev,
J. Hakala,
A. Jokinen,
V. S. Kolhinen,
J. Koponen,
N. Lalović,
Ch. Lorenz,
I. D. Moore,
P. Papadakis,
J. Reinikainen,
S. Rinta-Antila,
D. Rudolph,
L. G. Sarmiento,
A. Voss,
J. Äystö
Abstract:
Masses of $^{52}$Co, $^{52}$Co$^m$, $^{52}$Fe, $^{52}$Fe$^m$, and $^{52}$Mn have been measured with the JYFLTRAP double Penning trap mass spectrometer. Of these, $^{52}$Co and $^{52}$Co$^m$ have been experimentally determined for the first time and found to be more bound than predicted by extrapolations. The isobaric multiplet mass equation for the $T=2$ quintet at $A=52$ has been studied employin…
▽ More
Masses of $^{52}$Co, $^{52}$Co$^m$, $^{52}$Fe, $^{52}$Fe$^m$, and $^{52}$Mn have been measured with the JYFLTRAP double Penning trap mass spectrometer. Of these, $^{52}$Co and $^{52}$Co$^m$ have been experimentally determined for the first time and found to be more bound than predicted by extrapolations. The isobaric multiplet mass equation for the $T=2$ quintet at $A=52$ has been studied employing the new mass values. No significant breakdown (beyond the $3σ$ level) of the quadratic form of the IMME was observed ($χ^2/n=2.4$). The cubic coefficient was 6.0(32) keV ($χ^2/n=1.1$). The excitation energies for the isomer and the $T=2$ isobaric analogue state in $^{52}$Co have been determined to be 374(13) keV and 2922(13) keV, respectively. The $Q$ value for the proton decay from the $19/2^-$ isomer in $^{53}$Co has been determined with an unprecedented precision, $Q_{p} = 1558.8(17)$ keV. The proton separation energies of $^{52}$Co and $^{53}$Ni relevant for the astrophysical rapid proton capture process have been experimentally determined for the first time. \end{abstract}
△ Less
Submitted 15 January, 2017;
originally announced January 2017.
-
Partition-free approach to open quantum systems in harmonic environments: an exact stochastic Liouville equation
Authors:
G. M. G. McCaul,
C. D. Lorenz,
L. Kantorovich
Abstract:
We present a partition-free approach to the evolution of density matrices for open quantum systems coupled to a harmonic environment. The influence functional formalism combined with a two-time Hubbard-Stratonovich transformation allows us to derive a set of exact differential equations for the reduced density matrix of an open system, termed the Extended Stochastic Liouville-von Neumann equation.…
▽ More
We present a partition-free approach to the evolution of density matrices for open quantum systems coupled to a harmonic environment. The influence functional formalism combined with a two-time Hubbard-Stratonovich transformation allows us to derive a set of exact differential equations for the reduced density matrix of an open system, termed the Extended Stochastic Liouville-von Neumann equation. Our approach generalises previous work based on Caldeira-Leggett models and a partitioned initial density matrix. This provides a simple, yet exact, closed-form description for the evolution of open systems from equilibriated initial conditions. The applicability of this model and the potential for numerical implementations are also discussed.
△ Less
Submitted 3 January, 2017; v1 submitted 16 December, 2016;
originally announced December 2016.
-
Nonequilibrium Generalised Langevin Equation for the calculation of heat transport properties in model 1D atomic chains coupled to two 3D thermal baths
Authors:
H. Ness,
L. Stella,
C. D. Lorenz,
L. Kantorovich
Abstract:
We use a Generalised Langevin Equation (GLE) scheme to study the thermal transport of low dimensional systems. In this approach, the central classical region is connected to two realistic thermal baths kept at two different temperatures [H. Ness et al., Phys. Rev. B {\bf 93}, 174303 (2016)]. We consider model Al systems, i.e. one-dimensional atomic chains connected to three-dimensional baths. The…
▽ More
We use a Generalised Langevin Equation (GLE) scheme to study the thermal transport of low dimensional systems. In this approach, the central classical region is connected to two realistic thermal baths kept at two different temperatures [H. Ness et al., Phys. Rev. B {\bf 93}, 174303 (2016)]. We consider model Al systems, i.e. one-dimensional atomic chains connected to three-dimensional baths. The thermal transport properties are studied as a function of the chain length $N$ and the temperature difference $ΔT$ between the baths. We calculate the transport properties both in the linear response regime and in the non-linear regime. Two different laws are obtained for the linear conductance versus the length of the chains. For large temperatures ($T \gtrsim 500$ K) and temperature differences ($ΔT \gtrsim 500$ K), the chains, with $N > 18$ atoms, present a diffusive transport regime with the presence of a temperature gradient across the system. For lower temperatures($T \lesssim 500$ K) and temperature differences ($ΔT \lesssim 400$ K), a regime similar to the ballistic regime is observed. Such a ballistic-like regime is also obtained for shorter chains ($N \le 15 $). Our detailed analysis suggests that the behaviour at higher temperatures and temperature differences is mainly due to anharmonic effects within the long chains.
△ Less
Submitted 8 April, 2017; v1 submitted 3 December, 2016;
originally announced December 2016.
-
c-number Quantum Generalised Langevin Equation for an open system
Authors:
L. Kantorovich,
H. Ness,
L. Stella,
C. Lorenz
Abstract:
We derive a $c-$number Generalised Langevin Equation (GLE) describing the evolution of the expectation values $\left\langle x_{i}\right\rangle_{t}$ of the atomic position operators $x_{i}$ of an open system. The latter is coupled linearly to a harmonic bath kept at a fixed temperature. The equations of motion contain a non-Markovian friction term with the classical kernel {[}L. Kantorovich - PRB 7…
▽ More
We derive a $c-$number Generalised Langevin Equation (GLE) describing the evolution of the expectation values $\left\langle x_{i}\right\rangle_{t}$ of the atomic position operators $x_{i}$ of an open system. The latter is coupled linearly to a harmonic bath kept at a fixed temperature. The equations of motion contain a non-Markovian friction term with the classical kernel {[}L. Kantorovich - PRB 78, 094304 (2008){]} and a zero mean \emph{non-Gaussian} random force with correlation functions that depend on the initial preparation of the open system. We used a density operator formalism without assuming that initially the combined system was decoupled. The only approximation made in deriving quantum GLE consists in assuming that the Hamiltonian of the open system at time $t$ can be expanded up to the second order with respect to operators of atomic displacements $u_{i}=x_{i}-\left\langle x_{i}\right\rangle_{t}$ in the open system around their exact atomic positions $\left\langle x_{i}\right\rangle_{t}$ (the "harmonisation" approximation). The noise is introduced to ensure that sampling many quantum GLE trajectories yields exactly the average one. An explicit expression for the pair correlation function of the noise, consistent with the classical limit, is also proposed. Unlike the usually considered quantum operator GLE, the proposed $c-$number quantum GLE can be used in direct molecular dynamic simulations of open systems under general equilibrium or non-equilibrium conditions.
△ Less
Submitted 8 July, 2016;
originally announced July 2016.
-
Theano: A Python framework for fast computation of mathematical expressions
Authors:
The Theano Development Team,
Rami Al-Rfou,
Guillaume Alain,
Amjad Almahairi,
Christof Angermueller,
Dzmitry Bahdanau,
Nicolas Ballas,
Frédéric Bastien,
Justin Bayer,
Anatoly Belikov,
Alexander Belopolsky,
Yoshua Bengio,
Arnaud Bergeron,
James Bergstra,
Valentin Bisson,
Josh Bleecher Snyder,
Nicolas Bouchard,
Nicolas Boulanger-Lewandowski,
Xavier Bouthillier,
Alexandre de Brébisson,
Olivier Breuleux,
Pierre-Luc Carrier,
Kyunghyun Cho,
Jan Chorowski,
Paul Christiano
, et al. (88 additional authors not shown)
Abstract:
Theano is a Python library that allows to define, optimize, and evaluate mathematical expressions involving multi-dimensional arrays efficiently. Since its introduction, it has been one of the most used CPU and GPU mathematical compilers - especially in the machine learning community - and has shown steady performance improvements. Theano is being actively and continuously developed since 2008, mu…
▽ More
Theano is a Python library that allows to define, optimize, and evaluate mathematical expressions involving multi-dimensional arrays efficiently. Since its introduction, it has been one of the most used CPU and GPU mathematical compilers - especially in the machine learning community - and has shown steady performance improvements. Theano is being actively and continuously developed since 2008, multiple frameworks have been built on top of it and it has been used to produce many state-of-the-art machine learning models.
The present article is structured as follows. Section I provides an overview of the Theano software and its community. Section II presents the principal features of Theano and how to use them, and compares them with other similar projects. Section III focuses on recently-introduced functionalities and improvements. Section IV compares the performance of Theano against Torch7 and TensorFlow on several machine learning models. Section V discusses current limitations of Theano and potential ways of improving it.
△ Less
Submitted 9 May, 2016;
originally announced May 2016.
-
Nonequilibrium processes from Generalised Langevin Equations: realistic nanoscale systems connected to two thermal baths
Authors:
H. Ness,
A. Genina,
L. Stella,
C. D. Lorenz,
L. Kantorovich
Abstract:
We extend the Generalised Langevin Equation (GLE) method [Phys. Rev. B 89, 134303 (2014)] to model a central classical region connected to two realistic thermal baths at two different temperatures. In such nonequilibrium conditions a heat flow is established, via the central system, in between the two baths. The GLE-2B (GLE two baths) scheme permits us to have a realistic description of both the d…
▽ More
We extend the Generalised Langevin Equation (GLE) method [Phys. Rev. B 89, 134303 (2014)] to model a central classical region connected to two realistic thermal baths at two different temperatures. In such nonequilibrium conditions a heat flow is established, via the central system, in between the two baths. The GLE-2B (GLE two baths) scheme permits us to have a realistic description of both the dissipative central system and its surrounding baths. Following the original GLE approach, the extended Langevin dynamics scheme is modified to take into account two sets of auxiliary degrees of freedom corresponding to the mapping of the vibrational properties of each bath. These auxiliary variables are then used to solve the non-Markovian dissipative dynamics of the central region. The resulting algorithm is used to study a model of a short Al nanowire connected to two baths. The results of the simulations using the GLE-2B approach are compared to the results of other simulations that were carried out using standard thermostatting approaches (based on Markovian Langevin and Nose-Hoover thermostats). We concentrate on the steady state regime and study the establishment of a local temperature profile within the system. The conditions for obtaining a flat profile or a temperature gradient are examined in detail, in agreement with earlier studies. The results show that the GLE-2B approach is able to treat, within a single scheme, two widely different thermal transport regimes, i.e. ballistic systems, with no temperature gradient, and diffusive systems with a temperature gradient.
△ Less
Submitted 30 April, 2016; v1 submitted 11 December, 2015;
originally announced December 2015.
-
Applications of the Generalised Langevin Equation: towards a realistic description of the baths
Authors:
H. Ness,
L. Stella,
C. D. Lorenz,
L. Kantorovich
Abstract:
The Generalised Langevin Equation (GLE) method, as developed in Ref. [Phys. Rev. B 89, 134303 (2014)], is used to calculate the dissipative dynamics of systems described at the atomic level. The GLE scheme goes beyond the commonly used bilinear coupling between the central system and the bath, and permits us to have a realistic description of both the dissipative central system and its surrounding…
▽ More
The Generalised Langevin Equation (GLE) method, as developed in Ref. [Phys. Rev. B 89, 134303 (2014)], is used to calculate the dissipative dynamics of systems described at the atomic level. The GLE scheme goes beyond the commonly used bilinear coupling between the central system and the bath, and permits us to have a realistic description of both the dissipative central system and its surrounding bath. We show how to obtain the vibrational properties of a realistic bath and how to convey such properties into an extended Langevin dynamics by the use of the mapping of the bath vibrational properties onto a set of auxiliary variables. Our calculations for a model of a Lennard-Jones solid show that our GLE scheme provides a stable dynamics, with the dissipative/relaxation processes properly described. The total kinetic energy of the central system always thermalises toward the expected bath temperature, with appropriate fluctuation around the mean value. More importantly, we obtain a velocity distribution for the individual atoms in the central system which follows the expected canonical distribution at the corresponding temperature. This confirms that both our GLE scheme and our mapping procedure onto an extended Langevin dynamics provide the correct thermostat. We also examined the velocity autocorrelation functions and compare our results with more conventional Langevin dynamics.
△ Less
Submitted 18 December, 2014;
originally announced December 2014.
-
The Generalized Langevin Equation: An efficient approach to non-equilibrium molecular dynamics of open systems
Authors:
L. Stella,
C. D. Lorenz,
L. Kantorovich
Abstract:
The Generalized Langevin Equation (GLE) has been recently suggested to simulate the time evolution of classical solid and molecular systems when considering general non-equilibrium processes. In this approach, a part of the whole system (an open system), which interacts and exchanges energy with its dissipative environment, is studied. Because the GLE is derived by projecting out exactly the harmo…
▽ More
The Generalized Langevin Equation (GLE) has been recently suggested to simulate the time evolution of classical solid and molecular systems when considering general non-equilibrium processes. In this approach, a part of the whole system (an open system), which interacts and exchanges energy with its dissipative environment, is studied. Because the GLE is derived by projecting out exactly the harmonic environment, the coupling to it is realistic, while the equations of motion are non-Markovian. Although the GLE formalism has already found promising applications, e.g., in nanotribology and as a powerful thermostat for equilibration in classical molecular dynamics simulations, efficient algorithms to solve the GLE for realistic memory kernels are highly non-trivial, especially if the memory kernels decay non-exponentially. This is due to the fact that one has to generate a colored noise and take account of the memory effects in a consistent manner. In this contribution, we present a simple, yet efficient, algorithm for solving the GLE for practical memory kernels and we demonstrate its capability for the exactly solvable case of a harmonic oscillator coupled to a Debye bath.
△ Less
Submitted 18 March, 2014; v1 submitted 17 December, 2013;
originally announced December 2013.
-
Drift dependence of optimal trade execution strategies under transient price impact
Authors:
Christopher Lorenz,
Alexander Schied
Abstract:
We give a complete solution to the problem of minimizing the expected liquidity costs in presence of a general drift when the underlying market impact model has linear transient price impact with exponential resilience. It turns out that this problem is well-posed only if the drift is absolutely continuous. Optimal strategies often do not exist, and when they do, they depend strongly on the deriva…
▽ More
We give a complete solution to the problem of minimizing the expected liquidity costs in presence of a general drift when the underlying market impact model has linear transient price impact with exponential resilience. It turns out that this problem is well-posed only if the drift is absolutely continuous. Optimal strategies often do not exist, and when they do, they depend strongly on the derivative of the drift. Our approach uses elements from singular stochastic control, even though the problem is essentially non-Markovian due to the transience of price impact and the lack in Markovian structure of the underlying price process. As a corollary, we give a complete solution to the minimization of a certain cost-risk criterion in our setting.
△ Less
Submitted 3 March, 2013; v1 submitted 12 April, 2012;
originally announced April 2012.
-
On the estimation of the curvatures and bending rigidity of membrane networks via a local maximum-entropy approach
Authors:
Fernando Fraternali,
Christan D. Lorenz,
Gianluca Marcelli
Abstract:
We present a meshfree method for the curvature estimation of membrane networks based on the Local Maximum Entropy approach recently presented in (Arroyo and Ortiz, 2006). A continuum regularization of the network is carried out by balancing the maximization of the information entropy corresponding to the nodal data, with the minimization of the total width of the shape functions. The accuracy and…
▽ More
We present a meshfree method for the curvature estimation of membrane networks based on the Local Maximum Entropy approach recently presented in (Arroyo and Ortiz, 2006). A continuum regularization of the network is carried out by balancing the maximization of the information entropy corresponding to the nodal data, with the minimization of the total width of the shape functions. The accuracy and convergence properties of the given curvature prediction procedure are assessed through numerical applications to benchmark problems, which include coarse grained molecular dynamics simulations of the fluctuations of red blood cell membranes (Marcelli et al., 2005; Hale et al., 2009). We also provide an energetic discrete-to-continuum approach to the prediction of the zero-temperature bending rigidity of membrane networks, which is based on the integration of the local curvature estimates. The Local Maximum Entropy approach is easily applicable to the continuum regularization of fluctuating membranes, and the prediction of membrane and bending elasticities of molecular dynamics models.
△ Less
Submitted 16 May, 2011;
originally announced May 2011.
-
Primordial Fluctuations in String Cosmology
Authors:
Larissa C. Lorenz
Abstract:
This thesis is dedicated to the study of inflationary scenarios based on string theory. Inflation is a brief period of accelerated expansion in the very early Universe which gives natural explanations for problems of the Standard Big Bang Model (SBBM) of cosmology. The phenomenological success of inflation provides the motivation for numerous efforts to establish its connection with particle phy…
▽ More
This thesis is dedicated to the study of inflationary scenarios based on string theory. Inflation is a brief period of accelerated expansion in the very early Universe which gives natural explanations for problems of the Standard Big Bang Model (SBBM) of cosmology. The phenomenological success of inflation provides the motivation for numerous efforts to establish its connection with particle physics. As a promising contender for a unified theory at very high energy scales, string theory is an obvious framework to look for the theoretical foundations of inflation; we therefore try to determine the characteristic features of string-inspired inflationary scenarios, along with their observational consequences. In particular, we shall be interested in models of brane inflation, in which the inflaton field has non-canonical dymanics and where the era of reheating is interpreted as the mutual annihilation of (anti-)branes.
This document contains an introduction to cosmology and string theory as well as the model building efforts untertaken at their intersection in recent years. The scientific articles published during the thesis are represented by a short summary each and are available separately on the arXiv.
△ Less
Submitted 10 February, 2010;
originally announced February 2010.
-
Forces between functionalized silica nanoparticles in solution
Authors:
J. Matthew D. Lane,
Ahmed E. Ismail,
Michael Chandross,
Christian D. Lorenz,
Gary S. Grest
Abstract:
To prevent the flocculation and phase separation of nanoparticles in solution, nanoparticles are often functionalized with short chain surfactants. Here we present fully-atomistic molecular dynamics simulations which characterize how these functional coatings affect the interactions between nanoparticles and with the surrounding solvent. For 5 nm diameter silica nanoparticles coated with poly(et…
▽ More
To prevent the flocculation and phase separation of nanoparticles in solution, nanoparticles are often functionalized with short chain surfactants. Here we present fully-atomistic molecular dynamics simulations which characterize how these functional coatings affect the interactions between nanoparticles and with the surrounding solvent. For 5 nm diameter silica nanoparticles coated with poly(ethylene oxide) (PEO) oligomers in water, we determined the hydrodynamic drag on two approaching nanoparticles moving through solvent and on a single nanoparticle as it approaches a planar surface. In most circumstances, acroscale fluid theory accurately predicts the drag on these nano-scale particles. Good agreement is seen with Brenner's analytical solutions for wall separations larger than the soft nanoparticle radius. For two approaching coated nanoparticles, the solvent-mediated (velocity-independent) and lubrication (velocity-dependent) forces are purely repulsive and do not exhibit force oscillations that are typical of uncoated rigid spheres.
△ Less
Submitted 16 April, 2009;
originally announced April 2009.
-
Reheating in a Brane Monodromy Inflation Model
Authors:
Robert H. Brandenberger,
Anke Knauf,
Larissa C. Lorenz
Abstract:
We study reheating in a recently proposed brane "monodromy inflation" model in which the inflaton is the position of a D4 brane on a "twisted torus". Specifically, we study the repeated collisions between the D4 brane and a D6 brane (on which the Standard Model fields are assumed to be localized) at a fixed position along the monodromy direction as the D4 brane rolls down its potential. We find…
▽ More
We study reheating in a recently proposed brane "monodromy inflation" model in which the inflaton is the position of a D4 brane on a "twisted torus". Specifically, we study the repeated collisions between the D4 brane and a D6 brane (on which the Standard Model fields are assumed to be localized) at a fixed position along the monodromy direction as the D4 brane rolls down its potential. We find that there is no trapping of the rolling D4 brane until it reaches the bottom of its potential, and that reheating is entirely described by the last brane encounter. Previous collisions have negligible effect on the brane velocity and hence on the reheat temperature. In the context of our setup, reheating is efficient and the reheat temperature is therefore high.
△ Less
Submitted 17 September, 2008; v1 submitted 28 August, 2008;
originally announced August 2008.
-
Micellar Crystals in Solution from Molecular Dynamics Simulations
Authors:
J. A. Anderson,
C. D. Lorenz,
A. Travesset
Abstract:
Polymers with both soluble and insoluble blocks typically self-assemble into micelles, aggregates of a finite number of polymers where the soluble blocks shield the insoluble ones from contact with the solvent. Upon increasing concentration, these micelles often form gels that exhibit crystalline order in many systems. In this paper, we present a study of both the dynamics and the equilibrium pr…
▽ More
Polymers with both soluble and insoluble blocks typically self-assemble into micelles, aggregates of a finite number of polymers where the soluble blocks shield the insoluble ones from contact with the solvent. Upon increasing concentration, these micelles often form gels that exhibit crystalline order in many systems. In this paper, we present a study of both the dynamics and the equilibrium properties of micellar crystals of triblock polymers using molecular dynamics simulations. Our results show that equilibration of single micelle degrees of freedom and crystal formation occurs by polymer transfer between micelles, a process that is described by transition state theory. Near the disorder (or melting) transition, bcc lattices are favored for all triblocks studied. Lattices with fcc ordering are also found, but only at lower kinetic temperatures and for triblocks with short hydrophilic blocks. Our results lead to a number of theoretical considerations and suggest a range of implications to experimental systems with a particular emphasis on Pluronic polymers.
△ Less
Submitted 18 January, 2008;
originally announced January 2008.
-
Entropy Fluctuations in Brane Inflation Models
Authors:
Robert H. Brandenberger,
Andrew R. Frey,
Larissa C. Lorenz
Abstract:
We study the development of entropy fluctuations in brane inflation in a warped throat, including the brane-antibrane tachyon as the waterfall field. We find that there is a period at the end of inflation during which the entropy mode associated with the tachyon field increases exponentially. In turn, the induced entropy seeds a contribution to the curvature fluctuation on cosmological scales whic…
▽ More
We study the development of entropy fluctuations in brane inflation in a warped throat, including the brane-antibrane tachyon as the waterfall field. We find that there is a period at the end of inflation during which the entropy mode associated with the tachyon field increases exponentially. In turn, the induced entropy seeds a contribution to the curvature fluctuation on cosmological scales which grows rapidly and could exceed the primordial curvature perturbation. We identify parameter values for which in the absence of back-reaction the induced curvature fluctuations are larger than the primordial adiabatic ones. In the specific model we study, however, back-reaction limits the growth of the entropy fluctuations. We discuss situations in which back-reaction effects are less constraining. The lesson of our investigation is that the study of the development of entropy fluctuations at the end of the period of inflation can lead to constraints on models of brane inflation and suggests that the curvaton mechanism may contribute significantly to the spectrum of cosmological perturbations.
△ Less
Submitted 24 June, 2011; v1 submitted 13 December, 2007;
originally announced December 2007.
-
Charge Inversion of Divalent Ionic Solutions in Silica Channels
Authors:
Christian D. Lorenz,
Alex Travesset
Abstract:
Recent experiments (F.H.J. Van Der Heyden et al., PRL 96, 224502 (2006)) of streaming currents in silica nanochannels with divalent ions report charge inversion, i.e. interfacial charges attracting counterions in excess of their own nominal charge, in conflict with existing theoretical and simulation results. We reveal the mechanism of charge inversion by using all-atomic molecular dynamics simu…
▽ More
Recent experiments (F.H.J. Van Der Heyden et al., PRL 96, 224502 (2006)) of streaming currents in silica nanochannels with divalent ions report charge inversion, i.e. interfacial charges attracting counterions in excess of their own nominal charge, in conflict with existing theoretical and simulation results. We reveal the mechanism of charge inversion by using all-atomic molecular dynamics simulations. Our results show excellent agreement with experiments, both qualitatively and quantitatively. We further discuss the implications of our study for the general problem of ionic correlations in solutions as well as in regards of the properties of silica-water interfaces.
△ Less
Submitted 24 June, 2007;
originally announced June 2007.
-
Excess number of percolation clusters on the surface of a sphere
Authors:
Christian D. Lorenz,
Robert M. Ziff
Abstract:
Monte Carlo simulations were performed in order to determine the excess number of clusters b and the average density of clusters n_c for the two-dimensional "Swiss cheese" continuum percolation model on a planar L x L system and on the surface of a sphere. The excess number of clusters for the L x L system was confirmed to be a universal quantity with a value b = 0.8841 as previously predicted a…
▽ More
Monte Carlo simulations were performed in order to determine the excess number of clusters b and the average density of clusters n_c for the two-dimensional "Swiss cheese" continuum percolation model on a planar L x L system and on the surface of a sphere. The excess number of clusters for the L x L system was confirmed to be a universal quantity with a value b = 0.8841 as previously predicted and verified only for lattice percolation. The excess number of clusters on the surface of a sphere was found to have the value b = 1.215(1) for discs with the same coverage as the flat critical system. Finally, the average critical density of clusters was calculated for continuum systems n_c = 0.0408(1).
△ Less
Submitted 1 February, 2001;
originally announced February 2001.
-
Similarity of percolation thresholds on the hcp and fcc lattices
Authors:
Christian D. Lorenz,
Raechelle May,
Robert M. Ziff
Abstract:
Extensive Monte-Carlo simulations were performed in order to determine the precise values of the critical thresholds for site ($p^{hcp}_{c,S} = 0.199 255 5 \pm 0.000 001 0$) and bond ($p^{hcp}_{c,B} = 0.120 164 0 \pm 0.000 001 0$) percolation on the hcp lattice to compare with previous precise measuremens on the fcc lattice. Also, exact enumeration of the hcp and fcc lattices was performed and y…
▽ More
Extensive Monte-Carlo simulations were performed in order to determine the precise values of the critical thresholds for site ($p^{hcp}_{c,S} = 0.199 255 5 \pm 0.000 001 0$) and bond ($p^{hcp}_{c,B} = 0.120 164 0 \pm 0.000 001 0$) percolation on the hcp lattice to compare with previous precise measuremens on the fcc lattice. Also, exact enumeration of the hcp and fcc lattices was performed and yielded generating functions and series for the zeroth, first, and second moments of both lattices. When these series and the values of $p_c$ are compared to those for the fcc lattice, it is apparent that the site percolation thresholds are different; however, the bond percolation thresholds are equal within error bars, and the series only differ slightly in the higher order terms, suggesting the actual values are very close to each other, if not identical.
△ Less
Submitted 2 August, 1999;
originally announced August 1999.
-
Shape-dependent universality in percolation
Authors:
Robert M. Ziff,
Christian D. Lorenz,
Peter Kleban
Abstract:
The shape-dependent universality of the excess percolation cluster number and cross-configuration probability on a torus is discussed. Besides the aspect ratio of the torus, the universality class depends upon the twist in the periodic boundary conditions, which for example are generally introduced when triangular lattices are used in simulations.
The shape-dependent universality of the excess percolation cluster number and cross-configuration probability on a torus is discussed. Besides the aspect ratio of the torus, the universality class depends upon the twist in the periodic boundary conditions, which for example are generally introduced when triangular lattices are used in simulations.
△ Less
Submitted 9 November, 1998;
originally announced November 1998.
-
Universality of the excess number of clusters and the crossing probability function in three-dimensional percolation
Authors:
Christian D. Lorenz,
Robert M. Ziff
Abstract:
Extensive Monte-Carlo simulations were performed to evaluate the excess number of clusters and the crossing probability function for three-dimensional percolation on the simple cubic (s.c.), face-centered cubic (f.c.c.), and body-centered cubic (b.c.c.) lattices. Systems L x L x L' with L' >> L were studied for both bond (s.c., f.c.c., b.c.c.) and site (f.c.c.) percolation. The excess number of…
▽ More
Extensive Monte-Carlo simulations were performed to evaluate the excess number of clusters and the crossing probability function for three-dimensional percolation on the simple cubic (s.c.), face-centered cubic (f.c.c.), and body-centered cubic (b.c.c.) lattices. Systems L x L x L' with L' >> L were studied for both bond (s.c., f.c.c., b.c.c.) and site (f.c.c.) percolation. The excess number of clusters $\tilde {b}$ per unit length was confirmed to be a universal quantity with a value $\tilde {b} \approx 0.412$. Likewise, the critical crossing probability in the L' direction, with periodic boundary conditions in the L x L plane, was found to follow a universal exponential decay as a function of r = L'/L for large r. Simulations were also carried out to find new precise values of the critical thresholds for site percolation on the f.c.c. and b.c.c. lattices, yielding $p_c(f.c.c.)= 0.199 236 5 \pm 0.000 001 0$, $p_c(b.c.c.)= 0.245 961 5\pm 0.000 001 0$.
△ Less
Submitted 29 June, 1998; v1 submitted 18 June, 1998;
originally announced June 1998.
-
Precise determination of the bond percolation thresholds and finite-size scaling corrections for the s.c., f.c.c., and b.c.c. lattices
Authors:
Christian D. Lorenz,
Robert M. Ziff
Abstract:
Extensive Monte-Carlo simulations were performed to study bond percolation on the simple cubic (s.c.), face-centered cubic (f.c.c.), and body-centered cubic (b.c.c.) lattices, using an epidemic kind of approach. These simulations provide very precise values of the critical thresholds for each of the lattices: pc(s.c.) = 0.248 812 6(5), pc(f.c.c.) = 0.120 163 5(10), and pc(b.c.c.) = 0.180 287 5(1…
▽ More
Extensive Monte-Carlo simulations were performed to study bond percolation on the simple cubic (s.c.), face-centered cubic (f.c.c.), and body-centered cubic (b.c.c.) lattices, using an epidemic kind of approach. These simulations provide very precise values of the critical thresholds for each of the lattices: pc(s.c.) = 0.248 812 6(5), pc(f.c.c.) = 0.120 163 5(10), and pc(b.c.c.) = 0.180 287 5(10). For p close to pc, the results follow the expected finite-size and scaling behavior, with values for the Fisher exponent $tau$ (2.189(2)), the finite-size correction exponent $omega$ (0.64(2)), and the scaling function exponent $sigma$ (0.445(1)) confirmed to be universal.
△ Less
Submitted 3 October, 1997;
originally announced October 1997.