Implications of the multi-minima character of molecular crystal phases onto the free energy

Ehsan Rahmatizad Khajehpasha; Giuseppe Fisicaro; Jonas A. Finkler; Marco Krummenacher; Martin Sommer-J\"orgensen; Moritz Gubler; Stefan Goedecker

REVIEW 4 major objections 5 minor 53 references

Implications of the multi-minima character of molecular crystal phases onto the free energy

T0 review · 4 major / 5 minor · reviewed 2026-08-09 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Counting every crystal minimum shifts a phase transition by 200 K

desk verdict A credible, important demonstration that configurational entropy can shift polymorph transition temperatures by ~200 K, but the central number needs stronger validation and a check of the harmonic-basin assumption. read the letter →

arxiv 2501.18372 v1 pith:QWSNB7EM submitted 2025-01-30 cond-mat.mtrl-sci physics.comp-ph

Marco Krummenacher , Martin Sommer-Jörgensen , Moritz Gubler , Jonas A. Finkler , Ehsan Rahmatizad Khajehpasha , Giuseppe Fisicaro , Stefan Goedecker This is my paper · ORCID

classification cond-mat.mtrl-sciphysics.comp-ph

keywords molecularcrystalspolymorphismconfigurationalentropycrystalstructurepredictionfreeenergyquantumsuperpositionmethodpotentiallandscapesolid-solidphasetransition

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper sets out to show that a polymorph of a molecular crystal is not a single structure but an ensemble of many local minima on the potential energy surface, and that the configurational entropy of this ensemble changes which phase is thermodynamically stable at finite temperature. Using the molecule N-(4-methylbenzylidene)-4-methylalanine, the authors find 304 distinct minima belonging to its three known forms, with the dense, low-energy Form III having only 84 minima while the looser Forms I and II have many more. Combining these minima through the quantum superposition method, the Boltzmann crossing between Form II and Form III moves from roughly 470 K to about 270 K when all sampled minima are included instead of just the lowest-energy structure of each phase. The authors conclude that neglecting configurational entropy introduces errors of about 200 K even in a relatively simple molecular crystal, and that finite-temperature crystal structure prediction must therefore include it.

What carries the argument

The load-bearing object is the quantum superposition method, which builds the total partition function as a sum over all sampled local minima: $Z = \sum_j \exp(-\beta F_j(\beta))$, where $F_j(\beta)$ is the harmonic free energy of minimum $j$. This turns each phase into a weighted set of minima rather than a single structure. The machinery also includes: Minima Hopping at DFT level to find the 304 minima; a NequIP machine-learned potential fitted to DFT data, which makes molecular dynamics and vibrational frequencies affordable; two collective variables (angles between ring-link vectors and ring normal vectors) that classify the minima into the three experimental forms; and a configurational density of states, shown in Fig. 4, which quantifies how many minima each form has in each energy window. The method is valid when the system can visit all low-energy basins on experimental timescales but crossings between basins are rare events; the paper verifies this with MD timings.

What would settle it

Compute the density of states of each form more exhaustively, for example by enumerating all minima within 500 meV of the global minimum using longer enhanced-sampling dynamics, and recompute the Boltzmann crossing temperature; if the crossing stays at about 270 K, the 200 K shift is robust, whereas if it moves back toward 470 K the truncation of unsampled minima was the deciding factor.

Watch

Extended reading notes

Core claim

The central claim is that the multi-minima character of the phases is large enough to reverse free-energy rankings. For N-(4-methylbenzylidene)-4-methylalanine, Form III is the global minimum in potential energy with the highest density, but it is structurally rigid: rotating methyl groups or sliding layers is costly, so it has few nearby low-energy minima. Form II and Form I are structurally tolerant, generating 157 and 63 minima respectively at modest energy costs, and this difference in configurational density of states outweighs their higher potential energies at elevated temperature. Using the quantum superposition partition function, the probability of finding Form II overtakes Form III at about 270 K when all 304 minima are counted, whereas using only the lowest structures per phase puts the crossing at about 470 K. The paper argues that the entropy of the multi-minima ensemble is essential for accurate finite-temperature free energies and for explaining why 'lowest energy' ranking fails to predict the experimentally relevant phase.

Load-bearing premise

The 304 minima found by the search, and the assumption that no unsampled higher-energy minima contribute, faithfully represent the phase space of the three forms at the temperatures of interest.

