Long-range Dispersion Research Articles

Riders on the storm? A short note on the biology of Severinia turcomaniae (Saussure, 1872) (Mantodea: Toxoderidae)

Summary New data are presented on the nymphal biology of a praying mantis, Severinia turcomaniae (Saussure, 1872), which has recently spread to the deserts of the Lower Volga Region. Nymphs of 3rd and 4th instars can be encountered in masses in August, suggesting a very short egg development period. They display a clear habitat preference for the tumbleweeds of Ceratocarpus arenarius L. It is hypothesized that the latter may serve for S. turcomaniae nymphs as a unique and efficient means for long-range dispersal and colonization of new localities, which would explain the fast colonization of the territories west of the Volga by S. turcomaniae, but not by Rivetina nana Mistshenko, 1967, another recent invader.

Partial Order–Disorder Transformation of Interpenetrated Porous Coordination Polymers

Controlling interpenetration and flexibility behaviors is intriguing and fundamental for porous coordination polymers. We report exceptional interpenetration behaviors involving controllable partia...

Considerations of Recent All-Atom Lipid Force Field Development.

Molecular simulations of biological molecules require an accurate description of molecular interactions through a force field (FF). The focus of this Perspective is on all-atom lipid FFs. Recent additions to the CHARMM36 lipid FF continue to expand a researcher's ability to probe membrane structure and function with a wide variety of biologically important lipids. Currently, there is an effort to reduce the assumptions in all-atom lipid FFs. The inclusion of long-range dispersion interaction through particle-mesh Ewald is allowing for more accurate descriptions of lipid bilayer and monolayer properties without additional computational cost. Soon, simulations with lipid FFs will no longer depend on short-range cutoffs and will accurately represent long-range dispersion. This requires efficient FF parametrization with an automated approach due to FF complexity. In addition, polarizable FFs for lipids will be important for the next generation of simulations that accurately represent how molecule interactions respond to a varied environment.

Catalytic Conversion of CO and H2 into Hydrocarbons on the Cobalt Co(111) Surface: Implications for the Fischer–Tropsch Process

The Fischer–Tropsch (FT) process consists of the reaction of a synthesis gas (syngas) mixture containing carbon monoxide (CO) and hydrogen (H2), which are polymerized into liquid hydrocarbon chains, often using a cobalt catalyst, although the mechanistic pathway is not yet fully understood. Here, we have employed unrestricted density functional theory calculations with a Hubbard Hamiltonian and long-range dispersion corrections [DFT+U−D3−(BJ)] to investigate the reaction of syngas and the selectivity toward the hydrocarbons formed on the cobalt Co(111) surface. The single CO and dissociated H2 molecules prefer to adsorb at two different types of trigonal surface sites, and we discuss how the interatomic distances, fundamental vibrational modes, charge transfers, surface-free energies, and work functions are modified by the adsorbates. The coadsorption of the syngas molecules in close proximity provides enough energy for the system to cross the saddle points on the minimum energy pathway (MEP), leading to the catalytic hydrogenolysis of the C–O bond. The adsorbed CO, alongside the intermediates CH and OH, are further stabilized when the ratio of equilibrium coverage (C) is CH/CCO,CH,OH > 6:1 under the temperature conditions required for the FT process. We propose several mechanistic pathways to account for the formation of ethane (C2H6), as a model for long-chain hydrocarbons, as well as methane (CH4) which is an undesirable product. The MEPs for these processes show that the coupling of the C–C bond followed by hydrogenation is the most favorable process, which takes precedence over the production of CH4. The termination reaction suggests that water (H2O) remains weakly physisorbed to the surface, allowing the reutilization of its catalytic site. The simulated fundamental vibrational frequencies and scanning tunneling microscopy images of the surface-bound intermediates are in agreement with the available experimental data. Our findings are important in the interpretation of the elementary steps of the FT process on the Co(111) surface.

Mutualistic microalgae co-diversify with reef corals that acquire symbionts during egg development.

