The Forth Bridge. Widely believed to be painted continuously since 1890.

This list is kept against the development log. Items done since August 2026 stay here for a while, so that anyone who knew the old list can see what moved. Items in hand carry a name.

Done since August 2026

(By Mike Towler working much harder than usual, except where noted)

New methods and measurements

• DMC that stops itself: when the error bar reaches a requested target, or at the smallest error bar a given time budget can buy.
• popbias: the population-control bias of a finished DMC run, and the energy corrected for it, with error bars.
• EQUIL_DIAGNOSTICS: after a DMC run, where two standard rules would have ended equilibration, and how far the answer would then have moved.
• Walker lineages: DMC reports how many of the walkers it began its statistics with still have descendants, and since when the whole population has descended from one of them, if it has.
• LIMDMC 7 and 8: the reweighting of Anderson, Per and Umrigar, which its authors found reduces time-step errors at large time steps.
• VMC_STEP_PROFILE: smaller VMC steps near all-electron nuclei, making VMC of all-electron atoms and small molecules heavier than helium two to five times cheaper (Konjkov/MDT).
• Densities with error bars: a zero-variance estimator (ZV_DENSITY) for molecules and crystals, error bars on the ordinary density, and expval2cube to write them as cube files for VESTA and other viewers.
• The exchange-correlation hole and the on-top pair density of a molecule, with error bars (HOLD_PARTICLE), and hole2html to drag the electron through the hole.
• The localisation tensor of a correlated state, which tells an insulator from a metal.
• Reptation Monte Carlo repaired: pure expectation values with error bars, every process’s reptile combined, and reptiles that no longer freeze near a node.
• Reptation Monte Carlo with DMC’s time-step error, where its pure estimates for an all-electron atom had been several per cent off at usable time steps.
• Reptation Monte Carlo no longer marked as ‘experimental’: moved to intermediate, with its non-existent chapter of the manual written from scratch.
• Light nuclei, such as protons and muons, treated as quantum particles alongside the electrons in DMC.
• The exact DMC energy of the muonic molecular ion ppmu⁺ at zero time step, with a Morse tail in the Jastrow term for the heavy pair (UHARM).
• The probability of finding a given number of electrons in a region of space (DOMAIN_PROB).
• A uniform electric field in a finite system (ELECTRIC_FIELD), for energies, dipole moments and polarisabilities in the field.
• DMC estimates of the energy’s derivatives with respect to determinant coefficients (DMC_NODE_DERIV), a first estimate of the nodal gradient.
• Nodes made visible: nodeview draws a wave function’s nodal surface in a web page, nodecells counts its nodal cells and tests Ceperley’s tiling property, and wfnscope lets you drag the electrons and watch the node, the local energy and the conditional density move.
• nodeview draws the superposition of two wave functions evolving in time: its node is the lines where the two nodes cross, becoming a whole surface twice in each period.
• eosfit : equations of state fitted to QMC energies, with error bars on the equilibrium volume, the bulk modulus and its derivative.
• A Tersoff bond-order potential among the manual interactions (NDD).
• dmcsound: listen to a DMC run. It turns the histories of one or more runs into a WAV file in which every excited state that a fit to the run’s energy resolves is a tone, fading at the rate its energy gap sets, so that you hear the walkers forget their starting point; where no state is resolved it plays the measured equilibration. Each sound comes with a list of what in it is measured and what is convention.
• Pictures from a run: WALKER_FILM films the walkers of a DMC or VMC run, frame by frame, with the numbers that tie each walker to its parent, and the new utility walkerfilm makes of the film a trajectory for OVITO or Blender or a web page that plays it, in place of the old movie, which is withdrawn with movie2avi.
• webqmc: CASINO’s VMC and DMC in a web browser, one file with nothing to install that opens on a phone, drawing the walkers as the energy accumulates (ghosts at Hallowe’en).
• pilotdrive and RUNTYPE pilot_eval: de Broglie-Bohm trajectories in superpositions of CASINO wave functions, and the relaxation of their ensembles to quantum equilibrium.
• pilotdrive: trajectories under the pilot-wave law of the spin-1/2 state, with the spin term of the current; runs checkpoint and resume, and are faster, 3.6 times for a diamond cell with a pseudopotential.
• An exact DMC step for an electron near a node (DMC_NEAR_NODE bessel), which with branching off samples the trial density exactly at any time step.

