CrystalGrower: a generic computer program for Monte Carlo modelling of crystal growth.
Ontology highlight
ABSTRACT:
SUBMITTER: Hill AR
PROVIDER: S-EPMC8179067 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Interacting many-electron problems pose some of the greatest computational challenges in science, with essential applications across many fields. The solutions to these problems will offer accurate predictions of chemical reactivity and kinetics, and other properties of quantum systems1-4. Fermionic quantum Monte Carlo (QMC) methods5,6, which use a statistical sampling of the ground state, are among the most powerful approaches to these problems. Controlling the fermionic sign problem with constraints ensures the efficiency of QMC at the expense of potentially significant biases owing to the limited flexibility of classical computation. Here we propose an approach that combines constrained QMC with quantum computation to reduce such biases. We implement our scheme experimentally using up to 16 qubits to unbias constrained QMC calculations performed on chemical systems with as many as 120 orbitals. These experiments represent the largest chemistry simulations performed with the help of quantum computers, while achieving accuracy that is competitive with state-of-the-art classical methods without burdensome error mitigation. Compared with the popular variational quantum eigensolver7,8, our hybrid quantum-classical computational model offers an alternative path towards achieving a practical quantum advantage for the electronic structure problem without demanding exceedingly accurate preparation and measurement of the ground-state wavefunction.
| S-EPMC8930773 | biostudies-literature
Project description:The invention and advancement of biological microscopy depends critically on an ability to accurately simulate imaging of complex biological structures embedded within complex scattering media. Unfortunately no technique exists for rigorous simulation of the complete imaging process, including the source, instrument, sample and detector. Monte-Carlo modelling is the gold standard for the modelling of light propagation in tissue, but is somewhat laborious to implement and does not incorporate the rejection of scattered light by the microscope. On the other hand microscopes may be rigorously and rapidly modelled using commercial ray-tracing software, but excluding the interaction with the biological sample. We report a hybrid Monte-Carlo optical ray-tracing technique for modelling of complete imaging systems of arbitrary complexity. We make the software available to enable user-friendly and rigorous virtual prototyping of biological microscopy of arbitrary complexity involving light scattering, fluorescence, polarised light propagation, diffraction and coherence. Examples are presented for the modelling and optimisation of representative imaging of neural cells using light-sheet and micro-endoscopic fluorescence microscopy and imaging of retinal vasculature using confocal and non-confocal scanning-laser ophthalmoscopes.
| S-EPMC6825179 | biostudies-literature
Project description:MotivationInvestigating the 3D structure of chromatin provides new insights into transcriptional regulation. With the evolution of 3C next-generation sequencing methods like ChiA-PET and Hi-C, the surge in data volume has highlighted the need for more efficient chromatin spatial modelling algorithms. This study introduces the cudaMMC method, based on the Simulated Annealing Monte Carlo approach and enhanced by GPU-accelerated computing, to efficiently generate ensembles of chromatin 3D structures.ResultsThe cudaMMC calculations demonstrate significantly faster performance with better stability compared to our previous method on the same workstation. cudaMMC also substantially reduces the computation time required for generating ensembles of large chromatin models, making it an invaluable tool for studying chromatin spatial conformation.Availability and implementationOpen-source software and manual and sample data are freely available on https://github.com/SFGLab/cudaMMC.
| S-EPMC10568367 | biostudies-literature
Project description:PurposeThe purpose of this work is to calculate individualized dose distributions in patients undergoing 18 F-FDG PET/CT studies through a methodology based on full Monte Carlo (MC) simulations and PET/CT patient images, and to compare such values with those obtained by employing nonindividualized phantom-based methods.MethodsWe developed a MC-based methodology for individualized internal dose calculations, which relies on CT images (for organ segmentation and dose deposition), PET images (for organ segmentation and distributions of activities), and a biokinetic model (which works with information provided by PET and CT images) to obtain cumulated activities. The software vGATE version 8.1. was employed to carry out the Monte Carlo calculations. We also calculated deposited doses with nonindividualized phantom-based methods (Cristy-Eckerman, Stabin, and ICRP-133).ResultsMedian MC-calculated dose/activity values are within 0.01-0.03 mGy/MBq for most organs, with higher doses delivered especially to the bladder wall, major vessels, and brain (medians of 0.058, 0.060, 0.066 mGy/MBq, respectively). Comparison with values obtained with nonindividualized phantom-based methods has shown important differences in many cases (ranging from -80% to + 260%). These differences are significant (p < 0.05) for several organs/tissues, namely, remaining tissues, adrenals, bladder wall, bones, upper large intestine, heart, pancreas, skin, and stomach wall.ConclusionsThe methodology presented in this work is a viable and useful method to calculate internal dose distributions in patients undergoing medical procedures involving radiopharmaceuticals, individually, with higher accuracy than phantom-based methods, fulfilling the guidelines provided by the European Council directive 2013/59/Euratom.
| S-EPMC7586975 | biostudies-literature
Project description:The radionuclide copper-64 is a promising candidate for nuclear medicine, but its complex decay creates challenges in the primary standardization of its activity. Monte Carlo simulations of live-timed anticoincidence (LTAC) counting of 64Cu were used to calculate corrections to extrapolation intercepts, resulting in improved activity determinations. A small correction (-0.33%) to the linear extrapolation of LTAC data acquired with a ?-gate over the 1346keV gamma peak was determined. We discuss the physical origin of the correction. We also use experimental data to demonstrate a Monte Carlo scaling that allows for inclusion of data acquired with a ?-gate set over the annihilation photon peak(s).
| S-EPMC6237201 | biostudies-literature
Project description:The rate at which a given site in a gene sequence alignment evolves over time may vary. This phenomenon--known as heterotachy--can bias or distort phylogenetic trees inferred from models of sequence evolution that assume rates of evolution are constant. Here, we describe a phylogenetic mixture model designed to accommodate heterotachy. The method sums the likelihood of the data at each site over more than one set of branch lengths on the same tree topology. A branch-length set that is best for one site may differ from the branch-length set that is best for some other site, thereby allowing different sites to have different rates of change throughout the tree. Because rate variation may not be present in all branches, we use a reversible-jump Markov chain Monte Carlo algorithm to identify those branches in which reliable amounts of heterotachy occur. We implement the method in combination with our 'pattern-heterogeneity' mixture model, applying it to simulated data and five published datasets. We find that complex evolutionary signals of heterotachy are routinely present over and above variation in the rate or pattern of evolution across sites, that the reversible-jump method requires far fewer parameters than conventional mixture models to describe it, and serves to identify the regions of the tree in which heterotachy is most pronounced. The reversible-jump procedure also removes the need for a posteriori tests of 'significance' such as the Akaike or Bayesian information criterion tests, or Bayes factors. Heterotachy has important consequences for the correct reconstruction of phylogenies as well as for tests of hypotheses that rely on accurate branch-length information. These include molecular clocks, analyses of tempo and mode of evolution, comparative studies and ancestral state reconstruction. The model is available from the authors' website, and can be used for the analysis of both nucleotide and morphological data.
| S-EPMC2607421 | biostudies-literature
Project description:In the three years since SARS-CoV-2 was first detected in China, hundreds of millions of people have been infected and millions have died. Along with the immediate need for treatment solutions, the COVID-19 epidemic has reinforced the need for mathematical models that can predict the spread of the pandemic in an ever-changing environment. The susceptible-infectious-removed (SIR) model has been widely used to model COVID-19 transmission, however, with limited success. Here, we present a novel, dynamic Monte-Carlo Agent-based Model (MAM), which is based on the basic principles of statistical physics. Using public aggregative data from Israel on three major outbreaks, we compare predictions made by SIR and MAM, and show that MAM outperforms SIR in all aspects. Furthermore, MAM is a flexible model and allows to accurately examine the effects of vaccinations in different subgroups, and the effects of the introduction of new variants.
| S-EPMC10098313 | biostudies-literature
Project description:Modelling of non-adiabatic dynamics in extended molecular systems and solids is a next frontier of atomistic electronic structure theory. The underlying numerical algorithms should operate only with a few quantities (that can be efficiently obtained from quantum chemistry), provide a controlled approximation (which can be systematically improved) and capture important phenomena such as branching (multiple products), detailed balance and evolution of electronic coherences. Here we propose a new algorithm based on Monte-Carlo sampling of classical trajectories, which satisfies the above requirements and provides a general framework for existing surface hopping methods for non-adiabatic dynamics simulations. In particular, our algorithm can be viewed as a post-processing technique for analysing numerical results obtained from the conventional surface hopping approaches. Presented numerical tests for several model problems demonstrate efficiency and accuracy of the new method.
| S-EPMC3759042 | biostudies-literature
Project description:Molecular interactions across intercellular interfaces serve to convey information between cells and to trigger appropriate cell functions. Examples include cell development and growth in tissues, neuronal and immune synapses (ISs). Here, we introduce an agent-based Monte-Carlo simulation of user-defined cellular interfaces. The simulation allows for membrane molecules, embedded at intercellular contacts, to diffuse and interact, while capturing the topography and energetics of the plasma membranes of the interface. We provide a detailed example related to pattern formation in the early IS. Using simulation predictions and three-color single molecule localization microscopy (SMLM), we detected the intricate mutual patterning of T cell antigen receptors (TCRs), integrins and glycoproteins in early T cell contacts with stimulating coverslips. The simulation further captures the dynamics of the patterning under the experimental conditions and at the IS with antigen presenting cells (APCs). Thus, we provide a generic tool for simulating realistic cell-cell interfaces, which can be used for critical hypothesis testing and experimental design in an iterative manner.
| S-EPMC6141710 | biostudies-literature
Project description:The coronavirus disease 2019 (COVID-19) has now spread throughout most countries in the world causing heavy life losses and damaging social-economic impacts. Following a stochastic point process modelling approach, a Monte Carlo simulation model was developed to represent the COVID-19 spread dynamics. First, we examined various expected performances (theoretical properties) of the simulation model assuming a number of arbitrarily defined scenarios. Simulation studies were then performed on the real COVID-19 data reported (over the period of 1 March to 1 May) for Australia and United Kingdom (UK). Given the initial number of COVID-19 infection active cases were around 10 for both countries, the model estimated that the number of active cases would peak around 29 March in Australia (≈ 1,700 cases) and around 22 April in UK (≈ 22,860 cases); ultimately the total confirmed cases could sum to 6,790 for Australia in about 75 days and 206,480 for UK in about 105 days. The results of the estimated COVID-19 reproduction numbers were consistent with what was reported in the literature. This simulation model was considered an effective and adaptable decision making/what-if analysis tool in battling COVID-19 in the immediate need, and for modelling any other infectious diseases in the future.
| S-EPMC7403316 | biostudies-literature