CO2-philicity to CO2-phobicity Transition on Smooth and Stochastic Rough Cu-like Substrate Surfaces.

CO2 on metal substrates is essential to CO2 liquefaction and transportation of CO2, yet the manipulation of the wettability of the CO2 and the elucidation of its underlying mechanism have not been fully achieved. Here, using molecular dynamics simulations, we report CO2 wetting characteristics on both smooth and stochastic rough Cu-like substrate surfaces. The results indicate that the apparent contact angle (CA) of the CO2 droplet on the smooth surface decreases from 180° to 0° as the CO2-solid characteristic interaction energy increases from 0.002 to 0.016 eV. In addition, the CAs become greater with increasing the density of surface asperities, regardless of the intrinsic surface wettability. This is attributed to the capillary drying-out of liquid CO2 molecules in gaps between surface asperities at the three-phase contact line of the droplet, which is usually overlooked in previous theoretical studies. Notably, the intrinsically CO2-philic surface transforms to the CO2-phobic due to an increase in the density of surface rugosity. Moreover, we verify the range of applicability of the CA prediction models concerning the nanoscale asperities. This work is beneficial for fully understanding the influence of nanoscale surface topography on CO2 wettability and shedding light on the design of functionalized and patterned surfaces to manipulate CO2 wettability.

Kinetic Mechanism of Surfactant-Based Molecular Recognition: Selective Permeability across an Oil-Water Interface Regulated by Supramolecular Aggregates.

Lipids are known to play a vital role in the molecular organization of all cellular life. Molecular recognition is another fundamental biological process that is generally attributed to biological polymers, such as proteins and nucleic acids. However, there is evidence that aggregates of lipids and lipid-like molecules are also capable of selectively binding to or regulating the partitioning of other molecules. We previously demonstrated that a model two-phase octanol/water system can selectively partition Red 40 and Blue 1 dyes added to an aqueous phase, with the selectivity depending on the surfactant (e.g., cetyltrimethylammonium bromide) dissolved in the organic phase. Here, we elucidate the mechanism of molecular recognition in this system by using quantitative partitioning experiments and molecular dynamics (MD) simulations. Our results indicate that the selectivity for the red dye is thermodynamically favored at all surfactant concentrations, while selectivity for the blue dye is kinetically favored at high surfactant concentrations. The kinetic selectivity for the blue dye correlates with the presence of molecular aggregation at the oil-water interface. Coarse-grained MD simulations elucidate nanoscale supramolecular structures that can preferentially bind one small molecule rather than another at an interface, providing a selectively permeable barrier in the absence of proteins. The results suggest a new supramolecular mechanism for molecular recognition with potential applications in drug delivery, drug discovery, and biosensing.

Dissolution of Forsterite Surface in Brine at CO2 Geo-storage Conditions: Insights from Molecular Dynamic Simulations.

Evaluating the long-term security of geological deep saline aquifers to store CO2 requires a comprehensive understanding of mineral dissolution properties. Molecular dynamics simulations are performed to study the dissolution of forsterite in deep saline aquifers. The forsterite surface is found to be covered by three H2O molecular layers, hindering CO2 from directly contacting the surface. The dissolution rates at 350 K are increased by more than 1012 with the presence of Mg defects or salt ions in solutions. The more disordered surface in pure water caused by Mg defects accounts for the acceleration of dissolution, while absorbed Cl- ions on the surface in NaCl and KCl solutions accelerate the dissolution through electrostatic interactions. Comparatively, the frequent attacks from alkaline earth cations in MgCl2 and CaCl2 solutions to the surface contribute to the enhanced dissolution. In the acidic H3OCl solution, the electrostatic interactions between O atoms in H3O+ and the surface facilitate the dissolution. Interestingly, the ionic clusters of CO32-/HCO3- and Na+ in Na2CO3/NaHCO3 solution promote the dissolution process. This work provides molecular insights into forsterite dissolution in deep saline aquifers and guidance toward the optimization of CO2 geo-storage conditions.

Thermal batteries based on inverse barocaloric effects.

To harvest and reuse low-temperature waste heat, we propose and realize an emergent concept-barocaloric thermal batteries based on the large inverse barocaloric effect of ammonium thiocyanate (NH4SCN). Thermal charging is initialized upon pressurization through an order-to-disorder phase transition, and the discharging of 43 J g-1 takes place at depressurization, which is 11 times more than the input mechanical energy. The thermodynamic equilibrium nature of the pressure-restrained heat-carrying phase guarantees stable long-duration storage. The barocaloric thermal batteries reinforced by their solid microscopic mechanism are expected to substantially advance the ability to take advantage of waste heat.