Wave-function optimisation

• Energy minimisation now tries the full step it computes, which its search had replaced by a meaningless one, usually refused.
• Analytical Jastrow derivatives corrected: several had been wrong, and cycles with fixed cutoffs had scaled them by unrelated factors for over a decade.
• Variance and MAD minimisation no longer stop early on false or singular convergence after a rejected step.
• Analytical derivatives for backflow parameters and cutoffs, determinant coefficients and f-term cutoffs: optimising backflow in finite systems is about six to seven times faster.
• Reweighted variance minimisation much less prone to run away (VM_VARW_MAX), and a warning when an optimisation cycle leaves the energy or variance far above the best so far.
• Variance and MAD minimisation of backflow with a pseudopotential recompute the non-local energy by default, and no longer report an improvement while the wave function gets worse.
• CHECK_ADERIVS checks the derivatives the optimisers actually use.
• Optimisation of Slater-type orbitals repaired: in some inputs it changed nothing, and it tied unrelated coefficients together and left redundant rotations free.
• Geminals: open shells supported and optimisable, constraints applied, manual and help current.
• The Jastrow H term’s cutoff derivative analytical, no longer by finite differences.
• Variance and MAD minimisation no longer stop at a trial step that overflows a local energy, and energy minimisation no longer takes a far-wrong first step from a parameter given outside its limits.
• Variance minimisation with many configurations needs a third less memory, e.g. 35 GiB in place of 52 for 20 million configurations and 112 parameters.
• Optimisation with a pseudopotential under OPT_FIXNL F takes complete Jastrow and backflow derivatives, where the change in the non-local energy had been left out since 2012.
• Energy minimisation no longer steps along a spurious eigenvector of an ill-conditioned problem, which could make its line search fail.
• Variance and energy minimisation of a complex wave function with non-local pseudopotentials use the whole imaginary part of the local energy, not its kinetic part alone.
• Optimisation of a general Jastrow term whose electron-electron functions are odd under exchange, such as the W term, no longer steps from wrong derivatives.
• Optimisation of the Gaussian orbitals of an open-shell molecule no longer crashes when one set of orbitals serves both spins, as an ROHF radical’s does, nor optimises coefficients that were never set.

Backflow and Jastrow factors

• The three-electron backflow term omega keeps the wave function antisymmetric again, is sampled exactly, and has all its spin triplets optimised.
• The backflow Phi term’s cusp conditions, applied to the wrong coefficients since version 2.0 whenever its sets had different orders, are now right.
• Backflow works with periodic all-electron Gaussian orbitals and the cusp correction, where energies could be wrong by tens of hartrees a cell.
• A parameters.casl Jastrow term that excludes a nucleus no longer counts it: DMC_METHOD 2 energies of water had been 0.52 Ha too low.
• DMC with a geminal wave function and backflow now right; it had been wrong since the two were first combined in 2013.
• Jastrow terms no longer impose a cusp the orbitals already carry, which had made the local energy diverge or the variance many times larger.
• A general Jastrow term with a single channel no longer counts the electron pairs its rules exclude; a Jastrow factor holding one is worth optimising again.

Getting the right answer

