Difference between revisions of "Oparse.m"
(Created page with "A parser for ORCA calculation logs. ==Syntax== props=oparse(filename) ==Description== The function extracts all potentially useful information from the ORCA log. Exam...") |
(sync with Spinach main 3975f139: hfc.iso NaN when not printed) |
||
| (17 intermediate revisions by 3 users not shown) | |||
| Line 1: | Line 1: | ||
| − | A parser for ORCA | + | {{DISPLAYTITLE:oparse.m}} __NOTOC__ |
| + | A parser for ORCA text output logs, versions 2.6 to 6.1. Reads the geometry and every magnetic parameter that ORCA prints in the main output file. | ||
==Syntax== | ==Syntax== | ||
| − | props=oparse( | + | props=oparse(file_name) |
| − | == | + | ==Parameters== |
| − | |||
| − | + | file_name - a character string with the file path | |
| − | == | + | ==Outputs== |
| − | + | props.filename - log file name | |
| − | + | props.orca_version - ORCA version string | |
| − | + | props.symbols - atomic symbols, 1 x natoms cell | |
| − | + | props.atomic_numbers - atomic numbers, 1 x natoms | |
| − | + | props.std_geom - atomic coordinates, natoms x 3, Angstrom | |
| − | + | props.natoms - number of atoms | |
| − | + | props.charge - total charge | |
| + | props.multiplicity - spin multiplicity | ||
| + | props.energy - final single point energy, Hartree | ||
| + | props.dip_moment - electric dipole moment, a.u. | ||
| + | props.mulliken_chg - Mulliken atomic charges, natoms x 1 | ||
| + | props.mulliken_spin - Mulliken spin populations, natoms x 1 | ||
| + | props.g_tensor.raw - g-matrix as printed by ORCA | ||
| + | props.g_tensor.matrix - symmetrised g-matrix | ||
| + | props.g_tensor.eigvals - eigenvalues of the symmetrised g-matrix | ||
| + | props.g_tensor.eigvecs - eigenvectors of the symmetrised g-matrix | ||
| + | props.zfs.matrix - zero-field splitting tensor, cm^-1 | ||
| + | props.zfs.eigvals - ZFS tensor eigenvalues, cm^-1 | ||
| + | props.zfs.eigvecs - ZFS tensor eigenvectors | ||
| + | props.hfc.full.matrix - hyperfine tensors, Gauss, natoms cell | ||
| + | props.hfc.full.eigvals - hyperfine eigenvalues, Gauss, natoms cell | ||
| + | props.hfc.full.eigvecs - hyperfine eigenvectors, natoms cell | ||
| + | props.hfc.iso - isotropic hyperfine couplings, Gauss, NaN if not printed | ||
| + | props.efg - EFG tensors, a.u.^-3, natoms cell | ||
| + | props.nqi - quadrupolar tensors, Hz, natoms cell | ||
| + | props.isotopes - isotopes used by ORCA, natoms cell | ||
| + | props.cst - shielding tensors, ppm, natoms cell | ||
| + | props.j_couplings - isotropic J-couplings, Hz, natoms x natoms | ||
| + | props.chi_temps - susceptibility temperatures, K | ||
| + | props.chi_tensors - molar magnetic susceptibility tensors, | ||
| + | cm^3*K/mol, one cell per temperature | ||
| + | |||
| + | Only the fields that ORCA has actually printed are returned; the caller should test for their presence with isfield. | ||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
==Examples== | ==Examples== | ||
None at the moment. | None at the moment. | ||
| Line 43: | Line 52: | ||
# Chemical shielding is ''not'' the same as chemical shift. | # Chemical shielding is ''not'' the same as chemical shift. | ||
# This function parses ORCA logs. Use [[g2spinach.m]] to convert that information into ''Spinach'' input structures. | # This function parses ORCA logs. Use [[g2spinach.m]] to convert that information into ''Spinach'' input structures. | ||
| − | # | + | # ORCA prints magnetic parameters only for the nuclei that were requested in the input, and labels each of them with the zero based index of the atom in the Cartesian coordinate table. All per-atom outputs above are therefore indexed by the position of the atom in props.std_geom, and are left empty for atoms whose parameters were not printed. |
| + | # When a log contains multiple geometries or multiple property sections, for example a geometry optimisation or a relaxed surface scan, the last one printed is returned. | ||
==See also== | ==See also== | ||
| − | [[gparse.m]], [[g2spinach.m]] | + | [[gparse.m]], [[g2spinach.m]], [[c2spinach.m]], [[cyprinol.m]], [[fatty_acid.m]], [[gissmo2spinach.m]], [[karplus_fit.m]], [[killcross.m]], [[killdiag.m]], [[merge_inp.m]], [[methyl_group.m]], [[nuclacid.m]], [[ocparse.m]], [[parsexml.m]], [[protein.m]], [[read_bmrb.m]], [[read_pdb_nuc.m]], [[read_pdb_pro.m]], [[v2spinach.m]], [[weblab2nqi.m]], [[x2spinach.m]], [[Import,_export,_and_visualisation]] |
| − | |||
| − | '' | + | ''Version 1.9, authors: [[Ilya Kuprov]], [[Elizaveta Suturina]]'' |
Latest revision as of 10:58, 18 September 2026
A parser for ORCA text output logs, versions 2.6 to 6.1. Reads the geometry and every magnetic parameter that ORCA prints in the main output file.
Syntax
props=oparse(file_name)
Parameters
file_name - a character string with the file path
Outputs
props.filename - log file name
props.orca_version - ORCA version string
props.symbols - atomic symbols, 1 x natoms cell
props.atomic_numbers - atomic numbers, 1 x natoms
props.std_geom - atomic coordinates, natoms x 3, Angstrom
props.natoms - number of atoms
props.charge - total charge
props.multiplicity - spin multiplicity
props.energy - final single point energy, Hartree
props.dip_moment - electric dipole moment, a.u.
props.mulliken_chg - Mulliken atomic charges, natoms x 1
props.mulliken_spin - Mulliken spin populations, natoms x 1
props.g_tensor.raw - g-matrix as printed by ORCA
props.g_tensor.matrix - symmetrised g-matrix
props.g_tensor.eigvals - eigenvalues of the symmetrised g-matrix
props.g_tensor.eigvecs - eigenvectors of the symmetrised g-matrix
props.zfs.matrix - zero-field splitting tensor, cm^-1
props.zfs.eigvals - ZFS tensor eigenvalues, cm^-1
props.zfs.eigvecs - ZFS tensor eigenvectors
props.hfc.full.matrix - hyperfine tensors, Gauss, natoms cell
props.hfc.full.eigvals - hyperfine eigenvalues, Gauss, natoms cell
props.hfc.full.eigvecs - hyperfine eigenvectors, natoms cell
props.hfc.iso - isotropic hyperfine couplings, Gauss, NaN if not printed
props.efg - EFG tensors, a.u.^-3, natoms cell
props.nqi - quadrupolar tensors, Hz, natoms cell
props.isotopes - isotopes used by ORCA, natoms cell
props.cst - shielding tensors, ppm, natoms cell
props.j_couplings - isotropic J-couplings, Hz, natoms x natoms
props.chi_temps - susceptibility temperatures, K
props.chi_tensors - molar magnetic susceptibility tensors,
cm^3*K/mol, one cell per temperature
Only the fields that ORCA has actually printed are returned; the caller should test for their presence with isfield.
Examples
None at the moment.
Notes
- Chemical shielding is not the same as chemical shift.
- This function parses ORCA logs. Use g2spinach.m to convert that information into Spinach input structures.
- ORCA prints magnetic parameters only for the nuclei that were requested in the input, and labels each of them with the zero based index of the atom in the Cartesian coordinate table. All per-atom outputs above are therefore indexed by the position of the atom in props.std_geom, and are left empty for atoms whose parameters were not printed.
- When a log contains multiple geometries or multiple property sections, for example a geometry optimisation or a relaxed surface scan, the last one printed is returned.
See also
gparse.m, g2spinach.m, c2spinach.m, cyprinol.m, fatty_acid.m, gissmo2spinach.m, karplus_fit.m, killcross.m, killdiag.m, merge_inp.m, methyl_group.m, nuclacid.m, ocparse.m, parsexml.m, protein.m, read_bmrb.m, read_pdb_nuc.m, read_pdb_pro.m, v2spinach.m, weblab2nqi.m, x2spinach.m, Import,_export,_and_visualisation
Version 1.9, authors: Ilya Kuprov, Elizaveta Suturina