Editorial extensions

If this is right

Polymorph rankings based on the lowest-energy structure alone can be wrong at finite temperature; the correct ranking must sum over the minima of each phase.
Transition temperatures between polymorphs can be off by about 200 K when configurational entropy is neglected, so crystal structure prediction benchmarks that ignore it inherit this error.
Structurally tolerant, less dense phases are stabilized by configurational entropy at high temperature, acting together with their larger vibrational entropy.
The quantum superposition method, combined with a machine-learned potential, is a practical route to finite-temperature free energies for molecular crystals.
For more complex molecular crystals, where the multi-minima character is presumably more pronounced, the effect on phase stability should be even larger.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

This suggests 'over-prediction' in crystal structure prediction is partly a thermodynamic statement: many predicted minima are genuinely occupied at finite temperature, and what matters is the free energy of the ensemble, not the depth of one basin.
If correct, the rank-ordering of polymorphs by computed lattice energies should be accompanied by a measure of landscape 'tolerance' (e.g., density of low-energy minima) as a screening descriptor.
The same reasoning may apply to the phenomenon of disappearing polymorphism: a phase that is kinetically accessible but entropically disfavored could vanish from recrystallization even if it is the global energy minimum.
A testable extension is to run the same quantum-superposition analysis on other molecules with known polymorph pairs, checking whether the 200 K-scale shift generalizes.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Colleague,

The headline: this paper gives a concrete number for a real effect—configurational entropy from multiple minima shifts a polymorph transition temperature by about 200 K in a molecular crystal. The qualitative idea has been around (Dybeck, Butler, Francia), but the quantitative estimate for N-(4-methylbenzylidene)-4-methylalanine is new. If it holds up, it strengthens the case that finite-temperature CSP must include multi-minima contributions.

What's good: the workflow is logical. Minima hopping on DFT, re-relaxation with PBE-D4, a NequIP potential for frequencies and MD, classification by collective variables, and a quantum-superposition partition function. The density of states per form (Fig. 4) makes the entropic argument visually clear. The convergence check with 100 vs 300 minima (290 vs 270 K) is a reasonable attempt. Exact saddle points are a nice addition.

The soft spots are real, though. The biggest one is the validity of the harmonic-superposition sum at the claimed crossing. The paper's own condition is that transitions between basins are rare events. But their barriers between Form I/II minima are 10–100 meV per molecule, and kBT at 270 K is about 23 meV. Their MD shows hops every few picoseconds across many distinct minima (Table I). That does not look like a rare-event regime. The harmonic partition function may be overcounting entropy that an anharmonic treatment would handle correctly, and the 200 K could be partly an artifact of the basin decomposition. The stress-test note is on target.

Second, the density of states is truncated at high energies—acknowledged in the Fig. 4 caption. The 100-vs-300 convergence is not a bound on unsampled minima. The direction is probably to increase the effect, but it's unquantified. Third, no error bars on the 200 K. Fourth, the NequIP potential is not validated against DFT for the specific observables (frequencies, MD barriers), so the whole free-energy calculation rests on an unbenchmarked model. Fifth, no code or data.

None of this kills the qualitative conclusion. The tolerant forms have more minima and lower frequencies, so they should be stabilized at high T. But the exact 200 K is not yet a reliable number. The paper deserves a serious referee; it's a methodologically important demonstration. The referee should ask for ML-potential validation, a check of the rare-event assumption (maybe via anharmonic free energies or longer MD), uncertainty quantification, and, ideally, release of the minima set and code.

Referee Report

4 major / 5 minor