• All-electron DMC_METHOD 2 branches with its own factor, not the T-move one it had used since 2018.
• Parallel DMC_METHOD 2 runs no longer use stale or uninitialised Jastrow values on walkers moved between processes, a bias of about 0.1 mHa per atom in silicon.
• DMC walkers are reweighted correctly after a change of wave function, with space warping in DMC molecular dynamics, and at a change of twist.
• Gaussian shells are screened on what they can contribute, removing errors of several mHa per cell in crystals such as nickel oxide.
• Cusp corrections of the virtual orbitals of multideterminant expansions are right; wrong ones could change the fixed-node energy markedly.
• Compressed multideterminant expansions with zero coefficients, as at the end of large CIPSI expansions, no longer change the wave function.
• Gaussian f and g functions: wrong g-function gradients under screening and wrong f-function values under periodic backflow corrected.
• VMC jumps between ions no longer bias the energy.
• NaN energies stop the run in optimised builds, where the tests had been compiled away and NaN results reported.
• POPSTATS no longer prints a failed estimate as a number, and warns when a run is too short for it.
• Slater-type orbitals with n = 5 normalised correctly; n up to 10 supported.
• Forces: planarity test and per-atom charges fixed, periodic forces refused, all-electron forces flagged as unreliable.
• Positrons: each ion type has its own positron pseudopotential, and a positron in a molecule sees the ions correctly.
• Ewald sums correct on nearly symmetric cells in 2D and 3D; the 1D Coulomb sum accurate to 10^-10 hartree per pair.
• A core-polarisation potential’s l-dependent term no longer cut off at the pseudopotential’s radius.
• Momentum densities and the condensate fraction weighted correctly in DMC.
• Post-weighted expectation values in parallel DMC use every process’s walkers, not the master’s alone.
• Expectation values under an alternative VMC sampling distribution carry the trial wave function’s weights.
• Twist-averaged DMC keeps the earlier twists’ expectation values through a catastrophe reset in a later one.
• Plane-wave wave functions no longer give a wrong energy for a k point given as a grid point’s time reverse, or for degenerate bands with COMPLEX_WF F.
• Error bars on an atom’s spherical and intracule densities, which had been too small.
• Finite-size corrections with the Jastrow factor of parameters.casl: the interaction-energy correction is applied, where it had been refused.
• Reptation Monte Carlo’s error bars no longer understated.
• Smaller faults put right: SPLOT‘s origin, MAD minimisation’s filter, the variance and the correlation period of VMC, and the orbital-parameter derivatives of one particle.
• GPU builds (OpenACC): the Jastrow factor’s p term right in one- and two-dimensional periodic systems, with G vectors in any order, and in every cycle of an optimisation, not the first alone.
• The general-purpose cusp correction no longer lets the local energy diverge at the nuclei of a spin-polarised molecule with restricted orbitals.
• A run beside an existing expval.data no longer stops or crashes over its POPULATION block, present or absent, and PAIR_CORR refuses every other expectation value, which its fixed particle would spoil.
• CASINO’s own BLAS library no longer reports a finite norm for a vector holding a NaN in an optimised build.
• The test suite compares every file a test writes with the library’s, the history files, config.out and expval.data included, not only its out file.
• The test suite has a deck for each mode that no deck reached before, and a census of what every deck exercises, which shows what is still not reached.
• Numerical atomic orbitals read from awfn.data: one f orbital was used at half its proper size and four g orbitals with the wrong sign; now right.
• A continued run that would add to an accumulation of other quantities than its own, after a change of input or of version, is refused with the difference named, where before it mixed them without a word; config.in now records what each of its accumulations is of.

Speed

• Faster Ewald sums: the reciprocal-space part of the 3D electron-electron sum now grows linearly with the number of particles rather than as its square, and slabs share the work between ions in a plane.
• Gaussian orbitals evaluated faster: DMC 17 to 33 per cent faster for a cluster of 64 water molecules.
• Slater-type orbitals screened (STOTOL) and evaluated faster: VMC of C₃₂H₆₆ 2.5 times as fast from the screening, and 1.7 times faster again from the evaluation.
• Faster T-moves, multideterminant DMC and Jastrow F term.
• DMC can share walkers by measured core speed, for desktops with fast and slow cores (REDIST_SPEED_WEIGHT).
• The set-up of reciprocal-space expectation values on a long cell: from about three weeks to fifteen seconds for a hydrogen chain of 53 million G vectors, quite definitely the largest speed-up in the history of CASINO’s development.
• Faster wave-function optimisation, up to 2.7 times for variance minimisation with a general Jastrow factor.
• Geminal wave functions evaluate, where the orbital type allows, only the orbitals their pairing matrices use: DMC with a silane geminal using 4 of its 116 orbitals runs about 2.7 times faster.
• Reptation Monte Carlo about a quarter faster a move, and lighter in memory between its stages.
• The backflow Phi and omega terms do a quarter to a half less arithmetic, saving 6 to 9 per cent of the work of a DMC step with a Phi term.

