Valence Basis Sets Research Articles

Synthesis of tamsulosin precursor through palladium nanocatalytic hydrogenolysis: Crystal structure, DFT analysis, and ADME studies

The present work describes a potential method for the synthesis of industrially important Tamsulosine intermediate under milder reaction conditions. The major advantages include i) nano catalytic system and ii) ammonium formate is used as greener alternative during hydrogenolysis with an overall yield of 98 %. Density functional theory (DFT) studies in combination with B3LYP exchange-correlation functional method with split valence basis set 6–311 G is used to understand molecular structure of intermediate is also described. The crystal structure of tamsulosine precursor is elucidated by using single crystal x-ray crystallography is also reported. Therefore, this study is expected to bridge the need to optimize a greener alternative process in this hydrogenolysis reaction to synthesize tamsulosin precursor. ADME studies reveals that better oral bioavailability, aqueous solubility, Gastrointerstinal absorption and option “yes” in Lipinski, Ghose, Veber, Egan and Muegge rules.

Exchange-correlation kernel for perturbation dependent auxiliary functions in auxiliary density perturbation theory

ContextAnalytic exchange-correlation kernel formulations are of the outermost importance for density functional theory (DFT) perturbation calculations. In this paper, the working equation for the exchange-correlation kernel of the generalized gradient approximation (GGA) for perturbation dependent auxiliary functions is derived and discussed in the framework of auxiliary density functional theory (ADFT). The presented new formulation is extended to the unrestricted approach, too. A comprehensive discussion of the implementation of the GGA ADFT kernel, using either the native exchange-correlation functional implementations in deMon2k or the ones from the LibXC library, is given. Calculations with analytic exchange-correlation kernels are compared to their finite difference counterparts. The obtained results are in quantitative agreement. Nevertheless, analytic GGA ADFT kernel implementations show substantial improvement in the computational performance. Similar results are reported for analytic second derivatives of effective core potential (ECP) and model core potential (MCP) matrix elements when compared to their finite difference counterparts in molecular frequency analyses.MethodAll calculations are performed in the framework of ADFT as implemented in deMon2k. In the ADFT analytic frequency calculations, auxiliary density perturbation theory was used. The underlying two-center exchange-correlation kernel matrix elements are calculated by numerical integration either with analytic or finite difference kernel expressions. Validation calculations are performed with the VWN and PBE functionals employing DFT-optimized DZVP basis sets in conjunction with automatically generated GEN-A2 auxiliary density function sets. In the (Pt3Cu)n cluster benchmark calculations, the RPBE functional was used. For Pt atoms, the quasi-relativistic LANL2DZ effective core potential with the corresponding valence basis set was employed, whereas for Cu atoms, the all-electron DFT-optimized TZVP basis was applied. The auxiliary density was expanded by the automatically generated GEN-A2* auxiliary function set. We run all benchmark calculations in parallel on 24 cores.

Quantum-Chemical Simulation of 13C NMR Chemical Shifts of Fullerene C60 Exo-Derivatives

The 13C NMR chemical shifts of fullerene C60 exo-derivatives were calculated using quantum chemical hybrid functionals combined with Pople, Dunning correlation-consistent, and def2-TZVP split valence basis sets taking into account the solvent effect (polarizable continuum model). A relationship between theoretical and experimental 13C NMR chemical shifts (CSs) is assessed quantitatively to select a functional/basis set combination. It is found that the CAM-B3LYP/6-31G and M06L/6-31G combinations have the best convergence with experimental data in modeling the 13С NMR CSs of sp3 fullerene carbon atoms in С60 derivatives, whereas X3LYP/6-31G and CAM-B3LYP/6-31G(d) in modeling the 13С NMR CSs of their sp2 fullerene carbon atoms.

Probing the Potential Energy Profile of the I + (H2O)3 → HI + (H2O)2OH Forward and Reverse Reactions: High Level CCSD(T) Studies with Spin-Orbit Coupling Included.

