A segmented basis set of quadruple zeta valence quality plus polarization functions(QZP)for H through Xe was developed to be used in conjunction with the ZORA Hamiltonian.This set was augmented with diffuse functions ...A segmented basis set of quadruple zeta valence quality plus polarization functions(QZP)for H through Xe was developed to be used in conjunction with the ZORA Hamiltonian.This set was augmented with diffuse functions to describe electrons farther away from the nuclei adequately.Using the ZORA-CCSD(T)/QZP-ZORA theoretical model,atomic ionization energies and bond lengths,harmonic vibrational frequencies,and atomization energies of some molecules were calculated.The addition of core-valence corrections has been shown to improve the agreement between theoretical and experimental results for molecular properties.For atomization energies,a similar observation emerges when considering spin-orbit couplings.With the augmented QZP-ZORA set,static mean dipole polarizabilities of a set of atoms were calculated and compared with previously published recommended and experimental values.Performance evaluations of the ZORA and Douglas–Kroll–Hess Hamiltonians were made for each property studied.展开更多
A DFT conformational and vibrational analysis of a single molecule of cisplatin (cis-[Pt(NH3)2Cl2]) was performed by means of PW91 functional and LANL08 ECP basis set for the Pt atom. 3-21G and 3-21G* Basis sets were ...A DFT conformational and vibrational analysis of a single molecule of cisplatin (cis-[Pt(NH3)2Cl2]) was performed by means of PW91 functional and LANL08 ECP basis set for the Pt atom. 3-21G and 3-21G* Basis sets were used for the remaining atoms. All the initially chosen conformations were found to converge to the global minimum conformation of C2v symmetry with H atoms lying in the coordination plane and pointing to the Cl atoms. The computational results were compared with the newest experimental structural data and with the vibrational spectroscopic data for cisplatin, obtained by other workers. The chosen level of theory was found to describe satisfactory the molecular structure (r. m. s. of the relative deviations ≤ 6%) and the harmonic vibrational frequencies (r. m. s. of the relative deviations ≤ 5%) of cisplatin.展开更多
Segmented all-electron basis set of triple zeta valence quality plus polarization functions(TZP)for the elements of the fifth row to be used together with the zero-order regular approximation(ZORA)is carefully constru...Segmented all-electron basis set of triple zeta valence quality plus polarization functions(TZP)for the elements of the fifth row to be used together with the zero-order regular approximation(ZORA)is carefully constructed.To correctly describe electrons distant from atomic nuclei,the basis set is augmented with diffuse functions giving rise to a set designated as ATZP-ZORA.At the ZORA-B3LYP theoretical level,these sets are used to calculate the ionization energy and mean dipole polarizability of some atoms,bond length,dissociation energy,and harmonic vibrational frequency of diatomic molecules.Then,these results are compared with the theoretical and experimental data found in the literature.Even considering that our sets are relatively compact,they are sufficiently accurate and reliable to perform property calculations involving simultaneously electrons from the inner shell and outer shell.The performances of the ZORA and second-order Douglas-Kroll-Hess Hamiltonians are evaluated and the results are also discussed.展开更多
2-D and 3-D micro-architectured multiphase thermoelastic metamaterials are designed and analyzed using a parametric level set method for topology optimization and the finite element method.An asymptotic homogenization...2-D and 3-D micro-architectured multiphase thermoelastic metamaterials are designed and analyzed using a parametric level set method for topology optimization and the finite element method.An asymptotic homogenization approach is employed to obtain the effective thermoelastic properties of the multiphase metamaterials.Theε-constraint multi-objective optimization method is adopted in the formulation.The coefficient of thermal expansion(CTE)and Poisson’s ratio(PR)are chosen as two objective functions,with the CTE optimized and the PR treated as a constraint.The optimization problems are solved by using the method of moving asymptotes.Effective isotropic and anisotropic CTEs and stiffness constants are obtained for the topologically optimized metamaterials with prescribed values of PR under the constraints of specified effective bulk modulus,volume fractions and material symmetry.Two solid materials along with one additional void phase are involved in each of the 2-D and 3-D optimal design examples.The numerical results reveal that the newly proposed approach can integrate shape and topology optimizations and lead to optimal microstructures with distinct topological boundaries.The current method can topologically optimize metamaterials with a positive,negative or zero CTE and a positive,negative or zero Poisson’s ratio.展开更多
Density functional theory method has been employed to investigate the structures of the prototypical technetium-labeled diphosphonate complex 99mTc-MDP, where MDP represents methylenediphosphonic acid. A total of 14 t...Density functional theory method has been employed to investigate the structures of the prototypical technetium-labeled diphosphonate complex 99mTc-MDP, where MDP represents methylenediphosphonic acid. A total of 14 trial structures were generated by allowing for the geometric, conformational, charge, and spin isomerism. Based on the optimized structures and calculated energies at the B3LYP/LANL2DZ level, two stable isomers were determined for the title complex. And they were further studied systematically in comparison with the experimental structure. The basis sets 6-31G*(LANL2DZ for Tc), 6-31G*(cc-pVDZ-pp for Tc), and DGDZVP have also been employed in combination with the B3LYP functional to study the basis set effect on the geometries of isomers. The optimized structures agree well with the available experimental data, and the bond lengths are more sensitive to the basis set than the bond angles. The charge distributions were studied by the Mulliken population analysis and natural bond orbital analysis. The results reflect a significant ligand-to-metal electron donation.展开更多
A finite element method with boundary element method (FEM-BEM) is presented for computing electromagnetic induction. The features of an edge element method including the volume and surface edge element method are inve...A finite element method with boundary element method (FEM-BEM) is presented for computing electromagnetic induction. The features of an edge element method including the volume and surface edge element method are investigated in depth. Surface basis functions of edge elements to an arbitrary shape of target are derived according to the geometrical property of basis functions and applied to discretize the surface integral equation for 3-D general targets. The proposed model is presented to compute resonant frequencies and surface current of underground unexplored ordnance (UXO), and then the electromagnetic responses of single target with different frequencies and positions of sensor are simulated and results are validated by experiments.展开更多
A novel numerical method for eliminating the singular integral and boundary effect is processed. In the proposed method, the virtual boundaries corresponding to the numbers of the true boundary arguments are chosen to...A novel numerical method for eliminating the singular integral and boundary effect is processed. In the proposed method, the virtual boundaries corresponding to the numbers of the true boundary arguments are chosen to be as simple as possible. An indirect radial basis function network (IRBFN) constructed by functions resulting from the indeterminate integral is used to construct the approaching virtual source functions distributed along the virtual boundaries. By using the linear superposition method, the governing equations presented in the boundaries integral equations (BIE) can be established while the fundamental solutions to the problems are introduced. The singular value decomposition (SVD) method is used to solve the governing equations since an optimal solution in the least squares sense to the system equations is available. In addition, no elements are required, and the boundary conditions can be imposed easily because of the Kronecker delta function properties of the approaching functions. Three classical 2D elasticity problems have been examined to verify the performance of the method proposed. The results show that this method has faster convergence and higher accuracy than the conventional boundary type numerical methods.展开更多
A non-orthogonal predefined exponential basis set is used to handle half-bounded domains in multi domain spectral method (MDSM). This approach works extremely well for real-valued semi-infinite differential problems. ...A non-orthogonal predefined exponential basis set is used to handle half-bounded domains in multi domain spectral method (MDSM). This approach works extremely well for real-valued semi-infinite differential problems. It spans simultaneously wide range of exponential decay rates with multi scaling and does not suffer from zero crossing. These two conditions are necessary for many physical problems. For comparison, the method is used to solve different problems and compared with analytical and published results. The comparison exhibits the strengths and accuracy of the presented basis set.展开更多
The interactions between metal ions such as Zn2+, Pb2+, Mn2+, Hg2+, Cd2+, Ni2+ and chitosan have been investigated using the model cluster model method and density functional method. Full optimization and frequency an...The interactions between metal ions such as Zn2+, Pb2+, Mn2+, Hg2+, Cd2+, Ni2+ and chitosan have been investigated using the model cluster model method and density functional method. Full optimization and frequency analysis of all cluster models have been performed employing B3LYP hybrid method at 3-21G basis set level except metal ions which were invoked to use effective core potential (ECP) method. The energy changes, and the main structural parameters have been obtained during the theoretical study of the adsorption of metal ions on the chitosan. The calculations showed that the coordination modes of metal ions with chitosan models were different, the geometries of Mn2+, Zn2+, Cd2+, Hg2+, Pb2+ ions coordinated with two nitrogen atoms and two oxygen atoms were distorted tetrahedral, while the square planar structure of Ni2+ coordinated two nitrogen atoms and two oxygen atoms was observed. The heat of reaction between six metal ions and chitosan models showed the order: Mn2+ >Ni2+ >Zn2+ >Pb2+ >Hg2+ >Cd2+, this suggested that the coordination strength of Mn2+ >Ni2+ >Zn2+ >Pb2+ >Hg2+ >Cd2+.展开更多
In this study, the energy for the ground state of helium and a few helium-like ions (Z=1-6) is computed variationally by using a Hylleraas-like wavefunction. A four-parameters wavefunction, satisfying boundary condi...In this study, the energy for the ground state of helium and a few helium-like ions (Z=1-6) is computed variationally by using a Hylleraas-like wavefunction. A four-parameters wavefunction, satisfying boundary conditions for coalescence points, is combined with a Hylleraas-like basis set which explicitly incorporates r12 interelectronic distance. The main contribution of this work is the introduction of modified correlation terms leading to the definition of integral transforms which provide the calculation of expectation value of energy to be done analytically over single-particle coordinates instead of Hylleraas coordinates.展开更多
In this work we consider the problem of shape reconstruction from an unorganized data set which has many important applications in medical imaging, scientific computing, reverse engineering and geometric modelling. Th...In this work we consider the problem of shape reconstruction from an unorganized data set which has many important applications in medical imaging, scientific computing, reverse engineering and geometric modelling. The reconstructed surface is obtained by continuously deforming an initial surface following the Partial Differential Equation (PDE)-based diffusion model derived by a minimal volume-like variational formulation. The evolution is driven both by the distance from the data set and by the curvature analytically computed by it. The distance function is computed by implicit local interpolants defined in terms of radial basis functions. Space discretization of the PDE model is obtained by finite co-volume schemes and semi-implicit approach is used in time/scale. The use of a level set method for the numerical computation of the surface reconstruction allows us to handle complex geometry and even changing topology,without the need of user-interaction. Numerical examples demonstrate the ability of the proposed method to produce high quality reconstructions. Moreover, we show the effectiveness of the new approach to solve hole filling problems and Boolean operations between different data sets.展开更多
Many physical problems involve unbounded domains where the physical quantities vanish at infinities.Numerically,this has been handled using different techniques such as domain truncation,approximations using infinitel...Many physical problems involve unbounded domains where the physical quantities vanish at infinities.Numerically,this has been handled using different techniques such as domain truncation,approximations using infinitely extended and vanishing basis sets,and mapping bounded basis sets using some coordinate transformations.Each technique has its own advantages and disadvantages.Yet,approximating simultaneously and efficiently a wide range of decaying rates has persisted as major challenge.Also,coordinate transformation,if not carefully implemented,can result in non-orthogonal mapped basis sets.In this work,we revisited this issue with an emphasize on designing appropriate transformations using sine series as basis set.The transformations maintain both the orthogonality and the efficiency.Furthermore,using simple basis set(sine function)help avoid the expensive numerical integrations.In the calculations,four types of physically recurring decaying behaviors are considered,which are:non-oscillating and oscillating exponential decays,and non-oscillating and oscillating algebraic decays.The results and the analyses show that properly designed high-order mapped basis sets can be efficient tools to handle challenging physical problems on unbounded domains.Decay rate ranges as large of 6 orders of magnitudes can be approximated efficiently and concurrently.展开更多
基金the financial support of Conselho Nacional de Desenvolvimento Científico e Tecnológico and Coordenacao de Aperfeic oamento de Pessoal de Nível Superior (Brazilian Agencies)。
文摘A segmented basis set of quadruple zeta valence quality plus polarization functions(QZP)for H through Xe was developed to be used in conjunction with the ZORA Hamiltonian.This set was augmented with diffuse functions to describe electrons farther away from the nuclei adequately.Using the ZORA-CCSD(T)/QZP-ZORA theoretical model,atomic ionization energies and bond lengths,harmonic vibrational frequencies,and atomization energies of some molecules were calculated.The addition of core-valence corrections has been shown to improve the agreement between theoretical and experimental results for molecular properties.For atomization energies,a similar observation emerges when considering spin-orbit couplings.With the augmented QZP-ZORA set,static mean dipole polarizabilities of a set of atoms were calculated and compared with previously published recommended and experimental values.Performance evaluations of the ZORA and Douglas–Kroll–Hess Hamiltonians were made for each property studied.
文摘A DFT conformational and vibrational analysis of a single molecule of cisplatin (cis-[Pt(NH3)2Cl2]) was performed by means of PW91 functional and LANL08 ECP basis set for the Pt atom. 3-21G and 3-21G* Basis sets were used for the remaining atoms. All the initially chosen conformations were found to converge to the global minimum conformation of C2v symmetry with H atoms lying in the coordination plane and pointing to the Cl atoms. The computational results were compared with the newest experimental structural data and with the vibrational spectroscopic data for cisplatin, obtained by other workers. The chosen level of theory was found to describe satisfactory the molecular structure (r. m. s. of the relative deviations ≤ 6%) and the harmonic vibrational frequencies (r. m. s. of the relative deviations ≤ 5%) of cisplatin.
基金the Conselho Nacional de Desenvolvimento Científico Tecnológico(Brazilian Agency)。
文摘Segmented all-electron basis set of triple zeta valence quality plus polarization functions(TZP)for the elements of the fifth row to be used together with the zero-order regular approximation(ZORA)is carefully constructed.To correctly describe electrons distant from atomic nuclei,the basis set is augmented with diffuse functions giving rise to a set designated as ATZP-ZORA.At the ZORA-B3LYP theoretical level,these sets are used to calculate the ionization energy and mean dipole polarizability of some atoms,bond length,dissociation energy,and harmonic vibrational frequency of diatomic molecules.Then,these results are compared with the theoretical and experimental data found in the literature.Even considering that our sets are relatively compact,they are sufficiently accurate and reliable to perform property calculations involving simultaneously electrons from the inner shell and outer shell.The performances of the ZORA and second-order Douglas-Kroll-Hess Hamiltonians are evaluated and the results are also discussed.
文摘2-D and 3-D micro-architectured multiphase thermoelastic metamaterials are designed and analyzed using a parametric level set method for topology optimization and the finite element method.An asymptotic homogenization approach is employed to obtain the effective thermoelastic properties of the multiphase metamaterials.Theε-constraint multi-objective optimization method is adopted in the formulation.The coefficient of thermal expansion(CTE)and Poisson’s ratio(PR)are chosen as two objective functions,with the CTE optimized and the PR treated as a constraint.The optimization problems are solved by using the method of moving asymptotes.Effective isotropic and anisotropic CTEs and stiffness constants are obtained for the topologically optimized metamaterials with prescribed values of PR under the constraints of specified effective bulk modulus,volume fractions and material symmetry.Two solid materials along with one additional void phase are involved in each of the 2-D and 3-D optimal design examples.The numerical results reveal that the newly proposed approach can integrate shape and topology optimizations and lead to optimal microstructures with distinct topological boundaries.The current method can topologically optimize metamaterials with a positive,negative or zero CTE and a positive,negative or zero Poisson’s ratio.
基金This work was supported by the National Natural Science Foundation of China (No.20801024 and No.21001055), the Natural Science Foundation of Jiangsu Province (No.BK2009077), and the Science Foundation of Health Department of Jiangsu Province (No.H200963).
文摘Density functional theory method has been employed to investigate the structures of the prototypical technetium-labeled diphosphonate complex 99mTc-MDP, where MDP represents methylenediphosphonic acid. A total of 14 trial structures were generated by allowing for the geometric, conformational, charge, and spin isomerism. Based on the optimized structures and calculated energies at the B3LYP/LANL2DZ level, two stable isomers were determined for the title complex. And they were further studied systematically in comparison with the experimental structure. The basis sets 6-31G*(LANL2DZ for Tc), 6-31G*(cc-pVDZ-pp for Tc), and DGDZVP have also been employed in combination with the B3LYP functional to study the basis set effect on the geometries of isomers. The optimized structures agree well with the available experimental data, and the bond lengths are more sensitive to the basis set than the bond angles. The charge distributions were studied by the Mulliken population analysis and natural bond orbital analysis. The results reflect a significant ligand-to-metal electron donation.
文摘A finite element method with boundary element method (FEM-BEM) is presented for computing electromagnetic induction. The features of an edge element method including the volume and surface edge element method are investigated in depth. Surface basis functions of edge elements to an arbitrary shape of target are derived according to the geometrical property of basis functions and applied to discretize the surface integral equation for 3-D general targets. The proposed model is presented to compute resonant frequencies and surface current of underground unexplored ordnance (UXO), and then the electromagnetic responses of single target with different frequencies and positions of sensor are simulated and results are validated by experiments.
文摘A novel numerical method for eliminating the singular integral and boundary effect is processed. In the proposed method, the virtual boundaries corresponding to the numbers of the true boundary arguments are chosen to be as simple as possible. An indirect radial basis function network (IRBFN) constructed by functions resulting from the indeterminate integral is used to construct the approaching virtual source functions distributed along the virtual boundaries. By using the linear superposition method, the governing equations presented in the boundaries integral equations (BIE) can be established while the fundamental solutions to the problems are introduced. The singular value decomposition (SVD) method is used to solve the governing equations since an optimal solution in the least squares sense to the system equations is available. In addition, no elements are required, and the boundary conditions can be imposed easily because of the Kronecker delta function properties of the approaching functions. Three classical 2D elasticity problems have been examined to verify the performance of the method proposed. The results show that this method has faster convergence and higher accuracy than the conventional boundary type numerical methods.
文摘A non-orthogonal predefined exponential basis set is used to handle half-bounded domains in multi domain spectral method (MDSM). This approach works extremely well for real-valued semi-infinite differential problems. It spans simultaneously wide range of exponential decay rates with multi scaling and does not suffer from zero crossing. These two conditions are necessary for many physical problems. For comparison, the method is used to solve different problems and compared with analytical and published results. The comparison exhibits the strengths and accuracy of the presented basis set.
基金The Ph. D. Foundation (Y030426) and Post-doctoral Foundation of China University of Petroleum (East China).
文摘The interactions between metal ions such as Zn2+, Pb2+, Mn2+, Hg2+, Cd2+, Ni2+ and chitosan have been investigated using the model cluster model method and density functional method. Full optimization and frequency analysis of all cluster models have been performed employing B3LYP hybrid method at 3-21G basis set level except metal ions which were invoked to use effective core potential (ECP) method. The energy changes, and the main structural parameters have been obtained during the theoretical study of the adsorption of metal ions on the chitosan. The calculations showed that the coordination modes of metal ions with chitosan models were different, the geometries of Mn2+, Zn2+, Cd2+, Hg2+, Pb2+ ions coordinated with two nitrogen atoms and two oxygen atoms were distorted tetrahedral, while the square planar structure of Ni2+ coordinated two nitrogen atoms and two oxygen atoms was observed. The heat of reaction between six metal ions and chitosan models showed the order: Mn2+ >Ni2+ >Zn2+ >Pb2+ >Hg2+ >Cd2+, this suggested that the coordination strength of Mn2+ >Ni2+ >Zn2+ >Pb2+ >Hg2+ >Cd2+.
基金Project supported by the Cumhuriyet University National MOVPE Crystal Growth and Characterization Laboratory,DPT-K120,TUBITAK (Grant Nos TBAG 105T492,TBAG 107T012,and TBAG-108T015)
文摘In this study, the energy for the ground state of helium and a few helium-like ions (Z=1-6) is computed variationally by using a Hylleraas-like wavefunction. A four-parameters wavefunction, satisfying boundary conditions for coalescence points, is combined with a Hylleraas-like basis set which explicitly incorporates r12 interelectronic distance. The main contribution of this work is the introduction of modified correlation terms leading to the definition of integral transforms which provide the calculation of expectation value of energy to be done analytically over single-particle coordinates instead of Hylleraas coordinates.
基金supported by PRIN-MIUR-Cofin 2006,project,by"Progetti Strategici EF2006"University of Bologna,and by University of Bologna"Funds for selected research topics"
文摘In this work we consider the problem of shape reconstruction from an unorganized data set which has many important applications in medical imaging, scientific computing, reverse engineering and geometric modelling. The reconstructed surface is obtained by continuously deforming an initial surface following the Partial Differential Equation (PDE)-based diffusion model derived by a minimal volume-like variational formulation. The evolution is driven both by the distance from the data set and by the curvature analytically computed by it. The distance function is computed by implicit local interpolants defined in terms of radial basis functions. Space discretization of the PDE model is obtained by finite co-volume schemes and semi-implicit approach is used in time/scale. The use of a level set method for the numerical computation of the surface reconstruction allows us to handle complex geometry and even changing topology,without the need of user-interaction. Numerical examples demonstrate the ability of the proposed method to produce high quality reconstructions. Moreover, we show the effectiveness of the new approach to solve hole filling problems and Boolean operations between different data sets.
基金The author would like to thank Qatar National Research Fund as the work is partially supported by NPRP 7-317-1-055.
文摘Many physical problems involve unbounded domains where the physical quantities vanish at infinities.Numerically,this has been handled using different techniques such as domain truncation,approximations using infinitely extended and vanishing basis sets,and mapping bounded basis sets using some coordinate transformations.Each technique has its own advantages and disadvantages.Yet,approximating simultaneously and efficiently a wide range of decaying rates has persisted as major challenge.Also,coordinate transformation,if not carefully implemented,can result in non-orthogonal mapped basis sets.In this work,we revisited this issue with an emphasize on designing appropriate transformations using sine series as basis set.The transformations maintain both the orthogonality and the efficiency.Furthermore,using simple basis set(sine function)help avoid the expensive numerical integrations.In the calculations,four types of physically recurring decaying behaviors are considered,which are:non-oscillating and oscillating exponential decays,and non-oscillating and oscillating algebraic decays.The results and the analyses show that properly designed high-order mapped basis sets can be efficient tools to handle challenging physical problems on unbounded domains.Decay rate ranges as large of 6 orders of magnitudes can be approximated efficiently and concurrently.