Other codes

• Cartesian d, f and g Gaussian shells read natively, as in 6-31G* and other Pople basis sets; converters say what they can produce.
• CASINO checks at start-up that a molecule’s Gaussian orbitals are orthonormal, catching converters that got the basis wrong.
• MOLDEN route completed and tested for PSI4, DALTON, Molpro, ORCA, NWChem, OpenMolcas and TURBOMOLE.
• New converters: trexio2qmc for TREXIO files, Quantum Package’s CIPSI expansions included; qbox2pwfn and casino2qso for QBOX; crysgen23 for CRYSTAL23.
• GAMESS converter remade and tested with GAMESS 2024, f and g functions included.
• Files from GAMESS’s own CASINO writer, whose d functions and open-shell orbitals are wrong, are refused in favour of gamess2qmc.
• Q-Chem through molden2qmc: its MOLDEN files in spherical or Cartesian d, f and g functions, its restricted open-shell orbitals, and its own output read directly.
• AMS/ADF: the adf2stowf converter distributed with CASINO.
• CASTEP: runcastep runs CASTEP and converts its orbitals in one command, and twistav_castep resumes an interrupted calculation.
• Quantum ESPRESSO 7.6: runpwscf and runqmcmd read its files.
• CRYSTAL: runcrystal brought up to CRYSTAL23 with exit statuses; g functions from CRYSTAL in molecules; opt_crystal repaired, with new -kappa and -stall options; a manual appendix on choosing and optimising CRYSTAL basis sets for QMC.
• 2DHF numerical orbitals for diatomic molecules usable again, from x2dhf 3.0.
• ABINIT: how to run it for CASINO and give it a CASINO pseudopotential.
• LAPW orbitals: all-electron trial wave functions from WIEN2k, exciting, FLEUR and Elk, core states included, read directly for VMC, DMC and RMC; spin-polarised WIEN2k too, and writers for exciting, FLEUR and Elk on request.
• MCEXX writes CASINO’s pwfn.data itself (l_casino; the writer available from MDT on request with permission of the authors), with an example of silicon.
• nwchem_pp_convert correct for ccECP and for NWChem’s own libraries.
• FHI-aims: its all-electron numeric atom-centred orbitals used for molecules and crystals (ATOM_BASIS_TYPE nao, the file written by aims2qmc from FHI-aims’ TREXIO or band output), with examples and its chapter in the manual.
• Numeric atom-centred orbitals (nao) now take backflow, pseudopotentials and ORB_NORM, and the new awfn2nao converts an atom tabulated for numerical to them, determinant expansion included.
• runpwscf and runcastep refuse ultrasoft and PAW pseudopotentials, and those CASTEP generates itself, whose orbitals CASINO cannot use, instead of running on to wrong energies.

Running CASINO