Three different pathways for the atomic iodine plus water trimer reaction I + (H2O)3 → HI + (H2O)2OH were preliminarily examined by the DFT-MPW1K method. Related to previous predictions for the F/Cl/Br + (H2O)3 reactions, three pathways for the I + (H2O)3 reaction are linked in terms of geometry and energetics. To legitimize the results, the "gold standard" CCSD(T) method was employed to investigate the lowest-lying pathway with the correlation-consistent polarized valence basis set up to cc-pVQZ(-PP). According to the CCSD(T)/cc-pVQZ(-PP)//CCSD(T)/cc-pVTZ(-PP) results, the I + (H2O)3 → HI + (H2O)2OH reaction is predicted to be endothermic by 47.0 kcal mol-1. The submerged transition state is predicted to lie 43.7 kcal mol-1 above the separated reactants. The I···(H2O)3 entrance complex lies below the separated reactants by 4.1 kcal mol-1, and spin-orbit coupling has a significant impact on this dissociation energy. The HI···(H2O)2OH exit complex is bound by 4.3 kcal mol-1 in relation to the separated products. Compared with simpler I + (H2O)2 and I + H2O reactions, the I + (H2O)3 reaction is energetically between them in general. It is speculated that the reaction between the iodine atom and the larger water clusters may be energetically analogous to the I + (H2O)3 reaction. The iodine reaction I + (H2O)3 is connected with the analogous valence isoelectronic bromine/chlorine reactions Br/Cl + (H2O)3 but much different from the F + (H2O)3 reaction. Significant difference with other halogen systems, especially for barrier heights, are seen for the iodine systems.

In-silico docking, synthesis, structure analysis, DFT calculations, energy frameworks, and pharmacological intervention of [1,3,4]-thiadiazoles analogous as XO inhibitor and on multiple molecular inflammatory targets COX and LOX

Xanthine oxidase (XO) is an interesting target for the synergic treatment of several diseases such as gout, hypertension, type 2 diabetes, and kidney disease. Associated complication includes hyperuricemia, which is considered one of the most common metabolic lifestyle disorders worldwide. Overproduction and insufficient excretion of uric acid during purine metabolism results in the formation of uric acid crystals in kidneys, joints, and other tissues. Allopurinol and Febuxostat, a clinically available drug used to reduce XO activity and serum uric acid level, are demonstrated to produce adverse effects including gastrointestinal distress, hypersensitivity, and skin rash, raising safety concerns. Exploring newer and novel XO inhibitors with minimum or no side effects is the aim of the current study. [1,3,4]-Thaidiazole analogs (7a-j) scaffold plays an important role in the design of efficient and potent inhibitors. In the present work, (7a-j) were evaluated for their ability to inhibit XO and found that among the series compound (7i) is selected as the best compound for inhibition of XO. Further, cyclooxygenases (COXs) and 5-lipoxygenase (5-LOX), enzymes that generate biologically active lipid molecules termed eicosanoids, are considered inflammatory and play a vital role in inflammatory pathways in the human body. Apart from their relation with inflammation, the additional involvement of COX-2 enzyme with cancer activity was recently discovered. In some cancer types, the level of COX-2 enzyme is increased indicating that this enzyme could be a suitable target for cancer therapy. Based on these findings, we have synthesized some new novel [1,3,4]-thiadiazole analogous (7a-j) and tested them against XO, COX, and LOX. To gain insight into the activity of XO, COX, and LOX inhibition, molecular docking studies were carried out for XO, COX, and LOX enzymes utilizing the newly synthesized compounds (7a-j), and the results inferred that among the synthesized series of [1,3,4]-thiadiazoles (7a-j), the compound (7i) showed extremely good result in both In silico and evaluation studies of XO, COX, and LOX, which is in agreement with the experimental results. Furthermore, density functional theory (DFT) was applied to calculate HOMO-LUMO, energy gap, and other parameters under the hybrid exchange-correlation functional (B3LYP) with the triple split valence basis set (6-311G).

Structural Stability and Electronic Properties of Boron Phosphide Nanotubes: A Density Functional Theory Perspective

Based on the Density Functional Theory (DFT) calculations, we analyze the structural and electronic properties of boron phosphide nanotubes (BPNTs) as functions of chirality. The DFT calculations are performed using the M06-2X method in conjunction with the 6-31G(d) divided valence basis set. All nanostructures, (n,0) BPNT (n = 5–8, 10, 12, 14) and (n,n) BPNT (n = 3–11), were optimized minimizing the total energy, assuming a non-magnetic nature and a total charge neutrality. Results show that the BPNT diameter size increases linearly with the chiral index “n” for both chiralities. According to the global molecular descriptors, the (3,3) BPNT is the most stable structure provided that it shows the largest global hardness value. The low chirality (5,0) BPNT has a strong electrophilic character, and it is the most conductive system due to the small |HOMO-LUMO| energy gap. The chemical potential and electrophilicity index in the zigzag-type BPNTs show remarkable chirality-dependent behavior. The increase in diameter/chirality causes a gradual decrease in the |HOMO-LUMO| energy gap for the zigzag BPNTs; however, in the armchair-type BPNTs, a phase transition is generated from a semiconductor to a conductor system. Therefore, the nanostructures investigated in this work may be suggested for both electrical and biophysical applications.

