Using a different Optimization.jl compatible solver
In this example we perform the simplistic optimisation the bond length of a Hydrogen molecule using a trust region quasi-Newton method from NLopt.
We create a calculator employing the density-functional toolkit to compute energies and forces at using the LDA density functional.
using DFTK
using PseudoPotentialData
pseudopotentials = PseudoFamily("dojo.nc.sr.lda.v0_4_1.oncvpsp3.standard.upf")
model_kwargs = (; functionals=LDA(), pseudopotentials)
basis_kwargs = (; kgrid=(1, 1, 1), Ecut=20.0)
calc = DFTKCalculator(; model_kwargs, basis_kwargs)DFTKCalculator(functionals=Xc(lda_x, lda_c_pw), pseudopotentials=PseudoFamily("dojo.nc.sr.lda.v0_4_1.standard.upf"), Ecut=20.0, kgrid=(1, 1, 1))and we build the hydrogen molecular system, where we attach pseudopotential information for DFTK:
using AtomsBuilder
using Unitful
using UnitfulAtomic
bounding_box = [[10.0, 0.0, 0.0], [0.0, 10.0, 0.0], [0.0, 0.0, 10.0]]u"Å"
system = periodic_system([:H => [0, 0, 1.]u"bohr",
:H => [0, 0, 3.]u"bohr"],
bounding_box)
nothingWe now run GeometryOptimization.minimize_energy!, but notably pass the NLopt.LD_TNEWTON solver from NLopt as the third argument to employ this solver. Extra keyword argument to NLopt can be added, e.g. here the maxevel=100, which limits the solver to 100 function evaluations:
using GeometryOptimization
using OptimizationNLopt
solver = NLopt.LD_TNEWTON()
results = minimize_energy!(system, calc, solver;
tol_forces=1e-4u"eV/Å", verbosity=1,
maxeval=100)
nothingGeometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.105713124397 │ 0.04 │ 0.0820222 │ │ │ 21.5s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.105713128065 │ -8.44 │ 0.0820221 │ │ │ 24.8s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.105713126263 │ + -8.74 │ 0.0820222 │ │ │ 28.0s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.105713125008 │ + -8.90 │ 0.0820222 │ │ │ 31.1s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.104146106594 │ + -2.80 │ 0.0827083 │ │ │ 40.0s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ +258.8619812484 │ + 2.41 │ 68417.0 │ │ │ 48.4s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -0.847587721778 │ 2.41 │ 1.33866 │ │ │ 56.3s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.132167594536 │ -0.55 │ 0.0521148 │ │ │ 65.1s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.132167596866 │ -8.63 │ 0.0521148 │ │ │ 67.8s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135213860862 │ -2.52 │ 0.00717939 │ │ │ 74.2s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135213861183 │ -9.49 │ 0.00717937 │ │ │ 77.2s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288094661 │ -4.13 │ 0.000359232 │ │ │ 83.6s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288094677 │ -10.79 │ 0.000359217 │ │ │ 86.6s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288288806 │ -6.71 │ 2.54941e-6 │ │ │ 92.0s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288288807 │ -12.87 │ 2.53435e-6 │ │ │ 94.8s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288288807 │ + -12.92 │ 2.54848e-6 │ │ │ 97.9s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288288806 │ + -14.00 │ 2.55158e-6 │ │ │ 101s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288288807 │ -13.05 │ 2.53881e-6 │ │ │ 104s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288288807 │ + -13.14 │ 2.54816e-6 │ │ │ 107s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288288854 │ -10.32 │ 2.56088e-6 │ │ │ 111s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288302362 │ -7.87 │ 6.4624e-6 │ │ │ 117s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135204765100 │ + -4.08 │ 0.00759705 │ │ │ 123s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135287958776 │ -4.08 │ 0.00099143 │ │ │ 132s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288530754 │ -6.24 │ 0.000364815 │ │ │ 138s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288530770 │ -10.79 │ 0.000364799 │ │ │ 141s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288530755 │ + -10.81 │ 0.000364813 │ │ │ 144s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288731668 │ -6.70 │ 2.40073e-5 │ │ │ 149s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288731670 │ -11.98 │ 2.3998e-5 │ │ │ 152s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288731669 │ + -12.39 │ 2.40068e-5 │ │ │ 155s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288731669 │ + -12.88 │ 2.40087e-5 │ │ │ 158s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288731668 │ + -12.00 │ 2.40148e-5 │ │ │ 161s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288731668 │ -12.34 │ 2.40079e-5 │ │ │ 164s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288830762 │ -7.00 │ 0.000157053 │ │ │ 170s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288830769 │ -11.15 │ 0.000157038 │ │ │ 173s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288830762 │ + -11.20 │ 0.000157053 │ │ │ 176s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288830763 │ -12.40 │ 0.000157052 │ │ │ 179s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135289044494 │ -6.67 │ 2.143e-5 │ │ │ 185s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135289044495 │ -12.01 │ 2.14272e-5 │ │ │ 188s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135289044495 │ + -13.04 │ 2.14295e-5 │ │ │ 191s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135289048542 │ -8.39 │ 5.02077e-5 │ │ │ 196s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135289048544 │ -11.65 │ 5.0195e-5 │ │ │ 200s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135289048543 │ + -11.86 │ 5.02074e-5 │ │ │ 203s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135289048542 │ + -12.05 │ 5.0213e-5 │ │ │ 205s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135289048541 │ + -12.15 │ 5.02123e-5 │ │ │ 209s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135289048542 │ -12.16 │ 5.02077e-5 │ │ │ 212s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135289154878 │ -6.97 │ 6.62141e-5 │ │ │ 217s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135288399573 │ + -6.12 │ 0.00163011 │ │ │ 223s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135290098123 │ -5.77 │ 0.000643492 │ │ │ 231s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135290098152 │ -10.54 │ 0.000643476 │ │ │ 234s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135290098127 │ + -10.61 │ 0.000643489 │ │ │ 237s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135290719114 │ -6.21 │ 8.00638e-5 │ │ │ 242s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135290719117 │ -11.45 │ 8.00522e-5 │ │ │ 245s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135290719115 │ + -11.67 │ 8.00643e-5 │ │ │ 248s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135290719116 │ -12.03 │ 8.00605e-5 │ │ │ 251s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135290719112 │ + -11.36 │ 8.00754e-5 │ │ │ 254s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135290719115 │ -11.48 │ 8.00632e-5 │ │ │ 257s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135291072855 │ -6.45 │ 0.000126333 │ │ │ 263s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135291072861 │ -11.25 │ 0.00012632 │ │ │ 266s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135291072857 │ + -11.39 │ 0.000126333 │ │ │ 269s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135291072858 │ -11.92 │ 0.00012633 │ │ │ 272s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135291072857 │ + -11.94 │ 0.000126332 │ │ │ 275s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135291357306 │ -6.55 │ 0.000139297 │ │ │ 281s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135291357312 │ -11.21 │ 0.000139283 │ │ │ 283s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135291357307 │ + -11.34 │ 0.000139296 │ │ │ 287s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135291357304 │ + -11.49 │ 0.000139303 │ │ │ 290s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135291357307 │ -11.62 │ 0.000139297 │ │ │ 293s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135291764386 │ -6.39 │ 7.26115e-5 │ │ │ 298s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135293042876 │ -5.89 │ 0.00113549 │ │ │ 305s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135293042927 │ -10.29 │ 0.00113547 │ │ │ 308s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135293042879 │ + -10.32 │ 0.00113549 │ │ │ 311s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135295004040 │ -5.71 │ 5.1921e-5 │ │ │ 317s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135295004043 │ -11.64 │ 5.19199e-5 │ │ │ 320s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135295004042 │ + -12.98 │ 5.19213e-5 │ │ │ 323s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135295004042 │ -13.99 │ 5.19212e-5 │ │ │ 326s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135298514163 │ -5.45 │ 0.000305787 │ │ │ 332s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135298514177 │ -10.86 │ 0.000305773 │ │ │ 335s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135298514166 │ + -10.95 │ 0.000305787 │ │ │ 338s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135298514164 │ + -11.77 │ 0.000305789 │ │ │ 341s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135298514166 │ -11.83 │ 0.000305787 │ │ │ 344s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135298432620 │ + -7.09 │ 0.00146554 │ │ │ 350s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135299216538 │ -6.11 │ 0.000843454 │ │ │ 357s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135299216575 │ -10.42 │ 0.00084344 │ │ │ 360s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135299216540 │ + -10.45 │ 0.000843454 │ │ │ 363s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135299216541 │ -11.96 │ 0.000843454 │ │ │ 366s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135304107523 │ -5.31 │ 0.000434431 │ │ │ 373s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135304107543 │ -10.71 │ 0.000434415 │ │ │ 376s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135304107526 │ + -10.78 │ 0.000434429 │ │ │ 379s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135304397104 │ -6.54 │ 5.34663e-5 │ │ │ 384s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135304397106 │ -11.62 │ 5.34659e-5 │ │ │ 387s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135304397106 │ + -13.36 │ 5.34664e-5 │ │ │ 390s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135304397106 │ -13.41 │ 5.34661e-5 │ │ │ 393s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135304397106 │ + -13.79 │ 5.34662e-5 │ │ │ 396s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.134936212307 │ + -3.43 │ 0.0150655 │ │ │ 402s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135304654642 │ -3.43 │ 0.0014815 │ │ │ 411s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135304654708 │ -10.18 │ 0.00148149 │ │ │ 414s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135304654645 │ + -10.20 │ 0.0014815 │ │ │ 417s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135307990975 │ -5.48 │ 8.04448e-5 │ │ │ 423s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135307990978 │ -11.44 │ 8.0435e-5 │ │ │ 426s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135307990977 │ + -11.82 │ 8.04444e-5 │ │ │ 429s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135307990976 │ + -12.01 │ 8.04481e-5 │ │ │ 432s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135307990976 │ -12.98 │ 8.04459e-5 │ │ │ 435s │
Geometry optimisation convergence (in atomic units)
┌─────┬─────────────────┬───────────┬─────────────┬─────────────┬──────────┬────────┐
│ n │ Energy │ log10(ΔE) │ max(Force) │ max(Virial) │ Pressure │ Δtime │
├─────┼─────────────────┼───────────┼─────────────┼─────────────┼──────────┼────────┤
│ 0 │ -1.135308699085 │ -6.15 │ 0.000224056 │ │ │ 441s │
┌ Warning: Geometry optimisation not converged.
└ @ GeometryOptimization ~/work/GeometryOptimization.jl/GeometryOptimization.jl/src/optimization.jl:156The final hydrogen bond length is:
using LinearAlgebra
norm(position(results.system[1]) - position(results.system[2]))1.4531922767206555 a₀