• CASINO, runqmc and autotest exit with a failure status when a run fails, and parallel runs stop rather than hang when an error reaches only some processes.
• Wave-function optimisation and DMC in a single run.
• Many malformed or unsupported inputs, over-long input lines among them, now refused with a message instead of being silently misread.
• runqmc no longer appends a new run to an earlier run’s output (--keep-old keeps the earlier files).
• Builds with Intel ifx, an optional MPI-module build with argument checking, and OpenMP shared-memory groups fixed.
• Regression suite extended: RMC, geminals, backflow, analytical derivatives, optimisation with fixed cutoffs and more.
• reblock: unconverged quantities and short-run POPSTATS warnings.
• mpc.data generated automatically when a crystal’s MPC interaction needs it; WRITE_GEOMETRY writes the atoms for viewing.
• Installing without a terminal (./install --arch=<CASINO_ARCH> --compile --save), and casinocheck to test an installation.
• A built tree runs when moved or renamed, so that versions can be kept side by side, and ~/.bashrc.casino no longer touches the user’s other shell variables.
• Reptation Monte Carlo stops at a time limit and continues, as VMC and DMC do.
• Reptation Monte Carlo writes rmc.hist, so that reblock, graphdmc, plot_hist and endmc read its runs as they read VMC’s and DMC’s.
• runqmc clears a dead run’s lock itself, and a vmc_dmc run stopped at a time limit can be continued.
• runqmc --continue no longer restarts from the beginning a run that cannot be continued, nor loses a general Jastrow factor’s optimised parameters.
• runqmc no longer stops a twist-averaged VMC chain after its first twist.
• A parallel run that compresses the multideterminant expansion of its orbital file no longer hangs at the start, and a DMC run killed while writing its expectation values keeps those of the blocks before.
• CASINO compiles again with gfortran 8.4 and later.
• New utility dftqmc: from a PySCF calculation on a small molecule to a DMC energy with an honest error bar, in one command.

Documentation

• The manual completely rewritten for new users, with new chapters on QMC and CASINO and on reptation Monte Carlo, and the chapter on DMC and other stuff completely redone.
• Every input keyword’s help rewritten; the manual’s keyword reference is now generated from the code, so the two agree.
• New guides for developers: doc/PROGRAMMING_GUIDE, in place of 2001’s, and doc/GIT_GUIDE, in place of 2008’s.
• An appendix to the manual of the faults of earlier versions of CASINO, each with the versions that brought it in and put it right. This is embarrassingly extensive.
• The FAQ remade from the forum’s questions, in fourteen sections from installing CASINO to pilot waves.
• A feature map in the manual: for every kind of orbital and every route from another program, what works, what is refused and what is untested.
• Running CASINO on a GPU: the build step by step, where it gains and where it loses, and how to rent one by the hour.
• The neural-network Jastrow factor (NNJAS) described in the manual as the code computes it.
• Examples setting reptation Monte Carlo against DMC for helium and lithium hydride, and a section of the manual on when RMC is worth its cost.
• The examples’ READMEs say how to run them, and the recorded results are those the present inputs give.
• Correlated sampling with vmc_sampling 'stdtest' documented.

LOUIS

• LOUIS 1.3: a user manual, an ensemble mode with Nelson’s stochastic dynamics, and corrected 3D densities.
• LOUIS’s coarse-grained H-function right in a density run in which a backward trajectory fails.

To Do

What we would like done next, by subject, the larger jobs first in each. An item carries the name of whoever is on it, so that two people do not do the same thing; an item without a name is free. If you want one, say so on the forum or write to Mike Towler, and it is yours. The developers’ page says how a change reaches the code, and the development log shows what became of the last list.

Open: new methods and measurements

Things CASINO cannot yet measure, or cannot yet do, that people ask for.

