Browse
Submit Data
Databases
API
Help

Dataset Information

0 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Crystal Structure Prediction for Benzene Using Basin-Hopping Global Optimization.

ABSTRACT: Organic molecules can be stable in distinct crystalline forms, known as polymorphs, which have significant consequences for industrial applications. Here, we predict the polymorphs of crystalline benzene computationally for an accurate anisotropic model parametrized to reproduce electronic structure calculations. We adapt the basin-hopping global optimization procedure to the case of crystalline unit cells, simultaneously optimizing the molecular coordinates and unit cell parameters to locate multiple low-energy structures from a variety of crystal space groups. We rapidly locate all the well-established experimental polymorphs of benzene, each of which corresponds to a single local energy minimum of the model. Our results show that basin-hopping can be both an efficient and effective tool for polymorphic crystal structure prediction, requiring no a priori experimental knowledge of cell parameters or symmetry.

SUBMITTER: Banerjee A

PROVIDER: S-EPMC8279651 | biostudies-literature | 2021 May

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

Crystal Structure Prediction for Benzene Using Basin-Hopping Global Optimization.

Banerjee Atreyee A Jasrasaria Dipti D Niblett Samuel P SP Wales David J DJ

The journal of physical chemistry. A 20210421 17

Organic molecules can be stable in distinct crystalline forms, known as polymorphs, which have significant consequences for industrial applications. Here, we predict the polymorphs of crystalline benzene computationally for an accurate anisotropic model parametrized to reproduce electronic structure calculations. We adapt the basin-hopping global optimization procedure to the case of crystalline unit cells, simultaneously optimizing the molecular coordinates and unit cell parameters to locate mu ...[more]

PMID: 33881850

Similar Datasets

Protein structure prediction using basin-hopping.

Project description:Associative memory Hamiltonian structure prediction potentials are not overly rugged, thereby suggesting their landscapes are like those of actual proteins. In the present contribution we show how basin-hopping global optimization can identify low-lying minima for the corresponding mildly frustrated energy landscapes. For small systems the basin-hopping algorithm succeeds in locating both lower minima and conformations closer to the experimental structure than does molecular dynamics with simulated annealing. For large systems the efficiency of basin-hopping decreases for our initial implementation, where the steps consist of random perturbations to the Cartesian coordinates. We implemented umbrella sampling using basin-hopping to further confirm when the global minima are reached. We have also improved the energy surface by employing bioinformatic techniques for reducing the roughness or variance of the energy surface. Finally, the basin-hopping calculations have guided improvements in the excluded volume of the Hamiltonian, producing better structures. These results suggest a novel and transferable optimization scheme for future energy function development.

| S-EPMC2674628 | biostudies-literature

Systematic Comparison of Genetic Algorithm and Basin Hopping Approaches to the Global Optimization of Si(111) Surface Reconstructions.

Project description:We present a systematic study of two widely used material structure prediction methods, the Genetic Algorithm and Basin Hopping approaches to global optimization, in a search for the 3 × 3, 5 × 5, and 7 × 7 reconstructions of the Si(111) surface. The Si(111) 7 × 7 reconstruction is the largest and most complex surface reconstruction known, and finding it is a very exacting test for global optimization methods. In this paper, we introduce a modification to previous Genetic Algorithm work on structure search for periodic systems, to allow the efficient search for surface reconstructions, and present a rigorous study of the effect of the different parameters of the algorithm. We also perform a detailed comparison with the recently improved Basin Hopping algorithm using Delocalized Internal Coordinates. Both algorithms succeeded in either resolving the 3 × 3, 5 × 5, and 7 × 7 DAS surface reconstructions or getting "sufficiently close", i.e., identifying structures that only differ for the positions of a few atoms as well as thermally accessible structures within kBT/unit area of the global minimum, with T = 300 K. Overall, the Genetic Algorithm is more robust with respect to parameter choice and in success rate, while the Basin Hopping method occasionally exhibits some advantages in speed of convergence. In line with previous studies, the results confirm that robustness, success, and speed of convergence of either approach are strongly influenced by how much the trial moves tend to preserve favorable bonding patterns once these appear.

