Ejemplos de HGTO en Python

Ejemplo n.º 1

Archivo: hydrogen-p.py Proyecto: lisyarus/chembook

def make_basis(params):
    basis = []

    for sigma in params:
        basis.append(hgto.HGTO(sigma, [1, 0, 0]))

    return basis

Ejemplo n.º 2

Archivo: helium-hydride-r.py Proyecto: lisyarus/chembook

def make_basis(R, params):
    nuclei = [
        hgto.Nucleus(2.0, [-R / 2, 0.0, 0.0]),
        hgto.Nucleus(1.0, [R / 2, 0.0, 0.0]),
    ]

    N = len(params) // 2

    sigma_He = params[0:N]
    sigma_H = params[N:]

    basis = []

    for sigma in sigma_He:
        basis.append(hgto.HGTO(sigma, origin=nuclei[0].position))

    for sigma in sigma_H:
        basis.append(hgto.HGTO(sigma, origin=nuclei[1].position))

    return basis, nuclei

Ejemplo n.º 3

Archivo: dihydrogen-r.py Proyecto: lisyarus/chembook

def make_basis(R, params):
    nuclei = [
        hgto.Nucleus(1.0, [-R / 2, 0.0, 0.0]),
        hgto.Nucleus(1.0, [R / 2, 0.0, 0.0]),
    ]

    basis = []

    for sigma in params:
        for n in nuclei:
            basis.append(hgto.HGTO(sigma, origin=n.position))

    return basis, nuclei

Ejemplo n.º 4

Archivo: trihydrogen.py Proyecto: lisyarus/chembook

def make_basis(params):
    A = 1.7
    R = A / numpy.sqrt(3.0)

    nuclei = [
        hgto.Nucleus(1.0, [R, 0.0, 0.0]),
        hgto.Nucleus(1.0, [-R / 2, A / 2, 0.0]),
        hgto.Nucleus(1.0, [-R / 2, -A / 2, 0.0]),
    ]

    basis = []

    for sigma in params:
        for n in nuclei:
            basis.append(hgto.HGTO(sigma, origin=n.position))

    return basis, nuclei