The application of molecular genetics has reinvigorated and improved how species are defined and investigated scientifically, especially for morphologically cryptic micro-organisms. Here we show how species recognition improves understanding of the ecology and evolution of mutualisms between reef-building corals and their mutualistic dinoflagellates (i.e. Symbiodiniaceae). A combination of genetic, ecological, and morphological evidence defines two sibling species of Cladocopium (formerly Symbiodinium Clade C), specific only to host corals in the common genus Pocillopora, which transmit their obligate symbionts during oogenesis. Cladocopium latusorum sp. nov. is symbiotic with P. grandis/meandrina while the smaller-celled C. pacificum sp. nov. associates with P. verrucosa. Both symbiont species form mutualisms with Pocillopora that brood their young. Populations of each species, like their hosts, are genetically well connected across the tropical and subtropical Pacific Ocean, indicating a capacity for long-range dispersal. A molecular clock approximates their speciation during the late Pliocene or early Pleistocene as Earth underwent cycles of precipitous cooling and warming; and corresponds to when their hosts were also diversifying. The long temporal and spatial maintenance of high host fidelity, as well as genetic connectivity across thousands of kilometers, indicates that distinct ecological attributes and close evolutionary histories will restrain the adaptive responses of corals and their specialized symbionts to rapid climate warming.

Van der Waals corrected density functionals for cylindrical surfaces: Ammonia and nitrogen dioxide adsorbed on a single-walled carbon nanotube

We extend the damped Zaremba-Kohn model (dZK) for long-range dispersion interaction between a molecule and a planar surface [J. Tao, H. Tang, A. Patra, P. Bhattarai, and J. P. Perdew, Phys. Rev. B 97, 165403 (2018)] to molecules adsorbed on a curved cylindrical surface, and employ this extended model as an additive correction to the semilocal density functionals PBE (Perdew-Burke-Ernzerhof) and SCAN (strongly constrained and appropriately normed). The resulting PBE+vdW (van der Waals)-dZK and SCAN+vdW-dZK are applied to two systems, $\mathrm{N}{\mathrm{H}}_{3}$ and $\mathrm{N}{\mathrm{O}}_{2}$ molecules adsorbed on a single-wall carbon nanotube (CNT), for calculations of binding energies and equilibrium distances. For comparison, the results from vdW nonlocal functionals, such as SCAN+rVV10 and PBE+rVV10, are also presented. The binding energies from PBE+rVV10 (Vydrov and Van Voorhis), SCAN+rVV10, PBE+vdW-dZK, and SCAN+vdW-dZK are about 70--115 meV for the system of $\mathrm{CNT}+\mathrm{N}{\mathrm{H}}_{3}$ and 300--500 meV for the system of $\mathrm{CNT}+\mathrm{N}{\mathrm{O}}_{2}$. The results from PBE+vdW-dZK and SCAN+vdW-dZK are closer to each other than those from PBE+rVV10 and SCAN+rVV10 are. The relatively closer results from PBE+vdW-dZK and SCAN+vdW-dZK indicate the consistency of our developed vdW\ensuremath{-}dZK model for cylindrical surfaces. All methods, including PBE, SCAN, PBE+rVV10, SCAN+rVV10, PBE+vdW-dZK, and SCAN+vdW-dZK, give approximately the same binding energy differences between two adsorption configurations (types I and II) for the two systems. This implies that the two adsorption sites have approximately the same adsorption stability. The exponent of the vdW interaction power law from our vdW-dZK model for the two systems is about 0 at short distance, largely due to the damping factor, and tends slowly to \ensuremath{-}4 to \ensuremath{-}4.5 at distances $D$ about 20--50 \AA{}. At even larger distances, the vdW power-law exponent approaches \ensuremath{-}5. This feature is very similar to the one calculated with random-phase approximation and renormalization group approaches, supporting the applicability of our methods. Our developed vdW-dZK method provides a highly efficient and reliable method for large systems with cylindrical surfaces, such as vdW interactions with nanotubes.

Phylogenomics and Historical Biogeography of Seahorses, Dragonets, Goatfishes, and Allies (Teleostei: Syngnatharia): Assessing Factors Driving Uncertainty in Biogeographic Inferences.

