Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userOrthorhombic Phase of Crystalline Polyethylene: A Constant Pressure Path Integral Monte Carlo Study
111
0
0.0
(
0
)
Authors
R. Martonak
Ask ChatGPT about the research
No Arabic abstract
In this paper we present a Path Integral Monte Carlo (PIMC) simulation of the orthorhombic phase of crystalline polyethylene, using an explicit atom force field with unconstrained bond lengths and angles. This work represents a quantum extension of our recent classical simulation (J. Chem. Phys. 106, 8918 (1997)). It is aimed both at exploring the applicability of the PIMC method on such polymer crystal systems, as well as on a detailed assessment of the importance of quantum effects on different quantities. We used the $NpT$ ensemble and simulated the system at zero pressure in the temperature range 25 - 300 K, using Trotter numbers between 12 and 144. In order to investigate finite-size effects, we used chains of two different lengths, C_12 and C_24, corresponding to the total number of atoms in the super-cell being 432 and 864, respectively. We show here the results for structural parameters, like the orthorhombic lattice constants a,b,c, and also fluctuations of internal parameters of the chains, such as bond lengths and bond and torsional angles. We have also determined the internal energy and diagonal elastic constants c_11, c_22 and c_33. We discuss the temperature dependence of the measured quantities and compare to that obtained from the classical simulation. For some quantities, we discuss the way they are related to the torsional angle fluctuation. In case of the lattice parameters we compare our results to those obtained from other theoretical approaches as well as to some available experimental data. In order to study isotope effects, we simulated also a deuterated polyethylene crystal at a low temperature. We also suggest possible ways of extending this study and present some general considerations concerning modeling of polymer crystals.
rate research
Read More
In this paper we present a classical Monte Carlo simulation of the orthorhombic phase of crystalline polyethylene, using an explicit atom force field with unconstrained bond lengths and angles and periodic boundary conditions. We used a recently developed algorithm which apart from standard Metropolis local moves employs also global moves consisting of displacements of the center of mass of the whole chains in all three spatial directions as well as rotations of the chains around an axis parallel to the crystallographic c-direction. Our simulations are performed in the NpT ensemble, at zero pressure, and extend over the whole range of temperatures in which the orthorhombic phase is experimentally known to be stable (10 - 450 K). In order to investigate the finite-size effects in this extremely anisotropic crystal, we used different system sizes and different chain lengths, ranging from C_12 to C_96 chains, the total number of atoms in the super-cell being between 432 and 3456. We show here the results for structural parameters, such as the orthorhombic cell parameters a,b,c, and the setting angle of the chains, as well as internal parameters of the chains, such as the bond lengths and angles. Among thermodynamic quantities, we present results for thermal expansion coefficients, elastic constants and specific heat. We discuss the temperature dependence of the measured quantities as well as the related finite-size effects. In case of lattice parameters and thermal expansion coefficients, we compare our results to those obtained from other theoretical approaches as well as to some available experimental data. We also suggest some possible ways of extending this study.
Fractional derivatives are nonlocal differential operators of real order that often appear in models of anomalous diffusion and a variety of nonlocal phenomena. Recently, a version of the Schrodinger Equation containing a fractional Laplacian has been proposed. In this work, we develop a Fractional Path Integral Monte Carlo algorithm that can be used to study the finite temperature behavior of the time-independent Fractional Schrodinger Equation for a variety of potentials. In so doing, we derive an analytic form for the finite temperature fractional free particle density matrix and demonstrate how it can be sampled to acquire new sets of particle positions. We employ this algorithm to simulate both the free particle and $^{4}$He (Aziz) Hamiltonians. We find that the fractional Laplacian strongly encourages particle delocalization, even in the presence of interactions, suggesting that fractional Hamiltonians may manifest atypical forms of condensation. Our work opens the door to studying fractional Hamiltonians with arbitrarily complex potentials that escape analytical solutions.
Quantum Monte Carlo belongs to the most accurate simulation techniques for quantum many-particle systems. However, for fermions, these simulations are hampered by the sign problem that prohibits simulations in the regime of strong degeneracy. The situation changed with the development of configuration path integral Monte Carlo (CPIMC) by Schoof textit{et al.} [T. Schoof textit{et al.}, Contrib. Plasma Phys. textbf{51}, 687 (2011)] that allowed for the first textit{ab initio} simulations for dense quantum plasmas. CPIMC also has a sign problem that occurs when the density is lowered, i.e. in a parameter range that is complementary to traditional QMC formulated in coordinate space. Thus, CPIMC simulations for the warm dense electron gas are limited to small values of the Brueckner parameter -- the ratio of the interparticle distance to the Bohr radius -- $r_s=bar{r}/a_B lesssim 1$. In order to reach the regime of stronger coupling (lower density) with CPIMC, here we investigate additional restrictions on the Monte Carlo procedure. In particular, we introduce two differe
We use the diffusion quantum Monte Carlo (DMC) method to calculate the ground state phase diagram of solid molecular hydrogen and examine the stability of the most important insulating phases relative to metallic crystalline molecular hydrogen. We develop a new method to account for finite-size errors by combining the use of twist-averaged boundary conditions with corrections obtained using the Kwee-Zhang-Krakauer (KZK) functional in density functional theory. To study band-gap closure and find the metallization pressure, we perform accurate quasi-particle many-body calculations using the $GW$ method. In the static approximation, our DMC simulations indicate a transition from the insulating Cmca-12 structure to the metallic Cmca structure at around 375 GPa. The $GW$ band gap of Cmca-12 closes at roughly the same pressure. In the dynamic DMC phase diagram, which includes the effects of zero-point energy, the Cmca-12 structure remains stable up to 430 GPa, well above the pressure at which the $GW$ band gap closes. Our results predict that the semimetallic state observed experimentally at around 360 GPa [Phys. Rev. Lett. {bf 108}, 146402 (2012)] may correspond to the Cmca-12 structure near the pressure at which the band gap closes. The dynamic DMC phase diagram indicates that the hexagonal close packed $P6_3/m$ structure, which has the largest band gap of the insulating structures considered, is stable up to 220 GPa. This is consistent with recent X-ray data taken at pressures up to 183 GPa [Phys. Rev. B {bf 82}, 060101(R) (2010)], which also reported a hexagonal close packed arrangement of hydrogen molecules.
High order actions proposed by Chin have been used for the first time in path integral Monte Carlo simulations. Contrarily to the Takahashi-Imada action, which is accurate to fourth order only for the trace, the Chin action is fully fourth order, with the additional advantage that the leading fourth and sixth order error coefficients are finely tunable. By optimizing two free parameters entering in the new action we show that the time step error dependence achieved is best fitted with a sixth order law. The computational effort per bead is increased but the total number of beads is greatly reduced, and the efficiency improvement with respect to the primitive approximation is approximately a factor of ten. The Chin action is tested in a one-dimensional harmonic oscillator, a H$_2$ drop, and bulk liquid $^4$He. In all cases a sixth-order law is obtained with values of the number of beads that compare well with the pair action approximation in the stringent test of superfluid $^4$He.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