BSSE-corrected consistent Gaussian basis sets of triple-zeta valence with polarization quality of the fifth period for solid-state calculations.

Revised versions of our published pob-TZVP basis sets [Laun, J.; Vilela Oliveira, D. and Bredow, T., J. Comput. Chem., 2018, 39 (19), 1285-1290] have been derived for periodic quantum-chemical solid-state calculations. They complete our pob-TZVP-rev2 series [Vilela Oliveira, D.; Laun, J.; Peintinger, M. F. and Bredow, T., J. Comput. Chem., 2019, 40 (27), 2364-2376 and Laun, J. and Bredow, J. Comput. Chem. 2021; 42 (15), 1064-1072] for the elements of the fifth period and are based on the fully relativistic effective core potentials (ECPs) of the Stuttgart/Cologne group and the def2-TZVP valence basis of the Ahlrichs group. The pob-TZVP-rev2 basis sets are developed to minimize the basis set superposition error (BSSE) in crystalline systems. For the applied PW1PW hybrid functional, the overall performance, transferability, and SCF stability of the resulting pob-TZVP-rev2 basis sets are significantly improved compared to the original pob-TZVP basis sets. After augmentation with single diffuse s- and p-functions, reference plane-wave band structures of metals can be accurately reproduced.

Air-Stable Dy(III)-Macrocycle Enantiomers: From Chiral to Polar Space Group

Air-Stable Dy(III)-Macrocycle Enantiomers: From Chiral to Polar Space Group

Property-optimized Gaussian basis sets for lanthanides

Property-optimized Gaussian basis sets of split-valence, triple-zeta valence, and quadruple-zeta valence quality are developed for the lanthanides Ce-Lu for use with small-core relativistic effective core potentials. They are constructed in a systematic fashion by augmenting def2 orbital basis sets with diffuse basis functions and minimizing negative static isotropic polarizabilities of lanthanide atoms with respect to basis set exponents within the unrestricted Hartree-Fock method. The basis set quality is assessed using a test set of 70 molecules containing the lanthanides in their common oxidation states and f electron occupations. 5d orbital occupation turns out to be the determining factor for the basis set convergence of polarizabilities in lanthanide atoms and the molecular test set. Therefore, two series of property-optimized basis sets are defined. The augmented def2-SVPD, def2-TZVPPD, and def2-QZVPPD basis sets balance the accuracy of polarizabilities across lanthanide oxidation states. The relative errors in atomic and molecular polarizability calculations are ≤8% for augmented split-valence basis sets, ≤ 2.5% for augmented triple-zeta valence basis sets, and ≤1% for augmented quadruple-zeta valence basis sets. In addition, extended def2-TZVPPDD and def2-QZVPPDD are provided for accurate calculations of lanthanide atoms and neutral clusters. The property-optimized basis sets developed in this work are shown to accurately reproduce electronic absorption spectra of a series of LnCp3 '- complexes (Cp' = C5H4SiMe3, Ln = Ce-Nd, Sm) with time-dependent density functional theory.

An Aggregation-Induced Emission Optical Highlighter for the Studies of Endoplasmic Reticulum-Lipid Droplet Content Dynamics

An Aggregation-Induced Emission Optical Highlighter for the Studies of Endoplasmic Reticulum-Lipid Droplet Content Dynamics

A revisited study of the low-lying electronic states of HF molecule