The charismatic trumpetfishes, goatfishes, dragonets, flying gurnards, seahorses, and pipefishes encompass a recently defined yet extraordinarily diverse clade of percomorph fishes-the series Syngnatharia. This group is widely distributed in tropical and warm-temperate regions, with a great proportion of its extant diversity occurring in the Indo-Pacific. Because most syngnatharians feature long-range dispersal capabilities, tracing their biogeographic origins is challenging. Here, we applied an integrative phylogenomic approach to elucidate the evolutionary biogeography of syngnatharians. We built upon a recently published phylogenomic study that examined ultraconserved elements by adding 62 species (total 169 species) and one family (Draconettidae), to cover ca. 25% of the species diversity and all 10 families in the group. We inferred a set of time-calibrated trees and conducted ancestral range estimations. We also examined the sensitivity of these analyses to phylogenetic uncertainty (estimated from multiple genomic subsets), area delimitation, and biogeographic models that include or exclude the jump-dispersal parameter ($j)$. Of the three factors examined, we found that the $j$ parameter has the strongest effect in ancestral range estimates, followed by number of areas defined, and tree topology and divergence times. After accounting for these uncertainties, our results reveal that syngnatharians originated in the ancient Tethys Sea ca. 87 Ma (84-94 Ma; Late Cretaceous) and subsequently occupied the Indo-Pacific. Throughout syngnatharian history, multiple independent lineages colonized the eastern Pacific (6-8 times) and the Atlantic (6-14 times) from their center of origin, with most events taking place following an east-to-west route prior to the closure of the Tethys Seaway ca. 12-18 Ma. Ultimately, our study highlights the importance of accounting for different factors generating uncertainty in macroevolutionary and biogeographic inferences.[Historical biogeography; jump-dispersal parameter; macroevolutionary uncertainty; marine fishes; syngnathiformes; ultraconserved elements].

The primitive brain of early Homo.

The brains of modern humans differ from those of great apes in size, shape, and cortical organization, notably in frontal lobe areas involved in complex cognitive tasks, such as social cognition, tool use, and language. When these differences arose during human evolution is a question of ongoing debate. Here, we show that the brains of early Homo from Africa and Western Asia (Dmanisi) retained a primitive, great ape-like organization of the frontal lobe. By contrast, African Homo younger than 1.5 million years ago, as well as all Southeast Asian Homo erectus, exhibited a more derived, humanlike brain organization. Frontal lobe reorganization, once considered a hallmark of earliest Homo in Africa, thus evolved comparatively late, and long after Homo first dispersed from Africa.

The fragment molecular orbital method combined with density-functional tight-binding and periodic boundary conditions.

The density-functional tight-binding (DFTB) formulation of the fragment molecular orbital method is combined with periodic boundary conditions. Long-range electrostatics and dispersion are evaluated with the Ewald summation technique. The first analytic derivatives of the energy with respect to atomic coordinates and lattice parameters are formulated. The accuracy of the method is established in comparison to numerical gradients and DFTB without fragmentation. The largest elementary cell in this work has 1631 atoms. The method is applied to elucidate the polarization, charge transfer, and interactions in the solution.

Exploring the Limits of the XYG3-Type Doubly Hybrid Approximations for the Main-Group Chemistry: The xDH@B3LYP Model.

Despite being the most accurate class of density functional approximations for the main-group chemistry, doubly hybrid approximations (DHAs) are generally considered to be incomplete in describing the medium- to long-range dispersive interactions. The existing DHAs are often supplemented with empirical long-range dispersion corrections. By using the extensive and chemically diverse GMTKN55 database, we explore the limits of the XYG3-type DHAs using the B3LYP reference orbitals, namely, xDH@B3LYP, with a gradually relaxed constraint on the mixing parameters of DHAs. Our results demonstrate that the xDH@B3LYP model can provide a balanced description of both covalent and noncovalent interactions with the accuracy and robustness comparable to or even better than the very expensive composite methods in wave function theory. Such an accuracy can be achieved without resorting to the use of any long-range correction scheme, shedding new light on the development of DHAs.

A non-empirical intermolecular force-field for trinitrobenzene and its application in crystal structure prediction.

An anisotropic atom-atom distributed intermolecular force-field (DIFF) for rigid trinitrobenzene (TNB) is developed using distributed multipole moments, dipolar polarizabilities, and dispersion coefficients derived from the charge density of the isolated molecule. The short-range parameters of the force-field are fitted to first- and second-order symmetry-adapted perturbation theory dimer interaction energy calculations using the distributed density-overlap model to guide the parameterization of the short-range anisotropy. The second-order calculations are used for fitting the damping coefficients of the long-range dispersion and polarization and also for relaxing the isotropic short-range coefficients in the final model, DIFF-srL2(rel). We assess the accuracy of the unrelaxed model, DIFF-srL2(norel), and its equivalent without short-range anisotropy, DIFF-srL0(norel), as these models are easier to derive. The model potentials are contrasted with empirical models for the repulsion-dispersion fitted to organic crystal structures with multipoles of iterated stockholder atoms (ISAs), FIT(ISA,L4), and with Gaussian Distributed Analysis (GDMA) multipoles, FIT(GDMA,L4), commonly used in modeling organic crystals. The potentials are tested for their ability to model the solid state of TNB. The non-empirical models provide more reasonable relative lattice energies of the three polymorphs of TNB and propose more sensible hypothetical structures than the empirical force-field (FIT). The DIFF-srL2(rel) model successfully has the most stable structure as one of the many structures that match the coordination sphere of form III. The neglect of the conformational flexibility of the nitro-groups is a significant approximation. This methodology provides a step toward force-fields capable of representing all phases of a molecule in molecular dynamics simulations.