Summary. The paper argues that the multi-minima character of molecular crystal potential energy surfaces contributes a significant configurational entropy that must be included in finite-temperature free-energy rankings. For N-(4'-methylbenzylidene)-4-methylalanine, the authors use minima hopping at the DFT level to sample 304 local minima, classify them into three experimentally known forms using collective variables, and compute harmonic free energies. They then apply the quantum superposition method (Eq. 5) to sum Boltzmann weights over all minima and obtain phase probabilities as a function of temperature. The central claim is that the Form II–III transition temperature moves from about 470 K when only the lowest minimum of each phase is used to about 270 K when all 304 minima are included, implying a roughly 200 K error if configurational entropy is neglected. The paper also reports exact saddle-point barriers along transformation pathways and MD-based estimates of inter-basin hopping times.

Significance. If correct, the result is important for crystal structure prediction: it would demonstrate that configurational entropy from multiple local minima can change polymorphic free-energy rankings by hundreds of kelvin, an effect currently ignored in standard CSP protocols. The study is valuable for its DFT-level minima-hopping sampling, exact saddle-point calculations, and the explicit connection to a crystal with a known disappearing-polymorphism problem. The methodology is standard but cleanly applied, and the observation that different forms have markedly different configurational densities of states is an interesting qualitative finding. However, the quantitative claim lacks uncertainty estimates and its central approximation is applied in a regime where the paper's own data indicate it may not be valid.

major comments (4)

[Eq. (5) and text after it] The paper states that the superposition method is valid only if transitions between neighboring catchment basins are rare events, but the data in Table I and Fig. 6 contradict this condition at the temperatures of interest. Table I shows inter-basin hops on a timescale of about 10 ps, while vibrational periods are reported to extend to about 1 ps; the 10–100 meV barriers in Fig. 6 are comparable to kBT ≈ 23 meV at the claimed 270 K crossing. Under these conditions, the harmonic superposition sum in Eq. (5) can double-count anharmonic contributions that a proper anharmonic free-energy calculation would already include, and the 200 K shift in Fig. 5 may be an artifact of the basin decomposition. I recommend that the authors demonstrate the validity of the rare-event harmonic description at the crossing temperature, for example by comparing with an anharmonic free-energy calculation for the lowest-energy minima or by showing that barriers are large compared to kBT over the relevant temperature range.
[Fig. 4 and Eq. (5)] The configurational density of states shown in Fig. 4 is truncated to zero above the sampled energy range, as explicitly acknowledged in the caption. The convergence test between 100 and 300 structures (290 K vs 270 K) only varies the number of sampled minima and does not bound the contribution of unsampled higher-energy minima. Since the partition function in Eq. (5) is sensitive to the high-energy tail of the DOS, the reported 'about 200 K' error has no estimated uncertainty and may be an underestimate. The authors should either extrapolate the DOS, estimate the missing contribution, or provide a quantitative bound on the error introduced by the truncation.
[Machine-learned potential section] The NequIP machine-learned potential is used to compute vibrational frequencies and to run longer MD, but the manuscript presents no validation of this potential against DFT (e.g., energy/force errors on a held-out test set). Because the harmonic free energies and the number of distinct minima visited depend directly on the accuracy of this potential, the quantitative claim currently lacks support. Please add a validation of the NequIP potential and an estimate of how the remaining DFT/ML discrepancy propagates into the computed transition temperatures.
[Fig. 5] The headline result of a 200 K error is presented as a single number without any uncertainty or sensitivity analysis. The 470 K and 270 K crossing temperatures are computed with specific choices of functional (PBE+D4), harmonic approximation, and a finite sample of 304 minima, but no error bar is given. At minimum, the authors should provide a crude estimate of the uncertainty based on, for example, variations of the DOS truncation or the ML potential accuracy, so that the central number can be meaningfully interpreted.

minor comments (5)

[Eq. (1) and surrounding text] There is a notation inconsistency: the text describes a vector '®𝑤' spanned between the centers of mass of the two carbon rings, but Eq. (1) then uses '®𝑣' in the angle definition; please unify the notation.
[Table I caption] The caption contains a typo: 'The ration' should be 'The ratio'. Please correct.
[Title of molecule] The molecule name appears inconsistently as 'N-(4-Methylbenzylidene)-4-methylalanine' in the abstract and 'N-(4’-Methylbenzylidene)-4-methylalanine' in the body; please use one name consistently.
[Methods details] The manuscript does not specify how the harmonic free energies in Eq. (5) are computed, for example whether the vibrational frequencies are obtained at the Γ point only or with a k-point mesh, and whether corrections for imaginary modes are applied. This information is needed to reproduce the results.
[Fig. 1 caption] The caption has a typo: 'Methylkbenzylidene' should be 'Methylbenzylidene'.

Circularity Check

0 steps flagged · score 2.0 of 10

No circularity in the central derivation; the 200 K shift is a computed consequence of the sampled multi-minima density of states, with only minor non-load-bearing self-citations.

full rationale

The central derivation is self-contained: the paper samples 304 minima with Minima Hopping on a DFT-level potential energy surface, relaxes them with PBE-D4, trains a NequIP machine-learned potential on DFT data for dynamics and frequencies, computes harmonic free energies per catchment basin, forms the superposition partition function in Eq. 5, and obtains Boltzmann probabilities from Eq. 6. The reported crossing temperatures (470 K vs 270 K) are outputs of this sum, not parameters fitted to reproduce any target value; no term in Eqs. 5-6 is adjusted to force the crossing. The acknowledged truncation of the density of states at high energies (Fig. 4 caption) and the convergence check using 100 vs 300 minima are limitations of sampling, not circular inputs. Self-citations, including Ref. [29] for an analogous 200 K effect and Refs. [34-37,51,52] for Minima Hopping, COMPASS, and fingerprints, are methodological or motivational; they provide tools and context, but the quantitative claim rests on computed energies and basin counts. The harmonic/rare-event assumption is a physical correctness risk rather than a circular reduction. At most there are minor non-load-bearing self-citations, so the circularity score is low.

Assumptions & free parameters 1 free parameters · 5 assumptions · 0 invented entities

The central free energy calculation rests on the validity of the quantum superposition sum, on the completeness of the sampled minima, and on the accuracy of the DFT/ML potential energy surface. The paper checks convergence of the transition temperature with respect to the number of minima (100 vs 300) but does not quantify the truncation error in the density of states or the uncertainty from the ML potential, which is the main reason the numerical claim (200 K shift) carries moderate risk.

free parameters (1)

NequIP neural network potential weights = not provided
The machine learned potential is fitted to DFT data from Minima Hopping and is used to compute vibrational free energies and run MD; no accuracy metrics or training set details are given, so its numerical values are undisclosed inputs to the central free energy calculation.

assumptions (5)

domain assumption Harmonic approximation for the vibrational free energy of each local minimum.
Quantum superposition method computes F_j as harmonic free energy (Eq. 5 and surrounding text); anharmonic effects are neglected and not quantified.
domain assumption The system visits all low-energy catchment basins on experimental timescales, and inter-basin transitions are rare events.
Stated explicitly as validity conditions for the quantum superposition method in the paragraph after Eq. 6; MD at 300-400 K (Table I) is used as support but does not prove equilibration across all relevant basins.
ad hoc to paper The 304 sampled minima from Minima Hopping adequately represent the configurational density of states of all three phases.
Convergence is checked by comparing 100 vs 300 structures (290 K vs 270 K), but the DOS is truncated at high energies (Fig. 4 caption), so the true phase-space weight of unsampled minima is unknown.
domain assumption Dispersion-corrected DFT (PBE-D4, with LDA for structure search) provides the reference energies and forces that both the minima enumeration and the ML potential inherit.
All minima and the NequIP training data come from DFT; no comparison to experimental lattice energies or coupled-cluster benchmarks is given.
domain assumption The two collective variable fingerprints (alpha and beta angles) correctly classify structures into the three experimental forms.
Classification uses distances to one reference structure per phase (Eqs. 1-4); no sensitivity analysis to the choice of reference or threshold is presented.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Implications of the multi-minima character of molecular crystal phases onto the free energy." pith.science (2026). https://pith.science/paper/QWSNB7EM

@misc{pith2026250118372,
  author       = {Pith},
  title        = {Pith review of: Implications of the multi-minima character of molecular crystal phases onto the free energy},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/QWSNB7EM}},
  note         = {Machine review of arXiv:2501.18372}
}