| S-EPMC9126620 | biostudies-literature

Crystal structure prediction by combining graph network and optimization algorithm.

Project description:Crystal structure prediction is a long-standing challenge in condensed matter and chemical science. Here we report a machine-learning approach for crystal structure prediction, in which a graph network (GN) is employed to establish a correlation model between the crystal structure and formation enthalpies at the given database, and an optimization algorithm (OA) is used to accelerate the search for crystal structure with lowest formation enthalpy. The framework of the utilized approach (a database + a GN model + an optimization algorithm) is flexible. We implemented two benchmark databases, i.e., the open quantum materials database (OQMD) and Matbench (MatB), and three OAs, i.e., random searching (RAS), particle-swarm optimization (PSO) and Bayesian optimization (BO), that can predict crystal structures at a given number of atoms in a periodic cell. The comparative studies show that the GN model trained on MatB combined with BO, i.e., GN(MatB)-BO, exhibit the best performance for predicting crystal structures of 29 typical compounds with a computational cost three orders of magnitude less than that required for conventional approaches screening structures through density functional theory calculation. The flexible framework in combination with a materials database, a graph network, and an optimization algorithm may open new avenues for data-driven crystal structural predictions.

| S-EPMC8938491 | biostudies-literature

Grand and Semigrand Canonical Basin-Hopping.

Project description:We introduce grand and semigrand canonical global optimization approaches using basin-hopping with an acceptance criterion based on the local contribution of each potential energy minimum to the (semi)grand potential. The method is tested using local harmonic vibrational densities of states for atomic clusters as a function of temperature and chemical potential. The predicted global minima switch from dissociated states to clusters for larger values of the chemical potential and lower temperatures, in agreement with the predictions of a model fitted to heat capacity data for selected clusters. Semigrand canonical optimization allows us to identify particularly stable compositions in multicomponent nanoalloys as a function of increasing temperature, whereas the grand canonical potential can produce a useful survey of favorable structures as a byproduct of the global optimization search.

| S-EPMC4750084 | biostudies-other

Search for Global Minimum Structures of P2n+1+ (<i>n</i> = 1-15) Using xTB-Based Basin-Hopping Algorithm.

Project description:A new program for searching global minimum structures of atomic clusters using basin-hopping algorithm based on the xTB method was developed here. The program can be performed with a much higher speed than its replacement directly based on DFT methods. Considering the structural varieties and complexities in finding their global minimum structures, phosphorus cluster cations were studied by the program. The global minimum structures of cationic P2n+1+ (n = 1-15) clusters are determined through the unbiased structure searching method. In the last step, further DFT optimization was performed for the selected isomers. For P2n+1+ (n = 1-4, 7), the found global minimum structures are in consistent with the ones previously reported; while for P2n+1+ (n = 5, 6, 8-12), newly found isomers are more energy-favorable than those previously reported. And those for P2n+1+ (n = 13-15) are reported here for the first time. Among them, the most stable isomers of P2n+1+ (n = 4-6, 9) are characterized by their C3v, Cs, C2v and Cs symmetry, in turn. But those of P2n+1+ (n = 7, 8, 10-12), no symmetry has been identified. The most stable isomers of P29+ and P31+ are characterized by single P-P bonds bridging units inside the clusters. Further analysis shows that the pnicogen bonds play an important role in the stabilization of these clusters. These results show that the new developed program is effective and robust in searching global minimum structures for atom clusters, and it also provides new insights into the role of pnicogen bonds in phosphorus clusters.

| S-EPMC8350033 | biostudies-literature

Protein structure prediction by global optimization of a potential energy function.