CHARMM36 Lipid Force Field with Explicit Treatment of Long-Range Dispersion: Parametrization and Validation for Phosphatidylethanolamine, Phosphatidylglycerol, and Ether Lipids.

Long-range Lennard-Jones (LJ) interactions have been incorporated into the CHARMM36 (C36) lipid force field (FF) using the LJ particle-mesh Ewald (LJ-PME) method in order to remove the inconsistency of bilayer and monolayer properties arising from the exclusion of long-range dispersion [Yu, Y.; Semi-automated Optimization of the CHARMM36 Lipid Force Field to Include Explicit Treatment of Long-Range Dispersion. J. Chem. Theory Comput. 2021, 10.1021/acs.jctc.0c01326. (preceding article in this issue)]. The new FF is denoted C36/LJ-PME. While the first optimization was based on three phosphatidylcholines (PCs), this work extends the validation and parametrization to more lipids including PC, phosphatidylethanolamine (PE), phosphatidylglycerol (PG), and ether lipids. The agreement with experimental structure data is excellent for PC, PE, and ether lipids. C36/LJ-PME also compares favorably with scattering data of PG bilayers but less so with NMR deuterium order parameters of 1,2-dimyristoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (DMPG) at 303.15 K, indicating a need for future optimization regarding PG-specific parameters. Frequency dependence of NMR T1 spin-lattice relaxation times is well-described by C36/LJ-PME, and the overall agreement with experiment is comparable to C36. Lipid diffusion is slower than C36 due to the added long-range dispersion causing a higher viscosity, although it is still too fast compared to experiment after correction for periodic boundary conditions. When using a 10 Å real-space cutoff, the simulation speed of C36/LJ-PME is roughly equal to C36. While more lipids will be incorporated into the FF in the future, C36/LJ-PME can be readily used for common lipids and extends the capability of the CHARMM FF by supporting monolayers and eliminating the cutoff dependence.

Semi-automated Optimization of the CHARMM36 Lipid Force Field to Include Explicit Treatment of Long-Range Dispersion.

The development of the CHARMM lipid force field (FF) can be traced back to the early 1990s with its current version denoted CHARMM36 (C36). The parametrization of C36 utilized high-level quantum mechanical data and free energy calculations of model compounds before parameters were manually adjusted to yield agreement with experimental properties of lipid bilayers. While such manual fine-tuning of FF parameters is based on intuition and trial-and-error, automated methods can identify beneficial modifications of the parameters via their sensitivities and thereby guide the optimization process. This work introduces a semi-automated approach to reparametrize the CHARMM lipid FF with consistent inclusion of long-range dispersion through the Lennard-Jones particle-mesh Ewald (LJ-PME) approach. The optimization method is based on thermodynamic reweighting with regularization with respect to the C36 set. Two independent optimizations with different topology restrictions are presented. Targets of the optimizations are primarily liquid crystalline phase properties of lipid bilayers and the compression isotherm of monolayers. Pair correlation functions between water and lipid functional groups in aqueous solution are also included to address headgroup hydration. While the physics of the reweighting strategy itself is well-understood, applying it to heterogeneous, complex anisotropic systems poses additional challenges. These were overcome through careful selection of target properties and reweighting settings allowing for the successful incorporation of the explicit treatment of long-range dispersion, and we denote the newly optimized lipid force field as C36/LJ-PME. The current implementation of the optimization protocol will facilitate the future development of the CHARMM and related lipid force fields.

Genome sequences reveal global dispersal routes and suggest convergent genetic adaptations in seahorse evolution