Modeling the viscoelastic relaxation dynamics of soft particles via molecular dynamics simulation-informed multi-dimensional transition-state theory.

Viscoelastic soft colloidal particles have been widely explored in mechanical, chemical, pharmaceutical and other engineering applications due to their unique combination of viscosity and elasticity. The characteristic viscoelastic relaxation time shows an Arrhenius-type (or super-Arrhenius due to temperature-dependent transition attempts) thermally-activated behavior, but a holistic explanation from the relevant transition-state theory remains elusive. In this paper, the viscoelastic relaxation times of Lennard-Jones soft colloidal particle systems, including a single particle type system and a binary particle mixture based on the Kob-Andersen model, are determined using molecular dynamics (MD) simulations as the benchmark. First, the particle systems show a non-Maxwellian behavior after comparing the MD-predicted viscoelastic relaxation time and dynamic moduli (storage and loss modulus) to the classic Maxwell viscoelastic model and the recent particle local connectivity theory. Surprisingly, neither the Maxwell relaxation time τMaxwell (obtained from the static shear viscosity η and the high-frequency shear modulus G∞) nor the particle local connectivity lifetime τLC can capture the super-Arrhenius temperature-dependent behavior in the MD-predicted relaxation time τMD. Then, the particle dissociation and association transition kinetics, fractal dimensions of the particle systems, and neighbor particle structure (obtained from the radial distribution functions) are shown to collectively determine the viscoelastic relaxation time. These factors are embedded into a new multi-dimensional transition kinetics model to directly estimate the viscoelastic relaxation time τModel, which is found to agree with the MD-predicted τMD remarkably well. This work highlights the microscopic origin of viscoelastic relaxation dynamics of soft colloidal particles, and theoretically connects rheological dynamics and transition kinetics in soft matters.

Hybridly double-crosslinked carbon nanotube networks with combined strength and toughness via cooperative energy dissipation.

Although chemical crosslinking has been extensively explored to enhance the mechanical properties of network-type materials for structural and energy (electrochemical, thermal, etc.) applications, loading-induced energy dissipations usually occur through a single channel that either leads to network brittleness or low strength/stiffness. In this work, we apply coarse-grained molecular dynamics simulations to explore the potential of hybridly double-crosslinked carbon nanotube (CNT) networks as a light weight functional material with combined strength and toughness. While increasing the crosslinking density or strong crosslink composition may, in general, enhance the strength and toughness, further increasing the two parameters would surprisingly lead to deteriorated strength and toughness. We find that double-crosslinked networks can nicely achieve cooperative energy dissipation with minimal structural damage. In particular, the weak crosslinks serve as "sacrificial bonds" to dissipate elastic energies from external loading, while the strong crosslinks act as "structure holders" and break at a much later stage during the tensile test. Therefore, the combination of more than one type of crosslinking with hybrid potential energy landscapes and breaking time scales can prevent premature simultaneous breaking of multiple strong crosslinks. By deploying intermediate amounts of weak and strong crosslinks, we observe an outstanding density-normalized strength of 227-2130 kPa m3 kg-1 as compared to many structural materials and advanced nanocomposites. The crosslinking strategies developed here would pave new avenues for the rational design of functional network materials beyond CNTs, such as hydrogels, nanofibers, and nanocomposites.

Symmetry Breaking Induced Anisotropic Carrier Transport and Remarkable Thermoelectric Performance in Mixed Halide Perovskites CsPb(I1-xBrx)3.