Hydrogen fluoride (HF) has served as a benchmark molecule for ab initio methods. Consequently, it has been extensively investigated in the past few decades. The energy ordering of the low-lying triplet electronic states has proven controversial in previous computational studies. In this work, the Gaussian basis of F and H atoms are calibrated by adding diffuse functions to the aug-cc-pV5Z valence basis sets. The electronic structure of HF molecule is studied by the multireference configuration interaction method plus Davidson correction method (MRCI+Q). The potential energy curves (PECs) of 19 electronic states are computed, and the spectroscopic constants of the bound states are determined. The low-lying triplet Rydberg electronic states 23Σ+ (3Σ+R1), 23Π (3ΠR1) and 13Δ (3ΔR1) are characterized and compared with the previous computational results. The transition dipole moments (TDMs) and Franck-Condon factors (FCFs) between the electronic states are subsequently calculated. Finally, the spontaneous radiative lifetimes of four spin-allowed transitions are obtained. The results of this paper will improve the understanding on the electronic excited states of HF molecule, especially the electronic structure and spectra characteristics of the triplet excited states.

BSSE‐corrected consistent Gaussian basis sets of triple‐zeta valence with polarization quality of the sixth period for solid‐state calculations

Consistent basis sets of triple-zeta valence with polarization quality for the elements Cs-Po were derived for periodic quantum-chemical solid-state calculations. They are an extension of the pob-TZVP-rev2 [Vilela Oliveira, D.; Laun, J.; Peintinger, M. F. and Bredow, T., J. Comput. Chem., 2019, 40 (27), 2364-2376] basis sets and are based on the fully relativistic effective core potentials (ECPs) of the Stuttgart/Cologne group and on the def2-TZVP valence basis of the Ahlrichs group. The basis sets are constructed to minimize the basis set superposition error (BSSE) in crystalline systems. The contraction scheme, the orbital exponents, and contraction coefficients were optimized in order to ensure robust and stable self-consistent-field (SCF) convergence for a set of compounds and metals. For the applied PW1PW hybrid functional, the average deviations of the calculated lattice constants from experimental references are smaller with pob-TZVP-rev2 than with standard basis sets available from the CRYSTAL basis set database.

Effect of varying the TD-lc-DFTB range-separation parameter on charge and energy transfer in a model pentacene/buckminsterfullerene heterojunction.

Atomistic modeling of energy and charge transfer at the heterojunction of organic solar cells is an active field with many remaining outstanding questions owing, in part, to the difficulties in performing reliable photodynamics calculations on very large systems. One approach to being able to overcome these difficulties is to design and apply an appropriate simplified method. Density-functional tight binding (DFTB) has become a popular form of approximate density-functional theory based on a minimal valence basis set and neglect of all but two center integrals. We report the results of our tests of a recent long-range correction (lc) [A. Humeniuk and R. Mitrić, J. Chem. Phys. 143, 134120 (2015)] for time-dependent (TD) lc-DFTB by carrying out TD-lc-DFTB fewest switches surface hopping calculations of energy and charge transfer times using the relatively new DFTBABY [A. Humeniuk and R. Mitrić, Comput. Phys. Commun. 221, 174 (2017)] program. An advantage of this method is the ability to run enough trajectories to get meaningful ensemble averages. Our interest in the present work is less in determining exact energy and charge transfer rates than in understanding how the results of these calculations vary with the value of the range-separation parameter (Rlc = 1/μ) for a model organic solar cell heterojunction consisting of a gas-phase van der Waals complex P/F made up of a single pentacene (P) molecule together with a single buckminsterfullerene (F) molecule. The default value of Rlc = 3.03 a0 is found to be much too small as neither energy nor charge transfer is observed until Rlc ≈ 10 a0. Tests at a single geometry show that the best agreement with high-quality ab initio spectra is obtained in the limit of no lc (i.e., very large Rlc). A plot of energy and charge transfer rates as a function of Rlc is provided, which suggests that a value of Rlc ≈ 15 a0 yields the typical literature (condensed-phase) charge transfer time of about 100 fs. However, energy and charge transfer times become as high as ∼300 fs for Rlc ≈ 25 a0. A closer examination of the charge transfer process P*/F → P+/F- shows that the initial electron transfer is accompanied by a partial delocalization of the P hole onto F, which then relocalizes back onto P, consistent with a polaron-like picture in which the nuclei relax to stabilize the resultant redistribution of charges.

In silico Evaluation of Cucurbit[6]uril as a Potential Detector for Cocaine and Its Adulterants Lidocaine, Caffeine, and Procaine

