Gold catalysis for organic synthesis II
F. Toste
DOI: https://doi.org/10.3762/bjoc.9.241
2013-10-09
Beilstein Journal of Organic Chemistry
Abstract:A combined experimental and computational study on regioselective gold-catalyzed synthetic routes to 1,3-oxazinan-2-ones (kinetically controlled products) and 1,3-oxazin-2-one derivatives (thermodynamically favored) from easily accessible allenic carbamates has been carried out. Introduction The search for new synthetic routes to 1,3-oxazin-2-one derivatives [1] is of interest because of the biological activity of these molecules [2-7]. Carretero and colleagues have published the Au(I)-catalyzed cyclization of N-Boc-3-butyn-1-amines to afford six-membered 2-oxazinones, namely, 6-methylene-1,3oxazinan-2-ones involving a 6-exo-dig cyclization [1]. This interesting method avoids complicated prefunctionalization of starting materials and minimizes the formation of byproducts; however, it is used just for the preparation of four simple examples. 1,3-Oxazin-2-ones are also used as valuable intermediates in organic synthesis [8-14]. Recently, allenes have attracted much attention as they have been used for the preparation of both biologically relevant drugs as well as advanced materials [15-27]. However, regioselectivity problems are significant (endo-trig versus endo-dig versus exo-dig versus exo-trig cyclization). Great effort is currently being made in the search for a variety of reactions promoted by gold salts due to their impressive catalytic properties [28-41]. The Boc protective group has been widely used in allene chemistry, being an inert and recommended partner in goldand palladium-catalyzed Beilstein J. Org. Chem. 2013, 9, 818–826. 819 aminocyclizations of allenes [42]. On the other hand, reports of gold-catalyzed cyclizations leading to heterocycles that contain more than one heteroatom are rare [43-48]. Besides, it has been reported very recently that the Au(I)-catalyzed cyclization of a N-phenethyl-N-Boc-protected allenamide failed [49]. Despite the above precedents, but in continuation of our interest in heterocyclic and allene chemistry [50-55], we decided to examine the gold-catalyzed cyclization of N-Boc-allenes with the aim of establishing a protocol for the synthesis of 1,3oxazin-2-one derivatives in which the carbamate group should serve as the source of CO2. Results and Discussion To explore the effects of various substrates on gold-catalyzed oxycyclization reactions, a number of new allenic carbamates were synthesized as shown in Scheme 1. Starting materials, tert-butyl (prop-2-ynyl)carbamates 1a–j, were obtained both in the racemic form and in optically pure form by using standard methodologies. Thus, alkynylcarbamates 1a–g were prepared through reductive amination of the appropriate aldehyde with propargylamine, followed by Boc2O treatment of the corresponding N-substituted prop-2-yn-1-amine. Alkynylcarbamate 1h was prepared from Garner’s aldehyde following a literature report [56,57]. Alkynylcarbamate 1i was readily accessed from (S)-prolinol by using a modified known procedure [58]. Alkynylcarbamate 1j was achieved through the reaction of 3-bromo-1H-indole-2-carbaldehyde with the Ohira–Bestmann reagent followed by the addition of Boc2O. Terminal alkynes 1 were conveniently converted into allenic carbamates 2 by treatment with paraformaldehyde in the presence of diisopropylamine and copper(I) bromide (Crabbé reaction) [59,60]. We employed three different gold salts in our initial screening of catalysts for the model system, allenic carbamate 2a. Initially, the use of AuCl3 and AuCl were tested, but both failed to catalyze the reaction. Fortunately, we found that [AuClPPh3]/ AgOTf was an excellent catalyst for our purpose. To our delight, the reaction of allenic carbamate 2a at room temperature afforded 3-benzyl-6-methylene-1,3-oxazinan-2-one (3a) bearing an exocyclic double bond as the sole product (Scheme 2). Adding a catalytic amount of Brønsted acid (PTSA) into the reaction system did slightly improve the yield of 3a. Solvent screening demonstrated that dichloromethane was the best choice in the reaction. As revealed in Scheme 2, a variety of allenic carbamates 2 were also suitable for such heterocyclization reactions to afford 1,3oxazinan-2-ones 3. To increase the molecular diversity by incorporating more 1,3-oxazin-2-ones in the molecule, compound 2g having two allenic carbamate units was used. Notably, bis(allenic carbamate) 2g also undergoes this interesting transformation to give bis(6-methylene-1,3-oxazinan-2-one) 3g through a two-fold cyclization. This product particularly underlines the power of the present cyclization reaction, as none of the conventional methods would allow its synthesis with such great ease. Interestingly, as a first try, we were pleased to notice that the reaction of allenyl derivative 2a in dichloromethane at 90 °C, afforded 1,3-oxazin-2-one 4a bearing an endocyclic double bond as the major component, and 1,3-oxazinan-2-one 3a was also isolated as a minor component. Notably, starting from allenic carbamates 2a–j and performing the reaction in dichloromethane at 130 °C, a series of 6-methyl-3-substituted 3,4-dihydro-2H-1,3-oxazin-2-ones 4a–j were exclusively formed (Scheme 2) [61-68]. The observed regioselectivity is worthy of note, because under our reaction conditions only 1,3oxazinan-2-ones 3 (arising from 6-endo-dig cyclization) or 3,4dihydro-2H-1,3-oxazin-2-ones 4 (arising from 6-exo-dig cyclization) were achieved, with the nucleophilic oxygen attacking the central allene carbon atom in each case. This is an interesting result, because the available examples on related metal-catalyzed allene heterocyclizations usually lead to 5-exotrig cyclization [69,70]; only Hashmi et al. have recently reported an attack at the central position of the allene in allenylamides [44]. Thus, it is possible to suppress the formation of the 1,3-oxazinan-2-one ring by performing the reaction at higher temperature, yielding the 1,3-oxazin-2-one as the exclusive product. A general trend can be deduced on the basis of these results: heterocycle 4 is the thermodynamically controlled product while heterocycle 3 is the kinetically controlled product [7173]. Probably, double-bond migration in compounds 3 results in the formation of the 1,3-oxazin-2-one 4. In order to verify the role of the Au(I) catalyst in the double-bond migration process, we set up two experiments. Heating a mixture of 3a with Au(OTf)PPh3 at a loading of 2.5 mol % in dichloromethane for 1.5 h at 130 °C resulted in full conversion into 4a. Running the same reaction in the absence of any catalyst resulted in 30% conversion after two days, as determined by 1H NMR. Treatment of 1,3-oxazinan-2-one 3a with 5 mol % TfOH in CH2Cl2 at room temperature did not proceed to give an appreciable amount of 3,4-dihydro-2H-1,3-oxazin-2-one 4a after 2 h. This indicates that the Au(I) catalyst might participate in the doublebond migration process; being a possible intermediate, the π-allyl complex 5 is depicted in Scheme 3 [74]. Despite that, the isomerization process can be also viewed as an intramolecular 1,3-H shift assisted by gold. A possible pathway for the gold-catalyzed achievement of heterocycles 3 from allenyl-tethered carbamates 2 may initially Beilstein J. Org. Chem. 2013, 9, 818–826. 820 Scheme 1: Preparation of allenic carbamates 2a–j. Reagents and conditions: (i) Propargylamine, MgSO4, CH2Cl2, rt, 15 h. (ii) NaBH4, MeOH, rt, 0.5 h. (iii) Boc2O, Et3N, CH2Cl2, rt, 2–15 h. (iv) (CH2O)n, iPr2NH, CuBr, 1,4-dioxane, reflux, 1 h. (v) Ohira–Bestmann reagent, K2CO3, MeOH, rt, 15 h. (vi) Boc2O, DMAP, CH3CN, rt, 2 h. (vii) Dess–Martin periodinane, CH2Cl2, rt. PMP = 4-MeOC6H4. Beilstein J. Org. Chem. 2013, 9, 818–826. 821 Scheme 2: Controlled oxycyclization reactions of allenic carbamates 2 to 1,3-oxazinan-2-ones 3 and 1,3-oxazin-2-ones 4 under selective goldcatalyzed conditions. Reagents and conditions: (i) 2.5 mol % [AuClPPh3], 2.5 mol % AgOTf, 10 mol % PTSA, CH2Cl2, rt, 3a: 6 h; 3b: 7 h, 3d: 6 h; 3e: 8 h; 3f: 7 h; 3g: 5.5 h; 3h: 5 h; 3i: 7 h; 3j: 2 h. (ii) 2.5 mol % [AuClPPh3], 2.5 mol % AgOTf, 10 mol % PTSA, CH2Cl2, sealed tube, 130 °C, 4a: 1.5 h; 4b: 2 h, 4c: 4 h; 4d: 6 h; 4e: 4 h; 4f: 0.5 h; 4g: 5.5 h; 4i: 2 h. (iii) 2.5 mol % [AuClPPh3], 2.5 mol % AgOTf, 10 mol % PTSA, CH2Cl2, sealed tube, 80 °C, 4j: 2 h. PMP = 4-MeOC6H4. Beilstein J. Org. Chem. 2013, 9, 818–826. 822 Scheme 3: Possible explanation for the gold-catalyzed isomerization reaction of 1,3-oxazinan-2-one 3a into 1,3-oxazin-2-one 4a. Reagents and conditions: (i) 2.5 mol % [AuClPPh3], 2.5 mol % AgOTf, 10 mol % PTSA, CH2Cl2, sealed tube, 130 °C, 1.5 h. Scheme 4: Mechanistic explanation for the gold-catalyzed oxycyclization reactions of allenic carbamates 2 into 6-methylene-1,3-oxazinan-2-ones 3 or into 3-substituted-3,4-dihydro-2H-1,3-oxazin-2-ones 4. involve the formation of a complex 6 through coordination of the gold salt to the proximal allenic double bond. Next, chemoand regioselective 6-endo-dig oxyauration of the carbamate carbonyl moiety forms species 7. Attack of the carbamate carbonyl group occurs as a result of the stability of the intermediate ammonium cation type 7. Loss of proton linked to 2-methylprop-1-ene release [75-78], generates neutral species 8, which followed by protonolysis of the carbon–gold bond affords 6-methylene-1,3-oxazinan-2-ones 3 with concurrent regeneration of the gold catalyst (Scheme 4, left catalytic cycle). In line with the above mechanistic proposal, the easy breakage of the tert-butyl group at species 7 is essential for the formation of 1,3-oxazinan-2-ones 3. Besides, the replacement of the tertbutyl group in allenic carbamates 2 by other alkyl functions, such as methyl, did not allow the preparation of heterocycles 3. In addition to the double-bond isomerization that transforms products 3 into the thermodynamically more favored compounds 4, a mechanistic scenario involving the initial coordination of the gold to the distal allenic double bond leading to complex 9, followed by a 6-exo-dig oxyauration is likely for the achieve
What problem does this paper attempt to address?
-
Thermodynamic calculations using reverse Monte Carlo: Simultaneously tuning multiple short-range order parameters for 2D lattice adsorption problem
Suhail Haque,Abhijit Chatterjee
2023-07-22
Abstract:Lattice simulations are an important class of problems in crystalline solids, surface science, alloys, adsorption, absorption, separation, catalysis, to name a few. We describe a fast computational method for performing lattice thermodynamic calculations that is based on the use of the reverse Monte Carlo (RMC) technique and multiple short-range order (SRO) parameters. The approach is comparable in accuracy to the Metropolis Monte Carlo (MC) method. The equilibrium configuration is determined in 5-10 Newton-Raphson iterations by solving a system of coupled nonlinear algebraic flux equations. This makes the RMC-based method computationally more efficient than MC, given that MC typically requires sampling of millions of configurations. The technique is applied to the interacting 2D adsorption problem. Unlike grand canonical MC, RMC is found to be adept at tackling geometric frustration, as it is able to quickly and correctly provide the ordered c(2x2) adlayer configuration for Cl adsorbed on a Cu (100) surface.
Statistical Mechanics
-
Self-driven lattice-model Monte Carlo simulations of alloy thermodynamic
A. van de Walle,M. Asta
DOI: https://doi.org/10.1088/0965-0393/10/5/304
2002-01-29
Abstract:Monte Carlo (MC) simulations of lattice models are a widely used way to compute thermodynamic properties of substitutional alloys. A limitation to their more widespread use is the difficulty of driving a MC simulation in order to obtain the desired quantities. To address this problem, we have devised a variety of high-level algorithms that serve as an interface between the user and a traditional MC code. The user specifies the goals sought in a high-level form that our algorithms convert into elementary tasks to be performed by a standard MC code. For instance, our algorithms permit the determination of the free energy of an alloy phase over its entire region of stability within a specified accuracy, without requiring any user intervention during the calculations. Our algorithms also enable the direct determination of composition-temperature phase boundaries without requiring the calculation of the whole free energy surface of the alloy system.
Statistical Mechanics
-
Nonequilibrium candidate Monte Carlo: A new tool for efficient equilibrium simulation
Jerome P. Nilmeier,Gavin E. Crooks,David D. L. Minh,John D. Chodera
DOI: https://doi.org/10.1073/pnas.1106094108
2011-10-21
Abstract:Metropolis Monte Carlo simulation is a powerful tool for studying the equilibrium properties of matter. In complex condensed-phase systems, however, it is difficult to design Monte Carlo moves with high acceptance probabilities that also rapidly sample uncorrelated configurations. Here, we introduce a new class of moves based on nonequilibrium dynamics: candidate configurations are generated through a finite-time process in which a system is actively driven out of equilibrium, and accepted with criteria that preserve the equilibrium distribution. The acceptance rule is similar to the Metropolis acceptance probability, but related to the nonequilibrium work rather than the instantaneous energy difference. Our method is applicable to sampling from both a single thermodynamic state or a mixture of thermodynamic states, and allows both coordinates and thermodynamic parameters to be driven in nonequilibrium proposals. While generating finite-time switching trajectories incurs an additional cost, driving some degrees of freedom while allowing others to evolve naturally can lead to large enhancements in acceptance probabilities, greatly reducing structural correlation times. Using nonequilibrium driven processes vastly expands the repertoire of useful Monte Carlo proposals in simulations of dense solvated systems.
Statistical Mechanics
-
Off-Lattice Kinetic Monte Carlo Methods
Mickaël Trochet,Normand Mousseau,Laurent Karim Béland,Graeme Henkelman
DOI: https://doi.org/10.1007/978-3-319-42913-7_29-1
2018-01-01
Abstract:Exact modeling of the dynamics of chemical and material systems over experimentally relevant time scales still eludes us even with modern computational resources. Fortunately, many systems can be described as rare event systems where atoms vibrate around equilibrium positions for a long time before a transition is made to a new atomic state. For those systems, the kinetic Monte Carlo (KMC) algorithm provides a powerful solution. In traditional KMC, mechanism and rates are computed beforehand, limiting moves to discretized positions and largely ignoring strain. Many systems of interest, however, are not well-represented by such lattice-based models. Moreover, materials often evolve with complex and concerted mechanisms that cannot be anticipated before the start of a simulation. In this chapter, we describe a class of algorithms, called off-lattice or adaptive KMC, which relaxes both limitations of traditional KMC, with atomic configurations represented in the full configuration space and reaction events are calculated on-the-fly, with the possible use of catalogs to speed up calculations. We discuss a number of implementations of off-lattice KMC developed by different research groups, emphasizing the similarities between the approaches that open modeling to new classes of problems.
-
Accelerating the discovery of low-energy structure configurations: a computational approach that integrates first-principles calculations, Monte Carlo sampling, and Machine Learning
Md Rajib Khan Musa,Yichen Qian,Jie Peng,David Cereceda
2024-10-08
Abstract:Finding Minimum Energy Configurations (MECs) is essential in fields such as physics, chemistry, and materials science, as they represent the most stable states of the systems. In particular, identifying such MECs in multi-component alloys considered candidate PFMs is key because it determines the most stable arrangement of atoms within the alloy, directly influencing its phase stability, structural integrity, and thermo-mechanical properties. However, since the search space grows exponentially with the number of atoms considered, obtaining such MECs using computationally expensive first-principles DFT calculations often results in a cumbersome task. To escape the above compromise between physical fidelity and computational efficiency, we have developed a novel physics-based data-driven approach that combines Monte Carlo sampling, first-principles DFT calculations, and Machine Learning to accelerate the discovery of MECs in multi-component alloys. More specifically, we have leveraged well-established Cluster Expansion (CE) techniques with Local Outlier Factor models to establish strategies that enhance the reliability of the CE method. In this work, we demonstrated the capabilities of the proposed approach for the particular case of a tungsten-based quaternary high-entropy alloy. However, the method is applicable to other types of alloys and enables a wide range of applications.
Materials Science,Machine Learning
-
Multilevel Monte Carlo for quantum mechanics on a lattice
Karl Jansen,Eike Hermann Müller,Robert Scheichl
DOI: https://doi.org/10.1103/PhysRevD.102.114512
2020-11-17
Abstract:Monte Carlo simulations of quantum field theories on a lattice become increasingly expensive as the continuum limit is approached since the cost per independent sample grows with a high power of the inverse lattice spacing. Simulations on fine lattices suffer from critical slowdown, the rapid growth of autocorrelations in the Markov chain. This causes a strong increase in the number of lattice configurations that have to be generated to obtain statistically significant results. This paper discusses hierarchical sampling methods to tame the growth in autocorrelations. Combined with multilevel variance reduction, this significantly reduces the computational cost of simulations for given tolerances $\epsilon_{\text{disc}}$ on the discretisation error and $\epsilon_{\text{stat}}$ on the statistical error. For observables with lattice errors of order $\alpha$ and integrated autocorrelation times that grow like $\tau_{\mathrm{int}}\propto a^{-z}$, multilevel Monte Carlo (MLMC) reduces the cost from $\mathcal{O}(\epsilon_{\text{stat}}^{-2}\epsilon_{\text{disc}}^{-(1+z)/\alpha})$ to $\mathcal{O}(\epsilon_{\text{stat}}^{-2}\vert\log \epsilon_{\text{disc}} \vert^2+\epsilon_{\text{disc}}^{-1/\alpha})$ or $\mathcal{O}(\epsilon_{\text{stat}}^{-2}+\epsilon_{\text{disc}}^{-1/\alpha})$. Higher gains are expected for simulations of quantum field theories in $D$ dimensions. The efficiency of the approach is demonstrated on two model systems, including a topological oscillator that is badly affected by critical slowdown from topological charge freezing. On fine lattices, the new methods are orders of magnitude faster than standard Hybrid Monte Carlo sampling. For high resolutions, MLMC can be used to accelerate even the cluster algorithm for the topological oscillator. Performance is further improved through perturbative matching which guarantees efficient coupling of theories on the multilevel hierarchy.
High Energy Physics - Lattice,Numerical Analysis,Computational Physics
-
Adaptable Parallel Acceleration Strategy For Dynamic Monte Carlo Simulations Of Polymerization With Microscopic Resolution
Rui Liu,Antonios Armaou,Xi Chen
DOI: https://doi.org/10.1021/acs.iecr.0c05795
2021-01-01
Industrial & Engineering Chemistry Research
Abstract:Properties of polymer products are determined by their microscopic structures. Dynamic Monte Carlo (DMC) simulation is a powerful tool to capture detailed polymer microstructures. However, the heavy computational burden significantly limits the broad application of DMC in practice. In this work, a sub-box parallel strategy is proposed to accelerate the DMC simulation of polymerization processes with accuracy assurance. Different parallelizable task granularities are derived for different kinetic reaction mechanisms. The integration of reaction mechanism and hardware architecture is fully addressed by proposing the implementation of the simulations on a multicore processor platform or a graphics processing unit platform according to parallelizable task granularity. Key parallelization questions including random number generation, data structure, and communication strategy are purposely answered for different scenarios. Five case studies with kinetic mechanisms of increasing complexity, including linear and branching polymerizations, are presented to show the efficiency, accuracy, and reliability of the proposed parallelization strategy.
-
Calculations of excess free energies of precipitates via direct thermodynamic integration across phase boundaries
Babak Sadigh,Paul Erhart
DOI: https://doi.org/10.1103/PhysRevB.86.134204
2012-09-05
Abstract:We describe a technique for constraining macroscopic fluctuations in thermodynamic variables well-suited for Monte Carlo (MC) simulations of multiphase equilibria. In particular for multicomponent systems this amounts to a statistical ensemble that implements constraints on both the average composition as well as its fluctuations. The variance-constrained semi-grandcanonical (VC-SGC) ensemble allows for MC simulations, in which single-phase systems can be reversibly switched into multiphase equilibria allowing the calculation of excess free energies of precipitates of complex shapes by thermodynamic integration. The basic features as well as the scaling and convergence properties of this technique are demonstrated by application to an Ising model. Finally, the VC-SGC MC simulation technique is used to calculate alpha/alpha' interface free energies in Fe-Cr alloys as a function of orientation and temperature taking into account configurational, vibrational, and structural degrees of freedom.
Materials Science,Statistical Mechanics
-
A new algorithm for Monte Carlo simulation of Ising spin systems
A.B. Bortz,M.H. Kalos,J.L. Lebowitz
DOI: https://doi.org/10.1016/0021-9991(75)90060-1
IF: 4.645
1975-01-01
Journal of Computational Physics
Abstract:We describe a new algorithm for Monte Carlo simulation of Ising spin systems and present results of a study comparing the speed of the new technique to that of a standard technique applied to a square lattice of 6400 spins evolving via single spin flips. We find that at temperatures T < Tc, the critical temperature, the new technique is faster than the standard technique, being ten times faster at T = 0.588 Tc. We expect that the new technique will be especially valuable in Monte Carlo simulation of the time evolution of binary alloy systems. The new algorithm is essentially a reorganization of the standard algorithm. It accounts for the a priori probability of changing spins before, rather than after, choosing the spin or spins to change.
physics, mathematical,computer science, interdisciplinary applications
-
Fast and scalable quantum Monte Carlo simulations of electron-phonon models
Benjamin Cohen-Stead,Owen Bradley,Cole Miles,George Batrouni,Richard Scalettar,Kipton Barros
DOI: https://doi.org/10.1103/PhysRevE.105.065302
2022-07-15
Abstract:We introduce methodologies for highly scalable quantum Monte Carlo simulations of electron-phonon models, and report benchmark results for the Holstein model on the square lattice. The determinant quantum Monte Carlo (DQMC) method is a widely used tool for simulating simple electron-phonon models at finite temperatures, but incurs a computational cost that scales cubically with system size. Alternatively, near-linear scaling with system size can be achieved with the hybrid Monte Carlo (HMC) method and an integral representation of the Fermion determinant. Here, we introduce a collection of methodologies that make such simulations even faster. To combat "stiffness" arising from the bosonic action, we review how Fourier acceleration can be combined with time-step splitting. To overcome phonon sampling barriers associated with strongly-bound bipolaron formation, we design global Monte Carlo updates that approximately respect particle-hole symmetry. To accelerate the iterative linear solver, we introduce a preconditioner that becomes exact in the adiabatic limit of infinite atomic mass. Finally, we demonstrate how stochastic measurements can be accelerated using fast Fourier transforms. These methods are all complementary and, combined, may produce multiple orders of magnitude speedup, depending on model details.
Strongly Correlated Electrons,Statistical Mechanics,Computational Physics
-
Exact distributed kinetic Monte Carlo simulations for on-lattice chemical kinetics: lessons learnt from medium- and large-scale benchmarks
Giannis D Savva,Raz L Benson,Ilektra-Athanasia Christidi,Michail Stamatakis
DOI: https://doi.org/10.1098/rsta.2022.0235
2023-07-10
Abstract:Kinetic Monte Carlo (KMC) simulations have been instrumental in multiscale catalysis studies, enabling the elucidation of the complex dynamics of heterogeneous catalysts and the prediction of macroscopic performance metrics, such as activity and selectivity. However, the accessible length- and time-scales have been a limiting factor in such simulations. For instance, handling lattices containing millions of sites with 'traditional' sequential KMC implementations is prohibitive owing to large memory requirements and long simulation times. We have recently established an approach for exact, distributed, lattice-based simulations of catalytic kinetics which couples the Time-Warp algorithm with the Graph-Theoretical KMC framework, enabling the handling of complex adsorbate lateral interactions and reaction events within large lattices. In this work, we develop a lattice-based variant of the Brusselator system, a prototype chemical oscillator pioneered by Prigogine and Lefever in the late 60s, to benchmark and demonstrate our approach. This system can form spiral wave patterns, which would be computationally intractable with sequential KMC, while our distributed KMC approach can simulate such patterns 15 and 36 times faster with 625 and 1600 processors, respectively. The medium- and large-scale benchmarks thus conducted, demonstrate the robustness of the approach, and reveal computational bottlenecks that could be targeted in further development efforts. This article is part of a discussion meeting issue 'Supercomputing simulations of advanced materials'.
-
CASM Monte Carlo: Calculations of the thermodynamic and kinetic properties of complex multicomponent crystals
Brian Puchala,John C. Thomas,Anton Van der Ven
2023-09-21
Abstract:Monte Carlo techniques play a central role in statistical mechanics approaches for connecting macroscopic thermodynamic and kinetic properties to the electronic structure of a material. This paper describes the implementation of Monte Carlo techniques for the study multicomponent crystalline materials within the Clusters Approach to Statistical Mechanics (CASM) software suite, and demonstrates their use in model systems to calculate free energies and kinetic coefficients, study phase transitions, and construct first-principles based phase diagrams. Many crystal structures are complex, with multiple sublattices occupied by differing sets of chemical species, along with the presence of vacancies or interstitial species. This imposes constraints on concentration variables, the form of thermodynamic potentials, and the values of kinetic transport coefficients. The framework used by CASM to formulate thermodynamic potentials and kinetic transport coefficients accounting for arbitrarily complex crystal structures is presented and demonstrated with examples applying it to crystal systems of increasing complexity. Additionally, a new software package is introduced, casm-flow, which helps automate the setup, submission, management, and analysis of Monte Carlo simulations performed using CASM.
Materials Science
-
Micromagnetic Monte Carlo method with variable magnetization length based on the Landau-Lifshitz-Bloch equation for computation of large-scale thermodynamic equilibrium states
Serban Lepadatu
DOI: https://doi.org/10.1063/5.0059745
2021-09-14
Abstract:An efficient method for computing thermodynamic equilibrium states at the micromagnetic length scale is introduced, using the Markov chain Monte Carlo method. Trial moves include not only rotations of vectors, but also a change in their magnetization length. The method is parameterized using the longitudinal susceptibility, reproduces the same Maxwell-Boltzmann distribution as the stochastic Landau-Lifshitz-Bloch equation, and is applicable both below and above the Curie temperature. The algorithm is fully parallel, can be executed on graphical processing units, and efficiently includes the long range dipolar interaction. This method is generally useful for computing finite-temperature relaxation states both for uniform and non-uniform temperature profiles, and can be considered as complementary to zero-temperature micromagnetic energy minimization solvers, with comparable computation time. Compared to the dynamic approach it is shown the micromagnetic Monte Carlo method is up to almost 20 times faster. Moreover, unlike quasi-zero temperature approaches which do not take into account the magnetization length distribution and stochasticity, the method is better suited for structures with unbroken symmetry around the applied field axis, granular films, and at higher temperatures and fields. In particular, applications to finite-temperature hysteresis loop modelling, chiral magnetic thin films, granular magnetic media, and artificial spin ices are discussed.
Mesoscale and Nanoscale Physics
-
Highly efficient parallel grand canonical simulations of interstitial-driven deformation-diffusion processes
Seyed Sina Moeini-Ardakani,Seyedeh Mohadeseh Taheri-Mousavi,Ju Li
DOI: https://doi.org/10.1088/1361-651X/ac01b9
2021-02-14
Abstract:Diffusion of interstitial alloying elements like H, O, C, and N in metals and their continuous relocation and interactions with their microstructures have crucial influences on metals properties. However, besides limitations in experimental tools in capturing these mechanisms, the inefficiency of numerical tools also inhibits modeling efforts. Here, we present an efficient framework to perform hybrid grand canonical Monte Carlo and molecular dynamics simulations that allow for parallel insertion/deletion of Monte Carlo moves. A new methodology for calculation of the energy difference at trial moves that can be applied to many-body potentials as well as pair ones is a primary feature of our implementation. We study H diffusion in Fe (ferrite phase) and Ni polycrystalline samples to demonstrate the efficiency and scalability of the algorithm and its application. The computational cost of using our framework for half a million atoms is a factor of 250 less than the cost of using existing libraries.
Materials Science,Computational Physics
-
Adaptive tau-leaping methods for microscopic-lattice kinetic Monte Carlo simulations
Tianshi Che,Yang Zhou,Xiaoying Han,Habib N Najm
DOI: https://doi.org/10.1063/5.0218471
2024-08-28
Abstract:Traditional Kinetic Monte Carlo (KMC) approaches, rooted in Gillespie's stochastic simulation algorithm, become computationally demanding in systems with a large range of timescales. The goal of this work is to propose and study new adaptive lattice-KMC time integration strategies for spatially non-uniform systems. To that end, two novel adaptive tau-leaping methods and their corresponding time integration strategies are developed based on the idea of the "n-fold" direct KMC method. These strategies allow for the simultaneous execution of multiple reactions, advancing time by adaptively selected coarse increments. We present numerical experiments comparing the proposed methods with existing approaches in a catalytic surface kinetics application involving ammonia decomposition.
-
Reverse Monte Carlo Simulation: A New Technique for the Determination of Disordered Structures
R. L. McGreevy,L. Pusztai
DOI: https://doi.org/10.1080/08927028808080958
1988-12-01
Molecular Simulation
Abstract:We have developed a new technique, based on the standard Monte Carlo simulation method with Markov chain sampling, in which a set of three dimensional particle configurations are generated that are consistent with the experimentally measured structure factor. A(Q), and radial distribution function, g(r), of a liquid or other disordered system. Consistency is determined by a standard χ2 test using the experimental errors. No input potential is required, we present initial results for liquid argon. Since the technique can work directly from the structure factor it promises to be useful for modelling the structures of glasses or amorphous materials. It also has other advantages in multicomponent systems and as a tool for experimental data analysis.
chemistry, physical,physics, atomic, molecular & chemical
-
Speeding-up Hybrid Functional based Ab Initio Molecular Dynamics using Multiple Time-stepping and Resonance Free Thermostat
Ritama Kar,Sagarmoy Mandal,Vaishali Thakkur,Bernd Meyer,Nisanth N. Nair
2023-09-01
Abstract:Ab initio molecular dynamics (AIMD) based on density functional theory (DFT) has become a workhorse for studying the structure, dynamics, and reactions in condensed matter systems. Currently, AIMD simulations are primarily carried out at the level of generalized gradient approximation (GGA), which is at the 2nd rung of DFT-functionals in terms of accuracy. Hybrid DFT functionals which form the 4th rung in the accuracy ladder, are not commonly used in AIMD simulations as the computational cost involved is $100$ times or higher. To facilitate AIMD simulations with hybrid functionals, we propose here an approach that could speed up the calculations by ~30 times or more for systems with a few hundred of atoms. We demonstrate that, by achieving this significant speed up and making the compute time of hybrid functional based AIMD simulations at par with that of GGA functionals, we are able to study several complex condensed matter systems and model chemical reactions in solution with hybrid functionals that were earlier unthinkable to be performed.
Chemical Physics,Materials Science,Computational Physics
-
A quantitative study on the approximation error and speed-up of the multi-scale MCMC (Monte Carlo Markov chain) method for molecular dynamics
Jie Liu,Qinglin Tang,Jisheng Kou,Dingguo Xu,Tao Zhang,Shuyu Sun
DOI: https://doi.org/10.1016/j.jcp.2022.111491
IF: 4.645
2022-11-15
Journal of Computational Physics
Abstract:The past two decades have borne remarkable progress in the application of the molecular dynamics method in a number of engineering problems. However, the computational efficiency is limited by the massive-atoms system, and the study of rare dynamically-relevant events is challenging at the timescale of molecular dynamics. In this work, a multi-scale molecular simulation algorithm is proposed with a novel toy model that can mimic the state transitions in extensive scenarios. The algorithm consists of two scales, including producing the realistic particle trajectory and probability transition matrix in the molecular dynamics scale and calculating the state distribution and residence time in the Monte Carlo scale. A new state definition is proposed to take the velocity direction into consideration, and different coarsening models are established in the spatial and time scales. The accuracy, efficiency, and robustness of our proposed multi-scale method have been validated, and the general applicability is also demonstrated by explaining two practical applications in the shale gas adsorption and protein folding problems respectively.
physics, mathematical,computer science, interdisciplinary applications
-
Efficient swap algorithms for molecular dynamics simulations of equilibrium supercooled liquids
Ludovic Berthier,Elijah Flenner,Christopher J. Fullerton,Camille Scalliet,Murari Singh
DOI: https://doi.org/10.1088/1742-5468/ab1910
2018-11-30
Abstract:It was recently demonstrated that a simple Monte Carlo (MC) algorithm involving the swap of particle pairs dramatically accelerates the equilibrium sampling of simulated supercooled liquids. We propose two numerical schemes integrating the efficiency of particle swaps into equilibrium molecular dynamics (MD) simulations. We first develop a hybrid MD/MC scheme combining molecular dynamics with the original swap Monte Carlo. We implement this hybrid method in LAMMPS, a software package employed by a large community of users. Secondly, we define a continuous time version of the swap algorithm where both the positions and diameters of the particles evolve via Hamilton's equations of motion. For both algorithms, we discuss in detail various technical issues as well as the optimisation of simulation parameters. We compare the numerical efficiency of all available swap algorithms and discuss their relative merits.
Physics
-
Efficient and accurate simulation of vitrification in multi-component metallic liquids with neural-network potentials
Rui Su,Jieyi Yu,Pengfei Guan,Weihua Wang
DOI: https://doi.org/10.48550/arXiv.2211.03350
2022-11-07
Materials Science
Abstract:Constructing accurate interatomic potential and overcoming the exponential growth of structural equilibration time are challenges to the atomistic investigations of the composition-dependent structure and dynamics during the vitrification process of deeply supercooled multi-component metallic liquids. In this work, we describe a state-of-the-art strategy to address these challenges simultaneously. In the case of the representative Zr-Cu-Al system, in combination with a general algorithm for generating the neural-network potentials (NNP) of multi-component metallic glasses effectively and accurately, we propose a highly efficient atom-swapping hybrid Monte Carlo (SHMC) algorithm for accelerating the thermodynamic equilibration of deeply supercooled liquids. Extensive calculations demonstrate that the newly developed NNP faithfully reproduces the phase stabilities and structural characteristics obtained from the ab initio calculations and experiments. In the combined NNP-SHMC algorithm, the structure equilibration time in the deeply supercooled temperatures is accelerated by at least five orders of magnitudes, and the quenched glassy samples exhibit comparable stability to those prepared in the laboratory. Our results pave the way for the next-generation studies of the vitrification process and, thereby the composition-dependent glass-forming ability and physical properties of multi-component metallic glasses.