We present a combination of first-principles calculations and the Boltzmann transport theory to understand the carrier transport and thermoelectric performance of mixed halide perovskite alloys CsPb(I1-xBrx)3 with different Br compositions. Our computational results correlate the conduction band splitting in CsPb(I1-xBrx)3 to the significant anisotropy in their carrier transport properties, such as effective masses and deformation potential constants. Such band splitting originates from the symmetry-broken crystal structures of CsPb(I1-xBrx)3 polymorphs: with residue stresses/strains in asymmetric CsPb(I1-xBrx)3, nondegenerate orbitals reconstruct the conduction band and reduce the Pb-halide antibonding character along certain directions. While the Seebeck coefficient (S) and the relaxation time-normalized electrical conductivity (σ/τ) show weak directional anisotropy, the carrier relaxation time (τ) is highly direction-dependent. The reconstruction of the conduction band finally leads to significantly anisotropic and enhanced thermoelectric power factors (PF = S2σ) in CsPb(I1-xBrx)3 compared to those in pure CsPbI3 and CsPbBr3, showing anomalous nonlinear alloy behavior. A delicate balance between S2σ and combined measurement of the carrier effective mass and deformation potential constant, m*EDP, is confirmed. The lattice thermal conductivities of CsPb(I1-xBrx)3 are significantly suppressed compared to those of their pure counterparts due to strong mass disordering and strain fields upon halogen substitution. As a result, symmetry breaking in CsPb(I1-xBrx)3 leads to anisotropy in carrier transport, high PF, and scattered phonon transport (ultralow thermal conductivity), concurrently contributing to their promising thermoelectric figures of merit (ZT) up to 1.7 at room temperature. The principles behind the asymmetry-induced factors would serve as new design concepts to tailor the thermoelectric properties of alloys, mixtures, superlattices, and low-dimensional materials.

Competing Dissolution Pathways and Ligand Passivation-Enhanced Interfacial Stability of Hybrid Perovskites with Liquid Water.

Material instability issues, especially moisture degradation in ambient operating environments, limit the practical application of hybrid perovskite in photovoltaic and light-emitting devices. Very recent experiments demonstrate that ligand passivation can effectively improve the surface moisture tolerance of hybrid perovskites. In this work, the interfacial stability of as-synthesized pristine and alkylammonium-passivated methylammonium lead iodide (MAPbI3) with liquid water is systematically investigated using molecular dynamics simulations and reaction kinetics models. Interestingly, the more hydrophilic [PbI2]0 surface is more stable than the less hydrophilic [MAI]0 surface because of the higher polarity of the former surface. Linear alkylammoniums significantly stabilize the [MAI]0 surface with highly reduced (by 1-2 orders of magnitude) dissociation rates of both MA+ and ligands themselves, while branched ligands, surprisingly, lead to higher dissociation rates as the surface coverage increases. Such anomalous behavior is attributed to the aggregation-assisted dissolution of surfactant-like ligands as micelles during the degradation process. Short-chain linear alkylammonium at the full surface coverage is found to be the optimal ligand to stabilize the [MAI]0 surface. This work not only provides fundamental insights into the ionic dissolution pathways and mechanisms of hybrid perovskites in water but also inspires the design of highly stable hybrid perovskites with ligand passivation layers. The computational framework developed here is also transferrable to the investigation of surface passivation chemistry for weak ionic materials in general.

Colossal barocaloric effects in plastic crystals.

Refrigeration is of vital importance for modern society-for example, for food storage and air conditioning-and 25 to 30 per cent of the world's electricity is consumed for refrigeration1. Current refrigeration technology mostly involves the conventional vapour compression cycle, but the materials used in this technology are of growing environmental concern because of their large global warming potential2. As a promising alternative, refrigeration technologies based on solid-state caloric effects have been attracting attention in recent decades3-5. However, their application is restricted by the limited performance of current caloric materials, owing to small isothermal entropy changes and large driving magnetic fields. Here we report colossal barocaloric effects (CBCEs) (barocaloric effects are cooling effects of pressure-induced phase transitions) in a class of disordered solids called plastic crystals. The obtained entropy changes in a representative plastic crystal, neopentylglycol, are about 389 joules per kilogram per kelvin near room temperature. Pressure-dependent neutron scattering measurements reveal that CBCEs in plastic crystals can be attributed to the combination of extensive molecular orientational disorder, giant compressibility and highly anharmonic lattice dynamics of these materials. Our study establishes the microscopic mechanism of CBCEs in plastic crystals and paves the way to next-generation solid-state refrigeration technologies.

Suppressed phase separation of mixed-halide perovskites confined in endotaxial matrices.

The functionality and performance of a semiconductor is determined by its bandgap. Alloying, as for instance in InxGa1-xN, has been a mainstream strategy for tuning the bandgap. Keeping the semiconductor alloys in the miscibility gap (being homogeneous), however, is non-trivial. This challenge is now being extended to halide perovskites - an emerging class of photovoltaic materials. While the bandgap can be conveniently tuned by mixing different halogen ions, as in CsPb(BrxI1-x)3, the so-called mixed-halide perovskites suffer from severe phase separation under illumination. Here, we discover that such phase separation can be highly suppressed by embedding nanocrystals of mixed-halide perovskites in an endotaxial matrix. The tuned bandgap remains remarkably stable under extremely intensive illumination. The agreement between the experiments and a nucleation model suggests that the size of the nanocrystals and the host-guest interfaces are critical for the photo-stability. The stabilized bandgap will be essential for the development of perovskite-based optoelectronics, such as tandem solar cells and full-color LEDs.