• Correlated sampling in DMC. VMC has it (vmc_sampling stdtest); DMC does not, and it is what an adsorption energy, a reaction barrier or the small displacements of a geometry optimisation need, since two independent runs resolve a difference of a millihartree only at a cost that grows as the square of the ratio of the local energy’s spread to that difference.
• Forces for plane-wave and blip orbitals and in periodic systems, where they are refused today: the derivative of the ion-ion Ewald sum, the change of the orbitals with the nuclei (by space warping, or from the generating code), core-polarisation potentials and f and g pseudopotential channels. Then a geometry optimiser that knows its forces carry error bars.
• Excited states on the same footing as the ground state: state-specific optimisation (variance matching, and Zhao and Neuscamman’s variational principle for excited states) beside the promotion method, and the optical gaps of solids with twist averaging, as a documented workflow.
• Light nuclei in the rest of CASINO: DMC molecular dynamics, reptation, and a hydrogen atom moving among fixed nuclei in a metal.
• Nuclear quantum effects by averaging: electronic energies averaged over nuclear configurations drawn from a vibrational wave function or a path-integral run, with correlated sampling to take the differences cheaply.
• Neural-network trial wave functions: a plug-in evaluator, value, gradient and Laplacian at a configuration, called at every move, through which DMC runs on a FermiNet or Psiformer wave function or on a learned Jastrow factor. DMC on such functions exists elsewhere; the point is CASINO’s T-moves, twist averaging and expectation values on them. CASINO already has a small neural-network Jastrow term of its own, NNJAS (experimental); this item is for wave functions trained elsewhere.
• Spin-orbit coupling: spinor determinants from relativistic CRYSTAL and DIRAC calculations, with spin-orbit pseudopotentials, in VMC and DMC; complex then spinor Gaussian orbitals, then readers, projectors and spin sampling. Mike Towler, with the CRYSTAL developers.
• LAPW orbitals beyond their first release: pseudopotentials, backflow, the MPC interaction, localised blip orbitals between the spheres and relativistic WIEN2k calculations, all refused or unconverted today.
• Bader charges and the topology of a QMC density with error bars, from the cube files CASINO now writes, as a utility rather than a recipe.
• Electrons in separated subsystems as distinct particles, for bilayers and double quantum wells where antisymmetry between the layers costs much and buys nothing.
• A criterion that reports an unbound trial wave function for a small complex instead of a slowly drifting energy: an MDT test of the energy’s fall missed every unbound case; try how far the particles wander.
• The population-control bias inside the run: the first-order estimate that popbias now makes afterwards from a run’s own history, shown in the automatic stopping’s and EQUIL_DIAGNOSTICS‘ reports, so that a run stopped at its target error bar states the bias it carries.
• Embedding: a QMC region inside a classical or DFT environment, beginning from the external potentials CASINO already reads and going on to an embedding potential derived from a DFT density.

Open: wave functions

Jastrow factors, backflow, pairing and determinants: the forms CASINO can optimise and run.

• The general Jastrow factor and general backflow as the documented defaults, their cost against the standard forms measured on real systems before the switch.
• Open-shell geminals with backflow, and with a complex wave function, which cannot be combined today; then Pfaffians (PSI_S pfaffian is refused today), which hold geminals and determinants as special cases and allow triplet pairing.
• Hybrid orbitals, Gaussians plus blips, for a molecule on a surface or a defect in a crystal, where one basis serves neither part well.
• A library of starting Jastrow terms by element, the electron-nucleus and electron-electron-nucleus terms that carry from one system to another, so that an optimisation starts near its answer.
• The plane-wave p term of the Jastrow factor re-represented in blips, so that it costs what a blip orbital costs.

Open: wave-function optimisation

Getting the parameters right, and knowing when they are.

• An audit of energy minimisation’s convergence, with the gates a run must pass before its parameters are trusted: the condition of its matrices, the quality of its line-search fit and the stability of its lowest eigenvector, reported in the output rather than left to the user’s eye.
• Optimisers for very many parameters at once: a descent method, or the sample-space forms of the linear method and stochastic reconfiguration whose cost grows with the number of samples rather than of parameters, for the general Jastrow, backflow and the coefficients of a selected-CI expansion of a hundred thousand determinants together.
• An optimisation that stops when it has converged, as a DMC run now stops at its target error bar, instead of running the number of cycles it was given.
• Analytical derivatives for what still uses finite differences, the general Jastrow’s cutoffs and basis among them.

Open: getting the right answer

Correctness, and the checks that prove it.