read the original abstract

In recent years, significant advancements in computational methods have dramatically enhanced the precision in determining the energetic ranking of different phases of molecular crystals. The developments mainly focused on providing accurate dispersion corrected exchange correlation functionals and methods for describing the vibrational entropy contributions to the free energy at finite temperatures. Several molecular crystals phases were recently found to have of multi-minima character. For our investigations we highlight the multi-minima character in the example of the molecular crystal consisting of N-(4-Methylbenzylidene)-4-methylalanine. We explore its potential energy landscape on the full DFT level or with a machine learned potential that was fitted to DFT data. We calculate not only many local minima but also exact barriers along transformation pathways to demonstrate the multi-minima character of our system. Furthermore, we present a framework, based on the quantum superposition method, that includes both configurational and vibrational entropy. As an example, we show for our system that the transition temperature between two of its phases is afflicted by an error of about 200 K if the multi-minima character is not taken into account. This indicates that it is absolutely essential to consider configurational entropy to obtain reliable finite temperature free energy rankings for complex molecular crystals.

Figures

Figures reproduced from arXiv: 2501.18372 by the authors.

**Figure 2.** Vectors spanned to distinguish between the different poly [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗

**Figure 3.** Energy-density diagram of the 304 structures found with [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗

Figures from the paper (2 more)

**Figure 5.** Figure 5: Boltzmann probabilities of finding a certain form of N- [PITH_FULL_IMAGE:figures/full_fig_p004_5.png]

**Figure 6.** Figure 6: Representative reaction pathway between structures of the [PITH_FULL_IMAGE:figures/full_fig_p004_6.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

53 extracted references · 50 canonical work pages

[1]

A. J. Cruz-Cabeza, S. M. Reutzel-Edens, and J. Bernstein, Facts and fictions about polymorphism, Chem. Soc. Rev.44, 8619 (2015)

work page 2015
[2]

Q. Zhu, J. Johal, D. E. Widdowson, Z. Pang, B. Li, C. M. Kane, V. Kurlin, G. M. Day, M. A. Little, and A. I. Cooper, Analogy poweredbypredictionandstructuralinvariants: computationally led discovery of a mesoporous hydrogen-bonded organic cage crystal, J. Am. Chem. Soc.144, 9893 (2022)

work page 2022
[3]

P. Cui, D. P. McMahon, P. R. Spackman, B. M. Alston, M. A. Little, G. M. Day, and A. I. Cooper, Mining predicted crystal structure landscapes with high throughput crystallisation: old molecules, new insights, Chemical Science10, 9988 (2019)

work page 2019
[4]

Neumann, J

M. Neumann, J. Van De Streek, F. Fabbiani, P. Hidber, and O. Grassmann, Combined crystal structure prediction and high- pressure crystallization in rational pharmaceutical polymorph screening, Nature communications6, 7793 (2015)

work page 2015
[5]

A. G. Shtukenberg, Q. Zhu, D. J. Carter, L. Vogt, J. Hoja, E. Schneider, H. Song, B. Pokroy, I. Polishchuk, A. Tkatchenko, et al.,Powderdiffractionandcrystalstructurepredictionidentify four new coumarin polymorphs, Chemical Science8, 4926 (2017)

work page 2017
[6]

D. E. Braun, M. Ardid-Candel, E. D’Oria, P. G. Karamertzanis, J.-B.Arlin, A.J.Florence, A.G.Jones,andS.L.Price,Racemic naproxen: a multidisciplinary structural and thermodynamic comparison with the enantiopure form, Crystal growth & design 11, 5659 (2011)

work page 2011
[7]

G. M. Day, Current approaches to predicting molecular organic crystal structures, Crystallography Reviews17, 3 (2011)

work page 2011
[8]

Nyman and S

J. Nyman and S. M. Reutzel-Edens, Crystal structure prediction is changing from basic science to applied technology, Faraday Discussions 211, 459 (2018)

work page 2018

Show all 53 references

[9]

T.S.Thakur,R.Dubey,andG.R.Desiraju,Crystalstructureand prediction, Annu. Rev. Phys. Chem.66, 21 (2015)

2015
[10]

Mortazavi, J

M. Mortazavi, J. Hoja, L. Aerts, L. Quéré, J. van de Streek, M. A. Neumann, and A. Tkatchenko, Computational polymorph screening reveals late-appearing and poorly-soluble form of rotigotine, Communications Chemistry2, 70 (2019)

2019
[11]

Habgood, Form ii caffeine: a case study for confirming and predicting disorder in organic crystals, Crystal growth & design 11, 3600 (2011)

M. Habgood, Form ii caffeine: a case study for confirming and predicting disorder in organic crystals, Crystal growth & design 11, 3600 (2011)

2011
[12]

E. L. Stevenson, R. W. Lancaster, A. B. Buanz, L. S. Price, D. A. Tocher,andS.L.Price,Thesolidstateformsofthesexhormone 17-𝛽-estradiol, CrystEngComm21, 2154 (2019)

2019
[13]

D. H. Case, J. E. Campbell, P. J. Bygrave, and G. M. Day, Convergence properties of crystal structure prediction by quasi- random sampling, J. Chem. Theory Comput.12, 910 (2016)

2016
[14]

M. A. Neumann, F. J. Leusen, and J. Kendrick, A major advance incrystalstructureprediction,AngewandteChemieInternational Edition 47, 2427 (2008)

2008
[15]

Kronik and A

L. Kronik and A. Tkatchenko, Understanding molecular crystals with dispersion-inclusive density functional theory: pairwise corrections and beyond, Acc. Chem. Res.47, 3208 (2014)

2014
[16]

L. M. Hunnisett, J. Nyman, N. Francia, N. S. Abraham, C. S. Adjiman, S. Aitipamula, T. Alkhidir, M. Almehairbi, A. Anelli, D. M. Anstine,et al., The seventh blind test of crystal structure prediction: structure generation methods, Structural Science80, https://doi.org/10.1107/...

2024 doi
[17]

E. C. Dybeck, D. P. McMahon, G. M. Day, and M. R. Shirts, Exploring the multi-minima behavior of small molecule crystal polymorphs at finite temperature, Crystal Growth & Design19, 5568 (2019)

2019
[18]

Schön, H

J. Schön, H. Putz, and M. Jansen, Studying the energy hypersur- face of continuous systems-the threshold algorithm, J. Condens. Matter Phys8, 143 (1996)

1996
[19]

P.W.ButlerandG.M.Day,Reducingoverpredictionofmolecular crystal structures via threshold clustering, Proceedings of the National Academy of Sciences120, e2300516120 (2023)

2023
[20]

Yang and G

S. Yang and G. M. Day, Global analysis of the energy land- scapes of molecular crystal structures by applying the threshold algorithm, Commun. Chem.5, 86 (2022)

2022
[21]

Montis, M

R. Montis, M. B. Hursthouse, J. Kendrick, J. Howe, and R. J. Whitby, Combining structural rugosity and crystal packing com- parison: Aroutetomorepolymorphs?,CrystalGrowth&Design 22, 559 (2021)

2021
[22]

T.Beyer,G.M.Day,andS.L.Price,Theprediction,morphology, and mechanical properties of the polymorphs of paracetamol, J. Am. Chem. Soc.123, 5086 (2001)

2001
[23]

Montis, R

R. Montis, R. J. Davey, S. E. Wright, G. R. Woollam, and A. J. Cruz-Cabeza, Transforming computed energy landscapes into experimentalrealities: theroleofstructuralrugosity,Angewandte Chemie 132, 20537 (2020)

2020
[24]

N. F. Francia, L. S. Price, J. Nyman, S. L. Price, and M. Sal- valaglio, Systematic finite-temperature reduction of crystal en- ergy landscapes, Crystal Growth & Design20, 6847 (2020)

2020
[25]

Schneider, L

E. Schneider, L. Vogt, and M. E. Tuckerman, Exploring poly- morphism of benzene and naphthalene with free energy based enhanced molecular dynamics, Acta Crystallographica Section 6 B:StructuralScience,CrystalEngineeringandMaterials 72,542 (2016)

2016
[26]

I. J. Sugden, N. F. Francia, T. Jensen, C. S. Adjiman, and M.Salvalaglio,Rationalisingthedifferenceincrystallisabilityof two sulflowers using efficient in silico methods, CrystEngComm 24, 6830 (2022)

2022
[27]

H. Song, L. Vogt-Maranto, R. Wiscons, A. J. Matzger, and M. E. Tuckerman, Generating cocrystal polymorphs with information entropydrivenbymoleculardynamics-basedenhancedsampling, J. Phys. Chem. Lett.11, 9751 (2020)

2020
[28]

N. F. Francia, L. S. Price, and M. Salvalaglio, Reducing crystal structure overprediction of ibuprofen with large scale molecular dynamics simulations, CrystEngComm23, 5575 (2021)

2021
[29]

Krummenacher, O

M. Krummenacher, O. Tayfuroglu, J. Finkler, H. Huber, and S. Goedecker, Entropic stabilization of a structurally tolerant phase: The ionic phase of lithium alanate, Research Square Preprints https://doi.org/10.21203/rs.3.rs-4318358/v1 (2024)

2024 doi
[30]

Wales, Energy Landscapes: Applications to Clusters, Biomolecules and Glasses, 1st ed

D. Wales, Energy Landscapes: Applications to Clusters, Biomolecules and Glasses, 1st ed. (Cambridge University Press, 2004)

2004
[31]

Bar and J

I. Bar and J. Bernstein, Molecular conformation and electronic structure. vi. the structure of p-methyl-n-(p-methylbenzylidene) aniline (form i), Acta Crystallographica Section B: Structural Crystallography and Crystal Chemistry38, 121 (1982)

1982
[32]

Bar and J

I. Bar and J. Bernstein, Molecular conformation and elec- tronic structure. v. the crystal and molecular structure of n- (p-methylbenzylidine)-p-methylaniline (form ii), Acta Crystal- lographica Section B: Structural Crystallography and Crystal Chemistry 33, 1738 (1977)

1977
[33]

Bernstein, I

J. Bernstein, I. Bar, and A. Christensen, Molecular confor- mation and electronic structure. iv. p-(n-methylbenzylidene)- p-methylaniline (form iii), Acta Crystallographica Section B: Structural Crystallography and Crystal Chemistry 32, 1609 (1976)

1976
[34]

Goedecker, Minima hopping: An efficient search method for the global minimum of the potential energy surface of complex molecular systems, J

S. Goedecker, Minima hopping: An efficient search method for the global minimum of the potential energy surface of complex molecular systems, J. Chem. Phys.120, 9911 (2004)

2004
[35]

Amsler and S

M. Amsler and S. Goedecker, Crystal structure prediction us- ing the minima hopping method, The Journal of Chemical Physics 133, 224104 (2010), https://pubs.aip.org/aip/jcp/article- pdf/doi/10.1063/1.3512900/15433393/224104_1_online.pdf

2010 doi
[36]

Krummenacher, M

M. Krummenacher, M. Gubler, J. A. Finkler, H. Huber, M. Sommer-Jörgensen, and S. Goedecker, Performing highly efficient minima hopping structure predictions using the atomic simulation environment (ase), SoftwareX25, 101632 (2024)

2024
[37]

Gubler, M

M. Gubler, M. Krummenacher, H. Huber, and S. Goedecker, Efficient variable cell shape geometry optimization, J. Comput. Phys. X17, 100131 (2023)

2023
[38]

S. A. Ghasemi, M. Amsler, R. G. Hennig, S. Roy, S. Goedecker, T. J. Lenosky, C. J. Umrigar, L. Genovese, T. Morishita, and K.Nishio,Energylandscapeofsiliconsystemsanditsdescription byforcefields,tightbindingschemes,densityfunctionalmethods, and quantum monte carlo methods, Phy...

2010
[39]

A.Kozhevnikov,M.Taillefumier,andS.Pintarelli,Sirius, https: //github.com/electronic-structure/SIRIUS (2021)

2021
[40]

Goedecker, M

S. Goedecker, M. Teter, and J. Hutter, Separable dual-space gaussian pseudopotentials, Phys. Rev. B54, 1703 (1996)

1996
[41]

A. J. Price, K. R. Bryenton, and E. R. Johnson, Requirements for an accurate dispersion-corrected density functional, J. Chem. Phys. 154, https://doi.org/10.1063/5.0050993 (2021)

2021 doi
[42]

G. J. Beran, Modeling polymorphic molecular crystals with electronic structure theory, Chem. Rev.116, 5567 (2016)

2016
[43]

Bannwarth, and S

E.Caldeweyher,S.Ehlert,A.Hansen,H.Neugebauer,S.Spicher, C. Bannwarth, and S. Grimme, A generally applicable atomic- charge dependent london dispersion correction, J. Chem. Phys. 150, https://doi.org/10.1063/1.5090222 (2019)

2019 doi
[44]

Caldeweyher, J.-M

E. Caldeweyher, J.-M. Mewes, S. Ehlert, and S. Grimme, Ex- tension and evaluation of the d4 london-dispersion model for periodic systems, PCCP22, 8499 (2020)

2020
[45]

J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Phys. Rev. Lett.77, 3865 (1996)

1996
[46]

Batzner, A

S. Batzner, A. Musaelian, L. Sun, M. Geiger, J. P. Mailoa, M. Kornbluth, N. Molinari, T. E. Smidt, and B. Kozinsky, E (3)- equivariant graph neural networks for data-efficient and accurate interatomic potentials, Nature communications13, 2453 (2022)

2022
[47]

A. H. Larsen, J. J. Mortensen, J. Blomqvist, I. E. Castelli, R. Christensen, M. Dułak, J. Friis, M. N. Groves, B. Hammer, C. Hargus,et al., The atomic simulation environment—a python library for working with atoms, J. Condens. Matter Phys.29, 273002 (2017)

2017
[48]

S. L. Price, Control and prediction of the organic solid state: a challenge to theory and experiment, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences 474, 20180351 (2018)

2018
[49]

S. L. Price and S. M. Reutzel-Edens, The potential of computed crystal energy landscapes to aid solid-form development, Drug Discovery Today21, 912 (2016)

2016
[50]

Schaefer, S

B. Schaefer, S. Mohr, M. Amsler, and S. Goedecker, Min- ima hopping guided path search: An efficient method for find- ing complex chemical reaction pathways, J. Chem. Phys.140, https://doi.org/10.1063/1.4878944 (2014)

2014 doi
[51]

Sommer-Jörgensen and S

M. Sommer-Jörgensen and S. Goedecker, Compass: Double- endedsaddlepointsearchasaconstrainedoptimizationproblem, J. Chem. Phys.160, https://doi.org/10.1063/5.0186903 (2024)

2024 doi
[52]

L. Zhu, M. Amsler, T. Fuhrer, B. Schaefer, S. Faraji, S. Rostami, S.A.Ghasemi,A.Sadeghi,M.Grauzinyte,C.Wolverton, et al.,A fingerprint based metric for measuring similarities of crystalline structures, J. Chem. Phys.144, 034203 (2016)

2016
[53]

A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in ’t Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton, Lammps - a flexible simulation tool for particle-based mater...

2022

Pith tools

Reviewed August 9, 2026 · model on record in the stance chip above.

[1] [1]

A. J. Cruz-Cabeza, S. M. Reutzel-Edens, and J. Bernstein, Facts and fictions about polymorphism, Chem. Soc. Rev.44, 8619 (2015)

work page 2015

[2] [2]

Q. Zhu, J. Johal, D. E. Widdowson, Z. Pang, B. Li, C. M. Kane, V. Kurlin, G. M. Day, M. A. Little, and A. I. Cooper, Analogy poweredbypredictionandstructuralinvariants: computationally led discovery of a mesoporous hydrogen-bonded organic cage crystal, J. Am. Chem. Soc.144, 9893 (2022)

work page 2022

[3] [3]

P. Cui, D. P. McMahon, P. R. Spackman, B. M. Alston, M. A. Little, G. M. Day, and A. I. Cooper, Mining predicted crystal structure landscapes with high throughput crystallisation: old molecules, new insights, Chemical Science10, 9988 (2019)

work page 2019

[4] [4]

Neumann, J

M. Neumann, J. Van De Streek, F. Fabbiani, P. Hidber, and O. Grassmann, Combined crystal structure prediction and high- pressure crystallization in rational pharmaceutical polymorph screening, Nature communications6, 7793 (2015)

work page 2015

[5] [5]

A. G. Shtukenberg, Q. Zhu, D. J. Carter, L. Vogt, J. Hoja, E. Schneider, H. Song, B. Pokroy, I. Polishchuk, A. Tkatchenko, et al.,Powderdiffractionandcrystalstructurepredictionidentify four new coumarin polymorphs, Chemical Science8, 4926 (2017)

work page 2017

[6] [6]

D. E. Braun, M. Ardid-Candel, E. D’Oria, P. G. Karamertzanis, J.-B.Arlin, A.J.Florence, A.G.Jones,andS.L.Price,Racemic naproxen: a multidisciplinary structural and thermodynamic comparison with the enantiopure form, Crystal growth & design 11, 5659 (2011)

work page 2011

[7] [7]

G. M. Day, Current approaches to predicting molecular organic crystal structures, Crystallography Reviews17, 3 (2011)

work page 2011

[8] [8]

Nyman and S

J. Nyman and S. M. Reutzel-Edens, Crystal structure prediction is changing from basic science to applied technology, Faraday Discussions 211, 459 (2018)

work page 2018