Aggregation-induced emission in lamellar solids of colloidal perovskite quantum wells.

The outstanding excitonic properties, including photoluminescence quantum yield (ηPL), of individual, quantum-confined semiconductor nanoparticles are often significantly quenched upon aggregation, representing the main obstacle toward scalable photonic devices. We report aggregation-induced emission phenomena in lamellar solids containing layer-controlled colloidal quantum wells (QWs) of hybrid organic-inorganic lead bromide perovskites, resulting in anomalously high solid-state ηPL of up to 94%. Upon forming the QW solids, we observe an inverse correlation between exciton lifetime and ηPL, distinct from that in typical quantum dot solid systems. Our multiscale theoretical analysis reveals that, in a lamellar solid, the collective motion of the surface organic cations is more restricted to orient along the [100] direction, thereby inducing a more direct bandgap that facilitates radiative recombination. Using the QW solids, we demonstrate ultrapure green emission by completely downconverting a blue gallium nitride light-emitting diode at room temperature, with a luminous efficacy higher than 90 lumen W-1 at 5000 cd m-2, which has never been reached in any nanomaterial assemblies by far.

Doping-Driven Wettability of Two-Dimensional Materials: A Multiscale Theory.

Engineering molecular interactions at two-dimensional (2D) materials interfaces enables new technological opportunities in functional surfaces and molecular epitaxy. Understanding the wettability of 2D materials represents the crucial first step toward quantifying the interplay between the interfacial forces and electric potential of 2D materials interfaces. Here we develop the first theoretical framework to model the wettability of the doped 2D materials by properly bridging the multiscale physical phenomena at the 2D interfaces, including (i) the change of 2D materials surface energy (atomistic scale, several angstroms), (ii) the molecular reorientation of liquid molecules adjacent to the interface (molecular scale, 100-101 nm), and (iii) the electrical double layer (EDL) formed in the liquid phase (mesoscopic scales, 100-104 nm). The latter two effects are found to be the major mechanisms responsible for the contact angle change upon doping, while the surface energy change of a pure 2D material has no net effect on the wetting property. When the doping level is electrostatically tuned, we demonstrate that 2D materials with high quantum capacitances (e.g., transition metal dichalcogenides, TMDCs) possess a wider range of tunability in the interfacial tension, under the same applied gate voltage. Furthermore, practical considerations such as defects and airborne contamination are also quantitatively discussed. Our analysis implies that the doping level can be another variable to modulate the wettability at 2D materials interfaces, as well as the molecular packing behavior on a 2D material-coated surface, essentially facilitating the interfacial engineering of 2D materials.

Polar rotor scattering as atomic-level origin of low mobility and thermal conductivity of perovskite CH3NH3PbI3.

Perovskite CH3NH3PbI3 exhibits outstanding photovoltaic performances, but the understanding of the atomic motions remains inadequate even though they take a fundamental role in transport properties. Here, we present a complete atomic dynamic picture consisting of molecular jumping rotational modes and phonons, which is established by carrying out high-resolution time-of-flight quasi-elastic and inelastic neutron scattering measurements in a wide energy window ranging from 0.0036 to 54 meV on a large single crystal sample, respectively. The ultrafast orientational disorder of molecular dipoles, activated at ∼165 K, acts as an additional scattering source for optical phonons as well as for charge carriers. It is revealed that acoustic phonons dominate the thermal transport, rather than optical phonons due to sub-picosecond lifetimes. These microscopic insights provide a solid standing point, on which perovskite solar cells can be understood more accurately and their performances are perhaps further optimized.

Molecular Gibbs Surface Excess and CO2-Hydrate Density Determine the Strong Temperature- and Pressure-Dependent Supercritical CO2-Brine Interfacial Tension.