• All-electron forces: an estimator with an honest error bar, or a refusal. Today they print, flagged as unreliable, because the bare estimator’s variance is infinite; the renormalised estimator of Assaraf and Caffarel and the filtered one of Chiesa, Ceperley and Zhang are the known remedies.
• The DMC time step: Anderson and Umrigar’s 2021 form of the T-move (their reweighting is in, as LIMDMC 7 and 8), measured on pseudopotential systems; and a time-step extrapolation check in the analysis tools, since two extrapolations of the same data can disagree by millihartrees.
• expval.data does not record the atomic positions its expectation values were accumulated at, so a continued run with the ions moved adds to a density or a pair-correlation function without a word.
• A reference set that catches a wrong answer and not only a changed one: the regression suite compares a build with the outputs an earlier build wrote. Add benchmark energies with published values, atoms, molecules and solids, the hydrogen chain of Motta and co-workers (2017) among them, rerun at every release.
• Every expectation value checked under DMC weighting, twist averaging and VMC’s alternative sampling against an independent estimate: several weightings were found wrong in 2026 by reading rather than by any test.
• The coincidence limits, second part: the general Jastrow factor’s kinetic energy where electrons meet, and the local energy on the axis of numerical dimer orbitals, which give a gradient of exactly zero there.

Open: speed and scale

A run of measured speed-ups has shipped since August; these are the levers left, each to ship with its measured timing.

• CASINO on GPUs beyond the present OpenACC build: the Gaussian and blip evaluators, the Jastrow factor and the Ewald sums on the device, measured on the standard test systems against the processor at equal cost.
• Large multideterminant expansions from selected CI: the work done once per determinant shared over distinct spin strings, for real wave functions first, projected at several times faster at a hundred thousand determinants.
• Parallelisation over twist angles, each twist its own DMC population, so that a twist-averaged run fills a whole machine at once and the twists finish together; and no re-equilibration after a twist change that leaves the occupancy unchanged.
• Scaling measured again: MDT ran CASINO on many hundreds of thousands of cores in 2011 and nothing has been measured above a few hundred ranks since. A scaling study on a national machine with today’s code, and whatever it finds, fixed. (Having once had simultaneous accounts on four of the top-ten most powerful machines in the world, MDT has been brought low and is currently trying to prove that Vallico Academy is a real institute so he can apply for time on more than 64 cores. Embarrassing. If anyone can help let him know.)
• Fault tolerance for very large runs: a run that survives the loss of a node and continues from its last configurations.
• Configuration files written and read by MPI-IO, so that a checkpoint of a million walkers is not a serial bottleneck.
• More of the large arrays in shared memory within a node, as the orbital coefficients already are.
• OpenMP that pays: on a 64-core workstation two threads a rank are slower than two ranks, so the threading buys memory and nothing else.
• The periodic Gaussian contraction at complex k points masks inside its loops; masking once afterwards would cut a DMC step’s work by 3 to 12 per cent, moving the energies’ last digits. Mike Towler.
• Numerical orbitals are all evaluated and then masked, where only the masked ones are needed; and the periodic cusp correction zeroes large arrays on every call.
• FFT grids for mpc.data that are not powers of two, and FFTW as an alternative library.
• The cost per move at a thousand electrons profiled and attacked, as was done this year at a hundred, the localised orbitals and the Ewald sums first.

Open: external codes

Routes for trial wave functions into CASINO; the interfaces page says what works today. A new Gaussian-basis route is usually a MOLDEN or TREXIO file, which CASINO’s converters read already; otherwise a route needs a writer of pwfn.data, gwfn.data, bwfn.data or stowfn.data, and the interface guide says how. We help.