Illicit drugs and their trafficking require worldwide efforts in investigation, detection, and control. Colorimetric tests are often applied to identify drugs. Cocaine has some well-known adulterants that can provide a false positive response. Cucurbit[6]uril (CB[6]) has been suggested as a potential detector for cocaine and other illicit drugs. This work uses in silico methods to evaluate the use of CB[6] to detect cocaine and these interfering substances. More specifically, this work analyzes different possibilities of CB[6] complexation with cocaine, lidocaine, caffeine, and procaine and compares the results achieved for cocaine and its adulterants. Different methodologies were employed: quantum chemistry was investigated through DFT B3LYP/TZVP (density functional theory-Becke, three-parameter, Lee-Yang-Parr with triple zeta valence plus polarization basis set) and the semi-empirical methods Austin model 1 (AM1), parametric methods 3, 6, and 7 (PM3, PM6, PM7), and Recife model 1 (RM1). We used these methodologies intending to compare the reasonability and reproducibility of the results in the gas phase condition. Solvent influence was studied by molecular dynamics (MD) simulations. Results showed that CB[6] does not bind to these substances, as judged from the positive values of binding free energy obtained with all methods. DFT and MD were the most reliable methods whereas semiempirical ones were not reproductible in describing these systems. Results also showed that interactions are not specific, so CB[6] does not provide a good response for cocaine detection.

The electronic spectrum of AgI2. Ab initio benchmark CASSCF + Averaged Coupled Pair Functional studies on the ligand-field states including spin-orbit couplings

The X2Πg, 2Σg+ and 2Δg states of AgI2 have been studied through benchmark ab initio CASSCF + Averaged Coupled Pair Functional (ACPF) calculations using especially developed valence basis sets to study the transition energies and equilibrium geometries. The spin–orbit (SO) effects were included through the effective Hamiltonian formalism using the |ΛSΣ〉 ACPF energies as diagonal elements. AgI2 is a linear centrosymmetric molecule with a 2Πg ground state, in contradiction with ligand-field theory. The ACPF adiabatic excitation energies of the 2Σg+ and 2Δg states are 7344 and 26085 cm−1, respectively. A comparison of the present results with the previous ones for AgCl2 and AgBr2 shows a monotonic increase in the 2Πg-2Σg+ and 2Πg-2Δg transition energies due to a decrease of the ionic structures in the LF states as the atomic number of the halogen grows. The inclusion of the SO effects leads to a pure Ω = 3/2 X(2Πg) ground state, a Ω = 1/2 (83% 2Πg + 17% 2Σg+) A state, a Ω = 1/2 (17% 2Πg + 83% 2Σg+) B state, a Ω = 5/2 (2Δg) C state and a Ω = 3/2 (2Δg) D state. The X-A, X-B, X-C, X-D vertical transition energies are 4250, 11248, 27,127 and 30441 cm−1, respectively, nearly identical to the adiabatic Te values.

A Lattice Distortion Theory for Promotor Containing Clathrate Hydrates

A lattice distortion theory for promotor containing clathrate hydrates is formulated using the statistical thermodynamics based model of van der Waals and Platteeuw in association with the ab initio quantum mechanics to compute the cavity potentials. Despite of high degree of lattice distortion anticipated for large and polar molecules of liquid promotors, their variable lattice energy concept is unreported. With this intention, we estimate the lattice stabilization energy from spin-component scaled second order Møller-Plesset (SCS-MP2) perturbation theory applied with the augmented correlation-consistent polarized double zeta valence (aug-cc-pVDZ) basis set. Implementing this to compute cavity potential for different promotors, the reference properties of hydrates are harvested by regressing against the phase equilibrium conditions of their binary hydrates with methane. Our study confirms the exponential relation of reference chemical potential difference with van der Waals volume of the promotors. Moreover, using the excess Gibbs free energy theory, the higher order distortions for the multiple guests are captured. The proposed lattice distortion theory is attested with phase equilibrium conditions of eight promotors containing clathrate hydrate systems, namely propylene oxide, acetone, tetrahydrofuran, pyrrolidine, iso-butanaldehyde, cyclopentane, furan and thiophene, all having methane as a co-guest.

The Effects of Split Valence Basis Sets on Muon Hyperfine Interaction in Guanine Nucleobase and Nucleotide Structures