In CO2 geological storage, the interfacial tension (IFT) between supercritical CO2 and brine is critical for the storage capacitance design to prevent CO2 leakage. IFT relies not only on the interfacial molecule properties but also on the environmental conditions at different storage sites. In this paper, supercritical CO2-NaCl solution systems are modeled at 343-373 K and 6-35 MPa under the salinity of 1.89 mol/L using molecular dynamics simulations. After computing and comparing the molecular density profile across the interface, the atomic radial distribution function, the molecular orientation distribution, the molecular Gibbs surface excess (derived from the molecular density profile), and the CO2-hydrate number density under the above environmental conditions, we confirm that only the molecular Gibbs surface excess of CO2 molecules and the CO2-hydrate number density correlate strongly with the temperature- and pressure-dependent IFTs. We also compute the populations of two distinct CO2-hydrate structures (T-type and H-type) and attribute the observed dependence of IFTs to the dominance of the more stable, surfactant-like T-type CO2-hydrates at the interface. On the basis of these new molecular mechanisms behind IFT variations, this study could guide the rational design of suitable injecting environmental pressure and temperature conditions. We believe that the above two molecular-level metrics (Gibbs surface excess and hydrate number density) are of great fundamental importance for understanding the supercritical CO2-water interface and engineering applications in geological CO2 storage.

Low-Dimensional Organic Tin Bromide Perovskites and Their Photoinduced Structural Transformation.

Hybrid organic-inorganic metal halide perovskites possess exceptional structural tunability, with three- (3D), two- (2D), one- (1D), and zero-dimensional (0D) structures on the molecular level all possible. While remarkable progress has been realized in perovskite research in recent years, the focus has been mainly on 3D and 2D structures, with 1D and 0D structures significantly underexplored. The synthesis and characterization of a series of low-dimensional organic tin bromide perovskites with 1D and 0D structures is reported. Using the same organic and inorganic components, but at different ratios and reaction conditions, both 1D (C4 N2 H14 )SnBr4 and 0D (C4 N2 H14 Br)4 SnBr6 can be prepared in high yields. Moreover, photoinduced structural transformation from 1D to 0D was investigated experimentally and theoretically in which photodissociation of 1D metal halide chains followed by structural reorganization leads to the formation of a more thermodynamically stable 0D structure.

Electrochemical Doping of Halide Perovskites with Ion Intercalation.

Halide perovskites have recently been investigated for various solution-processed optoelectronic devices. The majority of studies have focused on using intrinsic halide perovskites, and the intentional incoporation of dopants has not been well explored. In this work, we discovered that small alkali ions, including lithium and sodium ions, could be electrochemically intercalated into a variety of halide and pseudohalide perovskites. The ion intercalation caused a lattice expansion of the perovskite crystals and resulted in an n-type doping of the perovskites. Such electrochemical doping improved the conductivity and changed the color of the perovskites, leading to an electrochromism with more than 40% reduction of transmittance in the 450-850 nm wavelength range. The doped perovskites exhibited improved electron injection efficiency into the pristine perovskite crystals, resulting in bright light-emitting diodes with a low turn-on voltage.

Understanding the colloidal dispersion stability of 1D and 2D materials: Perspectives from molecular simulations and theoretical modeling.

The colloidal dispersion stability of 1D and 2D materials in the liquid phase is critical for scalable nano-manufacturing, chemical modification, composites production, and deployment as conductive inks or nanofluids. Here, we review recent computational and theoretical studies carried out by our group to model the dispersion stability of 1D and 2D materials, including single-walled carbon nanotubes, graphene, and graphene oxide in aqueous surfactant solutions or organic solvents. All-atomistic (AA) molecular dynamics (MD) simulations can probe the molecular level details of the adsorption morphology of surfactants and solvents around these materials, as well as quantify the interaction energy between the nanomaterials mediated by surfactants or solvents. Utilizing concepts from reaction kinetics and diffusion, one can directly predict the rate constants for the aggregation kinetics and dispersion life times using MD outputs. Furthermore, the use of coarse-grained (CG) MD simulations allows quantitative prediction of surfactant adsorption isotherms. Combined with the Poisson-Boltzmann equation, the Langmuir isotherm, and the DLVO theory, one can directly use CGMD outputs to: (i) predict electrostatic potentials around the nanomaterial, (ii) correlate surfactant surface coverages with surfactant concentrations in the bulk dispersion medium, and (iii) determine energy barriers against coagulation. Finally, we discuss challenges associated with studying emerging 2D materials, such as, hexagonal boron nitride (h-BN), phosphorene, and transition metal dichalcogenides (TMDCs), including molybdenum disulfide (MoS2). An outlook is provided to address these challenges with plans to develop force-field parameters for MD simulations to enable predictive modeling of emerging 2D materials in the liquid phase.

Probing the Soft and Nanoductile Mechanical Nature of Single and Polycrystalline Organic-Inorganic Hybrid Perovskites for Flexible Functional Devices.