• Converters from OpenMX and DMol3 to naowfn.data, once their orbitals pass a DMC test against blips. SIESTA’s failed it; the manual says why.
• Converters for VASP (with norm-conserving potentials) and OpenMX. Aleksey Baglov, exploring.
• CFOUR: the corrected reorderdf.f that its MOLDEN output needs is distributed with CASINO; the route is untested.
• GAUSSIAN: the converter is twenty-five years old, spherical functions only, and unmaintained; a current route is wanted. Mike Towler.
• TREXIO as the common currency: CASINO reads orbitals and determinants from it today; next, pseudopotentials and Jastrow parameters, and writing CASINO’s own trial function to it, so that one wave function can go to CASINO, QMCPACK, TurboRVB, CHAMP or an AFQMC code and the methods be compared on equal terms.
• CP2K, Octopus and BigDFT: untried. CP2K writes MOLDEN files for molecules, so molden2qmc first; the real-space and wavelet codes need a conversion of their orbitals to blips.
• Blip converters inside CASTEP and ABINIT, as Quantum ESPRESSO has, so that the DFT code writes blips itself and no pwfn.data of tens of gigabytes is written and read.
• Consistent treatment of unoccupied orbitals across the plane-wave interfaces, so that an excitation means the same thing whatever code made the orbitals.
• Point-by-point verification carried to every route: orbitals and densities compared with the generator’s own at sample points, as the MOLDEN routes were in 2026, so that every entry on the interfaces page has a measured verdict.
• Gaussian basis sets built for QMC rather than borrowed from quantum chemistry: the successor to billy and opt_crystal (not the one embedded in CRYSTAL23!), and the study of what an optimised basis does to Hartree-Fock, DFT and DMC energies. Mike Towler.
• Pseudopotentials: spin-orbit potentials for the spinor route; every potential of the library in every code’s format from the one converter CASINO ships; and the library’s element pages brought up to date with the ccECP set added in October 2026.

Open: running CASINO

Installing it, driving it, and reading what it says.

• Defaults that let a DFT user get a trustworthy DMC number: a Jastrow factor chosen and optimised, the time step and population set from the system, the time-step and finite-size extrapolations done, and runqmc’s checks carrying the manual’s advice. The “make the answer better” button, built one default at a time.
• Installation: a download link that names its version, casinocheck extended from one system to several, and known-bad compiler versions refused at compile time.
• A Python layer: an ASE calculator that writes the input, runs runqmc and reads the out file, and a workflow for the time-step and finite-size extrapolations that carries every error bar through to the number quoted.
• In the output, a cost line, the product of run time and energy variance, which says what a smaller error bar will cost before the run is extended; and a line on the health of the branching, the spread of the weights and the rate of catastrophes, beyond what POPSTATS reports.
• Builds and architecture files for today’s machines kept current, ARM (Grace, Apple silicon), AMD, Intel ifx and NVHPC among them, with the suite run on each at every release.

Open: documentation

For users first, then for developers.

• A first week with CASINO: ten worked examples from H₂ to a crystal, each with its reference numbers and the time it takes on a laptop.
• A guide for DFT users: what DMC does for you, what it costs, how to read an error bar, and when a fixed-node answer can be trusted.
• A map of the source for the programming guide: which routine does what, and where a new expectation value, orbital type or interaction goes.
• Every example in the distribution with its reference output beside it, and the gaps in what the examples cover filled; the READMEs and recorded results are already right.

Open: LOUIS

The quantum trajectory code distributed with CASINO.

• The velocity loop: timing a shorter-chain table, and a discrepancy in one dimension. Mike Towler.
• LOUIS and CASINO’s pilot-wave driver reading and writing one trajectory format, so that a film of a two-electron atom and of an ensemble of four thousand atoms come from the same tools.

Open: infelicities

Small, well-defined faults: a first patch for a new developer.

• The FIXMEs in the source (grep finds them).

From the last list: the Scemama probability function is DOMAIN_PROB above; the DMC-MD Green’s function is part of the light-nuclei item; twistav_castep now resumes interrupted runs and its error traps go with it; complex orbitals with a Gaussian basis are not yet done and are now part of the spin-orbit item; symmetry operations in xwfn.data, high-throughput QMC (now the defaults item) and charges, spins and moments (now DOMAIN_PROB) are dropped as separate items; the Blue Gene/Q items and asynchronous I/O were dropped earlier as obsolete.