The DFT cluster method was employed to investigate the electronic structures and muonium hyperfine interactions in guanine nucleobase and nucleotide using three different basis sets. The total energy and Fermi contact values were calculated for muon trapped at carbon '8'. The three basis sets, 6-31G, 6-311G and 6-311G(d,p), were used in tandem with the B3LYP functional. There are significant quantitative differences in the calculated total energy. 6-311G(d,p) produced the lowest total energy as compared to the other basis sets. The lowering of the total energy is due to the increase in the number of basis functions to describe the atomic orbitals, which is consistent with the postulate on basis set completeness. The 6-31G basis set produced the muon Fermi contact value that is the closest to the experimental value. The calculated Fermi contact values for the nucleobase and nucleotide are significantly lowered in going from the double-zeta to the triple-zeta basis set by 5% and 4% respectively. The lowering of the Fermi contact value can be attributed to the extension of the triple-zeta basis set in describing the valence atomic orbitals. The presence of the sugar phosphate group in the nucleotide instead of the methyl group tends to lower the Fermi contact value. Thus, the sugar phosphate group should be taken into consideration when designing a calculation model.

Combined experimental and theoretical investigations on a half-sandwich organometallic Os(II) complex

The structural, spectral (FT-IR, Raman, 1H and 13C NMR chemical shift and UV–Vis. spectroscopies) and electronic properties for bromido(η6-1-isopropyl-4-methylbenzene)(7-oxocyclohepta-1,3,5-trien-1-olato-κ2O,Oʹ)osmium(II) were investigated using experimental and theoretical methods. The bond parameters and aromatic π-π stacking interactions in molecular structure of the half-sandwich organometallic Os(II)-arene complex under investigation were studied. Hirshfeld surface analysis was used to assess C-H⋯O and C-H⋯Br intermolecular interactions in the crystal packing. To support experimental vibrational wavenumbers, structural parameters, NMR chemical shifts and electronic absorption spectroscopic data, theoretical computations were performed with Density functional theory (DFT) using the B3LYP method with a double-zeta valence basis set (LanL2DZ). DFT was also applied to explore the HOMO, LUMO, MEP, NBO and NLO of the Os(II) complex compound at the mentioned computational level. Good agreement was found between the calculated results and experimental data.

Investigation on the electrical and optical properties of some zinc titanate ceramics

The ceramics of spinel zinc orthotitanate (SZT) and perovskite zinc metatitanate (PZT) sintered at $$800^{\circ }\hbox {C}$$ have been prepared using the high-energy ball milling technique. The structural and optical characterisation of the prepared samples has been performed using X-ray diffraction and Fourier transformation of infrared (FT-IR) spectroscopy analysis. The measurements of frequency-dependent dielectric constant of the PZT and SZT ceramics sintered at $$800^{\circ }\hbox {C}$$ have been performed in the $$(50{-}600)^{\circ }\hbox {C}$$ temperature range. Impedance spectroscopy studies have been reported for the prepared nanocrystalline ceramics. The empirical vibrational frequencies observed from the FT-IR spectra have been compared using the computational method of ORCA program package. Constraints for symmetry, bonds, angles or dihedral angles have not been applied in the geometry optimisation calculations. Hybrid three-parameter exchange-correlation functional of Becke, Lee, Yang, Parr with 20% amount of exact exchange and Ahlrichs triple-zeta valence (def2-TZVP) basis set with polarisation function have been investigated for all atoms without considering the relativistic effects.

Segmented all-electron basis sets of triple zeta quality for the lanthanides: application to structure calculations of lanthanide monoxides.

Nonrelativistic and relativistic (Douglas-Kroll-Hess, DKH) segmented all-electron Gaussian basis sets of valence triple zeta quality plus polarization functions (TZP) for the lanthanides were developed. As some atomic and molecular properties depend on a good description of the electrons far from the nuclei, these basis sets are augmented with diffuse functions, giving rise to the augmented TZP (ATZP) and ATZP-DKH basis sets. At the DKH level of theory, the B3LYP hybrid functional in conjunction with the TZP-DKH basis set were used to calculate the atomic charges and valence orbital populations of the lanthanide and oxygen atoms, the bond lengths, and the equilibrium dissociation energies of lanthanide monoxides. The DKH-B3LYP/ATZP-DKH polarizability of Yb and the DKH-M06/TZP-DKH first ionization energies of the lanthanides are also reported. Compared with the values obtained with a larger all-electron basis set, and theoretical and experimental data found in the literature, data obtained by our compact basis sets are verified to be accurate and reliable. Unlike effective core potential valence basis sets, our basis sets can also be employed in molecular property calculations that involve the simultaneous treatment of core and valence electrons. Graphical abstract ᅟ.