Although organic-inorganic hybrid perovskites have been extensively investigated for promising applications in energy-related devices, their mechanical properties, which restrict their practical deployment as flexible and wearable devices, have been largely unexplored at the atomistic level. Toward this level of understanding, we predict the elastic constant matrix and various elastic properties of CH3NH3PbI3 (MAPbI3) using atomistic simulations. We find that single-crystalline MAPbI3 is much stiffer and exhibits higher ultimate tensile strength than polycrystalline samples, but the latter exhibit unexpected, greatly enhanced nanoductility and fracture toughness, resulting from the extensive amorphization during the yielding process. More interestingly, polycrystalline MAPbI3 exhibits inverse Hall-Petch grain-boundary strengthening effect, in which the yield stress is reduced when decreasing the grain size, due to amorphous grain boundaries. By monitoring the centro-symmetry parameter and local stress evolution, we confirm the soft and nanoductile nature of defective MAPbI3 with a crack. By conducting atomic stress decomposition, we attribute such fracture toughness primarily to the strong electrostatic interactions between the ionic components. The observed limited brittle fracture behavior is attributed to the transformation of partial edge dislocations to disordered atoms (nanovoid formation). A significant plastic deformation region is observed when nanovoids enlarge and coalesce with adjacent ones, which ultimately leads to crack propagations via ionic-chain breaking. After comparing with traditional inorganic energy-related materials, we find that hybrid perovskites are more compressible and can absorb more strain energy before fracture, which makes them well suited for wearable functional devices with high mechanical flexibility and robustness.

Efficient Blue Electroluminescence Using Quantum-Confined Two-Dimensional Perovskites.

Solution-processed hybrid organic-inorganic lead halide perovskites are emerging as one of the most promising candidates for low-cost light-emitting diodes (LEDs). However, due to a small exciton binding energy, it is not yet possible to achieve an efficient electroluminescence within the blue wavelength region at room temperature, as is necessary for full-spectrum light sources. Here, we demonstrate efficient blue LEDs based on the colloidal, quantum-confined 2D perovskites, with precisely controlled stacking down to one-unit-cell thickness (n = 1). A variety of low-k organic host compounds are used to disperse the 2D perovskites, effectively creating a matrix of the dielectric quantum wells, which significantly boosts the exciton binding energy by the dielectric confinement effect. Through the Förster resonance energy transfer, the excitons down-convert and recombine radiatively in the 2D perovskites. We report room-temperature pure green (n = 7-10), sky blue (n = 5), pure blue (n = 3), and deep blue (n = 1) electroluminescence, with record-high external quantum efficiencies in the green-to-blue wavelength region.

Ionic Effects on Supercritical CO2-Brine Interfacial Tensions: Molecular Dynamics Simulations and a Universal Correlation with Ionic Strength, Temperature, and Pressure.

For geological CO2 storage in deep saline aquifers, the interfacial tension (IFT) between supercritical CO2 and brine is critical for the storage security and design of the storage capacitance. However, currently, no predictive model exists to determine the IFT of supercritical CO2 against complex electrolyte solutions involving various mixed salt species at different concentrations and compositions. In this paper, we use molecular dynamics (MD) simulations to investigate the effect of salt ions on the incremental IFT at the supercritical CO2-brine interface with respect to that at the reference supercritical CO2-water interface. Supercritical CO2-NaCl solution, CO2-CaCl2 solution and CO2-(NaCl+CaCl2) mixed solution systems are simulated at 343 K and 20 MPa under different salinities and salt compositions. We find that the valence of the cations is the primary contributor to the variation in IFT, while the Lennard-Jones potentials for the cations pose a smaller impact on the IFT. Interestingly, the incremental IFT exhibits a general linear correlation with the ionic strength in the above three electrolyte systems, and the slopes are almost identical and independent of the solution types. Based on this finding, a universal predictive formula for IFTs of CO2-complex electrolyte solution systems is established, as a function of ionic strength, temperature, and pressure. The predicted IFTs using the established formula agree perfectly (with a high statistical confidence level of ∼96%) with a wide range of experimental data for CO2 interfacing with different electrolyte solutions, such as those involving MgCl2 and Na2SO4. This work provides an efficient and accurate route to directly predict IFTs in supercritical CO2-complex electrolyte solution systems for practical engineering applications, such as geological CO2 sequestration in deep saline aquifers and other interfacial systems involving complex electrolyte solutions.