10
0
mirror of https://github.com/LCPQ/quantum_package synced 2024-06-26 07:02:14 +02:00
quantum_package/src/Bielec_integrals/README.rst

184 lines
8.8 KiB
ReStructuredText
Raw Normal View History

2014-04-17 23:50:51 +02:00
=============
BiInts Module
=============
Here, all bi-electronic integrals (:math:`1/r_{12}`) are computed. As they have
2014-04-18 00:17:22 +02:00
4 indices and many are zero, they are stored in a map, as defined in
``Utils/map_module.f90``. To fetch an AO integral, use the
``get_ao_bielec_integral(i,j,k,l,ao_integrals_map)`` function, and to fetch and
MO integral, use ``get_mo_bielec_integral(i,j,k,l,mo_integrals_map)`` or
``mo_bielec_integral(i,j,k,l)``.
2014-04-17 23:50:51 +02:00
Needed Modules
==============
.. Do not edit this section. It was auto-generated from the
.. NEEDED_MODULES file.
* `AOs <http://github.com/LCPQ/quantum_package/tree/master/src/AOs>`_
* `Bitmask <http://github.com/LCPQ/quantum_package/tree/master/src/Bitmask>`_
* `Electrons <http://github.com/LCPQ/quantum_package/tree/master/src/Electrons>`_
* `Ezfio_files <http://github.com/LCPQ/quantum_package/tree/master/src/Ezfio_files>`_
* `MOs <http://github.com/LCPQ/quantum_package/tree/master/src/MOs>`_
* `Nuclei <http://github.com/LCPQ/quantum_package/tree/master/src/Nuclei>`_
* `Output <http://github.com/LCPQ/quantum_package/tree/master/src/Output>`_
* `Utils <http://github.com/LCPQ/quantum_package/tree/master/src/Utils>`_
* `MonoInts <http://github.com/LCPQ/quantum_package/tree/master/src/MonoInts>`_
2014-05-13 13:57:58 +02:00
Documentation
=============
.. Do not edit this section. It was auto-generated from the
.. NEEDED_MODULES file.
2014-05-14 00:01:31 +02:00
`ao_bielec_integral <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/ao_bi_integrals.irp.f#L1>`_
2014-05-13 13:57:58 +02:00
integral of the AO basis <ik|jl> or (ij|kl)
i(r1) j(r1) 1/r12 k(r2) l(r2)
2014-07-14 18:52:32 +02:00
`ao_bielec_integral_schwartz <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/ao_bi_integrals.irp.f#L350>`_
2014-05-13 13:57:58 +02:00
Needed to compuet Schwartz inequalities
2014-06-25 14:58:58 +02:00
`ao_bielec_integrals_in_map <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/ao_bi_integrals.irp.f#L188>`_
2014-05-13 13:57:58 +02:00
Map of Atomic integrals
i(r1) j(r2) 1/r12 k(r1) l(r2)
2014-06-25 14:58:58 +02:00
`compute_ao_bielec_integrals <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/ao_bi_integrals.irp.f#L147>`_
2014-05-13 13:57:58 +02:00
Compute AO 1/r12 integrals for all i and fixed j,k,l
2014-07-14 18:52:32 +02:00
`eri <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/ao_bi_integrals.irp.f#L511>`_
2014-05-13 13:57:58 +02:00
ATOMIC PRIMTIVE bielectronic integral between the 4 primitives ::
primitive_1 = x1**(a_x) y1**(a_y) z1**(a_z) exp(-alpha * r1**2)
primitive_2 = x1**(b_x) y1**(b_y) z1**(b_z) exp(- beta * r1**2)
primitive_3 = x2**(c_x) y2**(c_y) z2**(c_z) exp(-delta * r2**2)
primitive_4 = x2**(d_x) y2**(d_y) z2**(d_z) exp(- gama * r2**2)
2014-07-14 18:52:32 +02:00
`general_primitive_integral <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/ao_bi_integrals.irp.f#L376>`_
2014-05-13 13:57:58 +02:00
Computes the integral <pq|rs> where p,q,r,s are Gaussian primitives
2014-07-14 18:52:32 +02:00
`give_polynom_mult_center_x <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/ao_bi_integrals.irp.f#L656>`_
2014-05-13 13:57:58 +02:00
subroutine that returns the explicit polynom in term of the "t"
variable of the following polynomw :
I_x1(a_x, d_x,p,q) * I_x1(a_y, d_y,p,q) * I_x1(a_z, d_z,p,q)
2014-07-14 18:52:32 +02:00
`i_x1_new <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/ao_bi_integrals.irp.f#L600>`_
2014-05-13 13:57:58 +02:00
recursive function involved in the bielectronic integral
2014-07-14 18:52:32 +02:00
`i_x1_pol_mult <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/ao_bi_integrals.irp.f#L719>`_
2014-05-13 13:57:58 +02:00
recursive function involved in the bielectronic integral
2014-07-14 18:52:32 +02:00
`i_x1_pol_mult_a1 <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/ao_bi_integrals.irp.f#L839>`_
2014-05-13 13:57:58 +02:00
recursive function involved in the bielectronic integral
2014-07-14 18:52:32 +02:00
`i_x1_pol_mult_a2 <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/ao_bi_integrals.irp.f#L893>`_
2014-05-13 13:57:58 +02:00
recursive function involved in the bielectronic integral
2014-07-14 18:52:32 +02:00
`i_x1_pol_mult_recurs <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/ao_bi_integrals.irp.f#L753>`_
2014-05-13 13:57:58 +02:00
recursive function involved in the bielectronic integral
2014-07-14 18:52:32 +02:00
`i_x2_new <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/ao_bi_integrals.irp.f#L623>`_
2014-05-13 13:57:58 +02:00
recursive function involved in the bielectronic integral
2014-07-14 18:52:32 +02:00
`i_x2_pol_mult <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/ao_bi_integrals.irp.f#L955>`_
2014-05-13 13:57:58 +02:00
recursive function involved in the bielectronic integral
2014-07-14 18:52:32 +02:00
`integrale_new <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/ao_bi_integrals.irp.f#L555>`_
2014-05-13 13:57:58 +02:00
calculate the integral of the polynom ::
I_x1(a_x+b_x, c_x+d_x,p,q) * I_x1(a_y+b_y, c_y+d_y,p,q) * I_x1(a_z+b_z, c_z+d_z,p,q)
between ( 0 ; 1)
2014-07-14 18:52:32 +02:00
`n_pt_sup <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/ao_bi_integrals.irp.f#L642>`_
2014-05-13 13:57:58 +02:00
Returns the upper boundary of the degree of the polynom involved in the
bielctronic integral :
Ix(a_x,b_x,c_x,d_x) * Iy(a_y,b_y,c_y,d_y) * Iz(a_z,b_z,c_z,d_z)
2014-05-14 00:01:31 +02:00
`gauleg <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/gauss_legendre.irp.f#L20>`_
2014-05-13 13:57:58 +02:00
Gauss-Legendre
2014-05-14 00:01:31 +02:00
`gauleg_t2 <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/gauss_legendre.irp.f#L1>`_
2014-05-13 13:57:58 +02:00
t_w(i,1,k) = w(i)
t_w(i,2,k) = t(i)
2014-05-14 00:01:31 +02:00
`gauleg_w <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/gauss_legendre.irp.f#L2>`_
2014-05-13 13:57:58 +02:00
t_w(i,1,k) = w(i)
t_w(i,2,k) = t(i)
2014-05-14 00:01:31 +02:00
`ao_integrals_map <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/map_integrals.irp.f#L6>`_
2014-05-13 13:57:58 +02:00
AO integrals
2014-05-14 00:01:31 +02:00
`bielec_integrals_index <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/map_integrals.irp.f#L17>`_
Undocumented
2014-05-14 00:01:31 +02:00
`clear_ao_map <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/map_integrals.irp.f#L128>`_
2014-05-13 13:57:58 +02:00
Frees the memory of the AO map
2014-06-25 14:58:58 +02:00
`clear_mo_map <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/map_integrals.irp.f#L285>`_
2014-05-13 13:57:58 +02:00
Frees the memory of the MO map
2014-05-14 00:01:31 +02:00
`get_ao_bielec_integral <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/map_integrals.irp.f#L33>`_
2014-05-13 13:57:58 +02:00
Gets one AO bi-electronic integral from the AO map
2014-05-14 00:01:31 +02:00
`get_ao_bielec_integrals <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/map_integrals.irp.f#L51>`_
2014-05-13 13:57:58 +02:00
Gets multiple AO bi-electronic integral from the AO map .
All i are retrieved for j,k,l fixed.
2014-05-14 00:01:31 +02:00
`get_ao_bielec_integrals_non_zero <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/map_integrals.irp.f#L84>`_
2014-05-13 13:57:58 +02:00
Gets multiple AO bi-electronic integral from the AO map .
All non-zero i are retrieved for j,k,l fixed.
2014-05-14 00:01:31 +02:00
`get_ao_map_size <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/map_integrals.irp.f#L120>`_
2014-05-13 13:57:58 +02:00
Returns the number of elements in the AO map
2014-05-14 00:01:31 +02:00
`get_mo_bielec_integral <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/map_integrals.irp.f#L184>`_
2014-05-13 13:57:58 +02:00
Returns one integral <ij|kl> in the MO basis
2014-05-14 00:01:31 +02:00
`get_mo_bielec_integrals <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/map_integrals.irp.f#L213>`_
2014-05-13 13:57:58 +02:00
Returns multiple integrals <ij|kl> in the MO basis, all
i for j,k,l fixed.
2014-06-25 14:58:58 +02:00
`get_mo_bielec_integrals_existing_ik <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/map_integrals.irp.f#L235>`_
Returns multiple integrals <ij|kl> in the MO basis, all
i for j,k,l fixed.
`get_mo_map_size <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/map_integrals.irp.f#L277>`_
2014-05-13 13:57:58 +02:00
Return the number of elements in the MO map
2014-05-14 00:01:31 +02:00
`insert_into_ao_integrals_map <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/map_integrals.irp.f#L153>`_
2014-05-13 13:57:58 +02:00
Create new entry into AO map
2014-05-14 00:01:31 +02:00
`insert_into_mo_integrals_map <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/map_integrals.irp.f#L168>`_
2014-05-13 13:57:58 +02:00
Create new entry into MO map, or accumulate in an existing entry
2014-05-14 00:01:31 +02:00
`mo_bielec_integral <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/map_integrals.irp.f#L201>`_
2014-05-13 13:57:58 +02:00
Returns one integral <ij|kl> in the MO basis
2014-05-14 00:01:31 +02:00
`mo_integrals_map <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/map_integrals.irp.f#L142>`_
2014-05-13 13:57:58 +02:00
MO integrals
2014-05-14 00:01:31 +02:00
`add_integrals_to_map <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/mo_bi_integrals.irp.f#L40>`_
2014-05-13 13:57:58 +02:00
Adds integrals to tha MO map according to some bitmask
2014-07-14 18:52:32 +02:00
`mo_bielec_integral_jj <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/mo_bi_integrals.irp.f#L314>`_
mo_bielec_integral_jj(i,j) = J_ij
mo_bielec_integral_jj_exchange(i,j) = J_ij
mo_bielec_integral_jj_anti(i,j) = J_ij - K_ij
2014-05-13 13:57:58 +02:00
2014-07-14 18:52:32 +02:00
`mo_bielec_integral_jj_anti <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/mo_bi_integrals.irp.f#L316>`_
mo_bielec_integral_jj(i,j) = J_ij
mo_bielec_integral_jj_exchange(i,j) = J_ij
mo_bielec_integral_jj_anti(i,j) = J_ij - K_ij
2014-05-13 13:57:58 +02:00
2014-07-14 18:52:32 +02:00
`mo_bielec_integral_jj_exchange <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/mo_bi_integrals.irp.f#L315>`_
mo_bielec_integral_jj(i,j) = J_ij
mo_bielec_integral_jj_exchange(i,j) = J_ij
mo_bielec_integral_jj_anti(i,j) = J_ij - K_ij
2014-05-13 13:57:58 +02:00
2014-05-14 00:01:31 +02:00
`mo_bielec_integrals_in_map <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/mo_bi_integrals.irp.f#L21>`_
2014-05-13 13:57:58 +02:00
If True, the map of MO bielectronic integrals is provided
2014-05-14 00:01:31 +02:00
`mo_bielec_integrals_index <http://github.com/LCPQ/quantum_package/tree/master/src/BiInts/mo_bi_integrals.irp.f#L1>`_
2014-05-13 13:57:58 +02:00
Computes an unique index for i,j,k,l integrals