Seahorses have a circum-global distribution in tropical to temperate coastal waters. Yet, seahorses show many adaptations for a sedentary, cryptic lifestyle: they require specific habitats, such as seagrass, kelp or coral reefs, lack pelvic and caudal fins, and give birth to directly developed offspring without pronounced pelagic larval stage, rendering long-range dispersal by conventional means inefficient. Here we investigate seahorses’ worldwide dispersal and biogeographic patterns based on a de novo genome assembly of Hippocampus erectus as well as 358 re-sequenced genomes from 21 species. Seahorses evolved in the late Oligocene and subsequent circum-global colonization routes are identified and linked to changing dynamics in ocean currents and paleo-temporal seaway openings. Furthermore, the genetic basis of the recurring “bony spines” adaptive phenotype is linked to independent substitutions in a key developmental gene. Analyses thus suggest that rafting via ocean currents compensates for poor dispersal and rapid adaptation facilitates colonizing new habitats.

FALL3D-8.0: a computational model for atmospheric transport and deposition of particles, aerosols and radionuclides – Part 2: Model validation

Abstract. This paper presents model validation results for the latest version release of the FALL3D atmospheric transport model. The code has been redesigned from scratch to incorporate different categories of species and to overcome legacy issues that precluded its preparation towards extreme-scale computing. The model validation is based on the new FALL3D-8.0 test suite, which comprises a set of four real case studies that encapsulate the major features of the model; namely, the simulation of long-range fine volcanic ash dispersal, volcanic SO2 dispersal, tephra fallout deposits and the dispersal and deposition of radionuclides. The first two test suite cases (i.e. the June 2011 Puyehue-Cordón Caulle ash cloud and the June 2019 Raikoke SO2 cloud) are validated against geostationary satellite retrievals and demonstrate the new FALL3D data insertion scheme. The metrics used to validate the volcanic ash and SO2 simulations are the structure, amplitude and location (SAL) metric and the figure of merit in space (FMS). The other two test suite cases (i.e. the February 2013 Mt. Etna ash cloud and associated tephra fallout deposit, and the dispersal of radionuclides resulting from the 1986 Chernobyl nuclear accident) are validated with scattered ground-based observations of deposit load and local particle grain size distributions and with measurements from the Radioactivity Environmental Monitoring database. For validation of tephra deposit loads and radionuclides, we use two variants of the normalised root-mean-square error metric. We find that FALL3D-8.0 simulations initialised with data insertion consistently improve agreement with satellite retrievals at all lead times up to 48 h for both volcanic ash and SO2 simulations. In general, SAL scores lower than 1.5 and FMS scores greater than 0.40 indicate acceptable agreement with satellite retrievals of volcanic ash and SO2. In addition, we show very good agreement, across several orders of magnitude, between the model and observations for the 2013 Mt. Etna and 1986 Chernobyl case studies. Our results, along with the validation datasets provided in the publicly available test suite, form the basis for future improvements to FALL3D (version 8 or later) and also allow for model intercomparison studies.

Ember risk modelling for improved wildfire risk management in the peri-urban fringes

Although embers are a leading cause of house loss and damage from wildfire, risk reduction activities such as vegetation management and planning requirements do not adequately account for the risk of embers. To help address this, a model of the potential ember risk is presented. It takes into account local vegetation and background wind conditions for both long-range and short-range ember dispersal. This model is developed to provide indications of the importance of embers to the household-level wildfire risk in communities at the wildland-urban interface. The model provides information for householders to improve understanding of the nature of ember risk within the community, and to assist planning in order to respond to that risk. A case study of the 2015 Warringine Park (Coastal Section) bushfire in Australia is presented, which demonstrates how the model could be used to assist with community planning. The utility of this outcome for community and household level wildfire planning and preparation is discussed.

Spatio-Temporal Representativeness of Air Quality Monitoring Stations in Mexico City: Implications for Public Health.

Assessment of the air quality in metropolitan areas is a major challenge in environmental sciences. Issues related include the distribution of monitoring stations, their spatial range, or missing information. In Mexico City, stations have been located spanning the entire Metropolitan zone for pollutants, such as CO, NO2, O3, SO2, PM2.5, PM10, NO, NOx, and PMCO. A fundamental question is whether the number and location of such stations are adequate to optimally cover the city. By analyzing spatio-temporal correlations for pollutant measurements, we evaluated the distribution and performance of monitoring stations in Mexico City from 2009 to 2018. Based on our analysis, air quality evaluation of those contaminants is adequate to cover the 16 boroughs of Mexico City, with the exception of SO2, since its spatial range is shorter than the one needed to cover the whole surface of the city. We observed that NO and NOx concentrations must be taken into account since their long-range dispersion may have relevant consequences for public health. With this approach, we may be able to propose policy based on systematic criteria to locate new monitoring stations.

Local embedding of coupled cluster theory into the random phase approximation using plane waves.