Project description:An approach based exclusively on finding the global minimum of an appropriate potential energy function has been used to predict the unknown structures of five globular proteins with sizes ranging from 89 to 140 amino acid residues. Comparison of the computed lowest-energy structures of two of them (HDEA and MarA) with the crystal structures, released by the Protein Data Bank after the predictions were made, shows that large fragments (61 residues) of both proteins were predicted with rms deviations of 4.2 and 6.0 A for the Calpha atoms, for HDEA and MarA, respectively. This represents 80% and 53% of the observed structures of HDEA and MarA, respectively. Similar rms deviations were obtained for approximately 60-residue fragments of the other three proteins. These results constitute an important step toward the prediction of protein structure based solely on global optimization of a potential energy function for a given amino acid sequence.

| S-EPMC21885 | biostudies-literature

An improved hybrid global optimization method for protein tertiary structure prediction.

Project description:First principles approaches to the protein structure prediction problem must search through an enormous conformational space to identify low-energy, near-native structures. In this paper, we describe the formulation of the tertiary structure prediction problem as a nonlinear constrained minimization problem, where the goal is to minimize the energy of a protein conformation subject to constraints on torsion angles and interatomic distances. The core of the proposed algorithm is a hybrid global optimization method that combines the benefits of the αBB deterministic global optimization approach with conformational space annealing. These global optimization techniques employ a local minimization strategy that combines torsion angle dynamics and rotamer optimization to identify and improve the selection of initial conformations and then applies a sequential quadratic programming approach to further minimize the energy of the protein conformations subject to constraints. The proposed algorithm demonstrates the ability to identify both lower energy protein structures, as well as larger ensembles of low-energy conformations.

| S-EPMC2847311 | biostudies-literature

mRNA secondary structure optimization using a correlated stem-loop prediction.

Project description:Secondary structure of messenger RNA plays an important role in the bio-synthesis of proteins. Its negative impact on translation can reduce the yield of protein by slowing or blocking the initiation and movement of ribosomes along the mRNA, becoming a major factor in the regulation of gene expression. Several algorithms can predict the formation of secondary structures by calculating the minimum free energy of RNA sequences, or perform the inverse process of obtaining an RNA sequence for a given structure. However, there is still no approach to redesign an mRNA to achieve minimal secondary structure without affecting the amino acid sequence. Here we present the first strategy to optimize mRNA secondary structures, to increase (or decrease) the minimum free energy of a nucleotide sequence, without changing its resulting polypeptide, in a time-efficient manner, through a simplistic approximation to hairpin formation. Our data show that this approach can efficiently increase the minimum free energy by >40%, strongly reducing the strength of secondary structures. Applications of this technique range from multi-objective optimization of genes by controlling minimum free energy together with CAI and other gene expression variables, to optimization of secondary structures at the genomic level.

| S-EPMC3616703 | biostudies-literature

Crystal structure of 1,3,5-trimethyl-2,4-di-nitro-benzene.

Project description:In the title compound, C9H10N2O4, the planes of the nitro groups subtend dihedral angles of 55.04 (15) and 63.23 (15)° with that of the aromatic ring. These tilts are in opposite senses and the mol-ecule possesses approximate mirror symmetry about a plane normal to the mol-ecule. In the crystal, mol-ecules form stacks in the [100] direction with adjacent mol-ecules related by translation, although the centroid-centroid separation of 4.136 (5) Å is probably too long to regard as a significant aromatic π-π stacking inter-action. An extremely weak C-H⋯O inter-action is also present.

| S-EPMC4555366 | biostudies-literature

Crystal structure of the bora-benzene-2,6-lutidine adduct.

Project description:In the title compound, C12H14BN, the complete mol-ecule is generated by a crystallographic twofold axis, with two C atoms, the B atom and the N atom lying on the rotation axis. The dihedral angle between the bora-benzene and pyridine rings is 81.20 (6)°. As well as dative electron donation from the N atom to the B atom [B-N = 1.5659 (18) Å], the methyl substituents on the lutidine ring shield the B atom, which further stabilizes the mol-ecule. In the crystal, weak aromatic π-π stacking between the pyridine rings [centroid-centroid separation = 3.6268 (9) Å] is observed, which generates [001] columns of mol-ecules.

| S-EPMC4719914 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data