We present an embedding approach to treat local electron correlation effects in periodic environments. In a single consistent framework, our plane wave based scheme embeds a local high-level correlation calculation [here, Coupled Cluster (CC) theory], employing localized orbitals, into a low-level correlation calculation [here, the direct Random Phase Approximation (RPA)]. This choice allows for an accurate and efficient treatment of long-range dispersion effects. Accelerated convergence with respect to the local fragment size can be observed if the low-level and high-level long-range dispersions are quantitatively similar, as is the case for CC in RPA. To demonstrate the capabilities of the introduced embedding approach, we calculate adsorption energies of molecules on a surface and in a chabazite crystal cage, as well as the formation energy of a lattice impurity in a solid at the level of highly accurate many-electron perturbation theories. The absorption energy of a methane molecule in a zeolite chabazite is converged with an error well below 20 meV at the CC level. As our largest periodic benchmark system, we apply our scheme to the adsorption of a water molecule on titania in a supercell containing more than 1000 electrons.

CO2 reduction to acetic acid on the greigite Fe3S4{111} surface.

Acetic acid (CH3-COOH) is an important commodity chemical widely used in a myriad of industrial processes, whose production still largely depends on homogeneous catalysts based on expensive rare metals. Here, we report a computational study on the formation of CH3-COOH from carbon dioxide (CO2) as an alternative chemical feedstock on the {111} surface of the low-cost greigite Fe3S4 catalyst. We have used density functional theory calculations with a Hubbard Hamiltonian approach and long-range dispersion corrections (DFT+U-D2) to simulate the various stages of the direct combination of C1 species of different composition to produce glyoxylic acid (CHO-COOH) as a key intermediate in the formation of CH3-COOH. Three reaction mechanisms are considered: (i) the main pathway where the direct formation of the C-C bond takes place spontaneously, followed by a step-wise reduction of CHO-CHOO to CH3-COOH; and the competitive pathways for the non-promoted and H-promoted elimination of hydroxy groups (OH) and water (H2O), respectively from (ii) the carboxyl; and (iii) the carbonyl end of the glyoxylate intermediates. The thermodynamic and kinetic profiles show that the energies for the intermediates on the main pathway are very similar for the two catalytic sites considered, although the activation energies are somewhat larger for the exposed tetrahedral iron (FeA) ion. In most cases, the intermediates for the deoxygenation of the carboxylic acid are less stable than the intermediates on the main pathway, which suggests that the molecule prefers to lose the carbonylic oxygen. The suitable surface properties of the Fe3S4{111} surface show that this material could be a promising sustainable catalyst in future technologies for the conversion of CO2 into organic acid molecules of commercial interest.

Assessing cathode property prediction via exchange-correlation functionals with and without long-range dispersion corrections.

We benchmark calculated interlayer spacings, average topotactic voltages, thermodynamic stabilities, and band gaps in layered lithium transition-metal oxides (TMOs) and their de-lithiated counterparts, which are used in lithium-ion batteries as positive electrode materials, against available experimental data. Specifically, we examine the accuracy of properties calculated within density functional theory (DFT) using eight different treatments of electron exchange-correlation: the strongly constrained and appropriately normed (SCAN) and Perdew-Burke-Ernzerhof (PBE) density functionals, Hubbard-U-corrected SCAN and PBE (i.e., SCAN+U and PBE+U), and SCAN(+U) and PBE(+U) with added long-range dispersion (D) interactions (i.e., DFT(+U)+D). van der Waals interactions are included respectively via the revised Vydrov-Van Voorhis (rVV10) for SCAN(+U) and the DFT-D3 for PBE(+U). We find that SCAN-based functionals predict larger voltages due to an underestimation of stability of the MO2 systems, while also predicting smaller interlayer spacings compared to their PBE-based counterparts. Furthermore, adding dispersion corrections to PBE has a greater effect on voltage predictions and interlayer spacings than with SCAN, indicating that DFT-SCAN - despite being a ground-state theory - fortuitously captures some short and medium-range dispersion interactions better than PBE. While SCAN-based and PBE-based functionals yield qualitatively similar band gap predictions, there is no significant quantitative improvement of SCAN-based functionals over the corresponding PBE-based versions. Finally, we expect SCAN-based functionals to yield more accurate property predictions than the respective PBE-based functionals for most TMOs, given SCAN's stronger theoretical underpinning and better predictions of systematic trends in interlayer spacings, intercalation voltages, and band gaps obtained in this work.