psi4 / psi4

Open-Source Quantum Chemistry – an electronic structure package in C++ driven by Python
http://psicode.org
GNU Lesser General Public License v3.0
955 stars 441 forks source link

v1.7 PSIO_ERROR + forrtl: error + core dumped #2970

Open njzjz opened 1 year ago

njzjz commented 1 year ago

I use psi4 v1.7 to perform gradient calculations under the WB97M-D3BJ/DEF2-TZVPPD level. Most structures work fine, but I get the following error in some structures without a clear error message.

  ==> Integral Setup <==

PSIO_ERROR: unit = 97, errval = 12
forrtl: error (76): Abort trap signal
Image              PC                Routine            Line        Source
libgdma.so         00007F6F7C68B10F  for__signal_handl     Unknown  Unknown
libpthread-2.28.s  00007F6F9C906B20  Unknown               Unknown  Unknown
libc-2.28.so       00007F6F9BDDC37F  gsignal               Unknown  Unknown
libc-2.28.so       00007F6F9BDC6DB5  abort                 Unknown  Unknown
libstdc++.so.6.0.  00007F6F8C6B535A  Unknown               Unknown  Unknown
libstdc++.so.6.0.  00007F6F8C6B53C5  Unknown               Unknown  Unknown
libstdc++.so.6.0.  00007F6F8C6B5658  Unknown               Unknown  Unknown
core.cpython-310-  00007F6F88B5A129  Unknown               Unknown  Unknown
core.cpython-310-  00007F6F88B69438  _ZN3psi4PSIO2rwEm     Unknown  Unknown
core.cpython-310-  00007F6F88B6B963  _ZN3psi4PSIO5writ     Unknown  Unknown
core.cpython-310-  00007F6F88B566E9  Unknown               Unknown  Unknown
libstdc++.so.6.0.  00007F6F8C6DFBF4  Unknown               Unknown  Unknown
libpthread-2.28.s  00007F6F9C8FC14A  Unknown               Unknown  Unknown
libc-2.28.so       00007F6F9BEA1DC3  clone                 Unknown  Unknown
/var/spool/slurm/slurmd/job83930/slurm_script: line 27: 2245656 Aborted                 (core dumped)

I used conda to install psi4 in a new environment.

psi4                      1.7+6ce35a5     py310h769327c_1    psi4

Below is the full log, which includes the geometry:

``` Memory set to 7.451 GiB by Python driver. Threads set to 4 by Python driver. Scratch directory: /scratch/njzjz/83930/ gradient() will perform analytic gradient computation. *** tstart() called on aerosol-020 *** at Fri May 26 20:10:38 2023 => Loading Basis Set <= Name: DEF2-TZVPPD Role: ORBITAL Keyword: BASIS atoms 1-33 entry C line 144 file /aerosol/users/njzjz/anaconda3/envs/p4env/share/psi4/basis/def2-tzvppd.gbs atoms 34-89 entry H line 14 file /aerosol/users/njzjz/anaconda3/envs/p4env/share/psi4/basis/def2-tzvppd.gbs atoms 90-101 entry N line 181 file /aerosol/users/njzjz/anaconda3/envs/p4env/share/psi4/basis/def2-tzvppd.gbs atoms 102-106 entry O line 218 file /aerosol/users/njzjz/anaconda3/envs/p4env/share/psi4/basis/def2-tzvppd.gbs => WB97M-D3BJ: Empirical Dispersion <= Grimme's -D3 (BJ-damping) Dispersion Correction Grimme S.; Ehrlich S.; Goerigk L. (2011), J. Comput. Chem., 32: 1456 Parametrisation from: A. Najib, L. Goerigk, J. Comput. Theory Chem.,14, 5725, 2018 s6 = 1.000000 s8 = 0.390800 a1 = 0.566000 a2 = 3.128000 --------------------------------------------------------- SCF by Justin Turney, Rob Parrish, Andy Simmonett and Daniel G. A. Smith RKS Reference 4 Threads, 7629 MiB Core --------------------------------------------------------- ==> Geometry <== Molecular point group: c1 Full point group: C1 Geometry (in Angstrom), charge = 0, multiplicity = 1: Center X Y Z Mass ------------ ----------------- ----------------- ----------------- ----------------- C -3.986536731580 -1.858264338730 0.058762714511 12.000000000000 C 0.798074669820 0.990515156370 -3.875762257089 12.000000000000 C 0.017832511920 -3.811459520430 -1.287345979189 12.000000000000 C 0.173626034620 2.867679137670 -2.668652350389 12.000000000000 C 1.516908962420 -3.483999515230 -1.216553356289 12.000000000000 C 5.992839611420 -1.026305827230 -0.727874868189 12.000000000000 C 3.017409895220 -3.046872095730 0.725792373811 12.000000000000 C -1.125684508680 2.202468723670 -2.198266685089 12.000000000000 C 3.752417747320 -2.017670362830 -0.145942986389 12.000000000000 C 2.625927257520 -0.162195960130 1.108527837311 12.000000000000 C 2.049668656220 1.018254628770 0.683856523611 12.000000000000 C 1.406278426320 2.114882744470 1.477238232711 12.000000000000 C -0.252334991380 -4.016660671530 1.099960464211 12.000000000000 C 2.325525208620 3.352293930170 1.620006393311 12.000000000000 C 3.704698435020 2.984894890570 2.165178250911 12.000000000000 C 1.173584381620 -4.581121691030 0.954011933511 12.000000000000 C 2.431556537520 4.093718609670 0.285153063211 12.000000000000 C 3.257587384320 5.371619339270 0.370658903411 12.000000000000 C 6.703150146120 1.001690449270 0.649806092711 12.000000000000 C 5.203542340120 -1.826786734930 0.310957935811 12.000000000000 C -0.828203820280 -1.948650900430 -0.035813584789 12.000000000000 C -0.467415497580 0.236497941770 -3.448468542789 12.000000000000 C 7.161038878020 -0.248602982430 -0.120716922289 12.000000000000 C -6.534117912180 -0.041136960130 1.763394298811 12.000000000000 C -7.040379234280 -0.106198240430 -0.536480292689 12.000000000000 C -1.008267407080 0.037590715470 -1.074714631289 12.000000000000 C -6.042461309680 0.895000830170 2.873852988911 12.000000000000 C -3.973926668980 1.640858796870 2.005172149011 12.000000000000 C -6.620264830580 0.512521693870 -1.874786852489 12.000000000000 C -2.487781161280 1.699687470070 2.358288612011 12.000000000000 C -1.496924334280 -0.143197204530 1.110276498111 12.000000000000 C -4.111966726280 -1.021753630730 -0.788303496889 12.000000000000 C -4.297279254080 -0.041549042130 -1.847183022189 12.000000000000 H 5.655603749820 -2.807208432730 0.464992325511 1.007825032230 H 5.201063217220 -1.318310892230 1.275493607811 1.007825032230 H -6.435591512680 1.903316193570 2.678134795611 1.007825032230 H -7.598598689480 -0.262057246930 1.937253906411 1.007825032230 H 6.368254448120 -1.707196721930 -1.494856337789 1.007825032230 H 5.320665048720 -0.323924208830 -1.229193601189 1.007825032230 H 7.722840764920 -0.898092332830 0.552778297411 1.007825032230 H -4.382652052880 2.659093431870 1.926273744011 1.007825032230 H 5.983929857420 0.718965391170 1.423232172611 1.007825032230 H 7.567087729120 1.447140715970 1.150601776411 1.007825032230 H 5.271581718820 1.691129301270 -0.641074156989 1.007825032230 H 6.740175270120 2.378355068770 -0.854873144789 1.007825032230 H -0.265043214380 -0.824964875630 -3.316711869289 1.007825032230 H -0.958234964980 -4.838498815530 1.265296774511 1.007825032230 H 1.655270784320 -4.672476721130 1.928170705211 1.007825032230 H -0.321256210380 -3.315793018530 1.933385934611 1.007825032230 H -1.973731587680 2.348813071670 1.649690443911 1.007825032230 H 1.157510957220 -5.564392103030 0.473682293511 1.007825032230 H 1.627549171620 -2.437462251930 -1.506128084489 1.007825032230 H -2.364277417280 2.085763061570 3.371579814811 1.007825032230 H -2.150448301980 -0.265490338430 2.991977227411 1.007825032230 H 2.088475342920 -4.116827186330 -1.907643331389 1.007825032230 H -1.937550982380 2.452578387770 -2.892926055089 1.007825032230 H 0.966209705420 2.655137431770 -1.938755781189 1.007825032230 H 1.602225912420 0.784184209570 -3.157164741189 1.007825032230 H 0.040231465620 3.944776444770 -2.781698926989 1.007825032230 H -1.246704760280 0.378097808570 -4.208468393589 1.007825032230 H 1.105761063120 0.687314982270 -4.878224780789 1.007825032230 H -4.102474244780 1.133612013170 1.043202319711 1.007825032230 H -1.386756860080 2.517921213070 -1.189495861789 1.007825032230 H 4.301329506720 5.146585664670 0.569491375311 1.007825032230 H 3.745380311120 -2.302604301930 -1.201110224389 1.007825032230 H -3.543648590980 -0.226374997630 -2.621427033889 1.007825032230 H -4.147722297280 0.968174695270 -1.435427673389 1.007825032230 H -6.385345850780 1.572080401770 -1.693640161989 1.007825032230 H -8.129348911680 -0.008353047530 -0.402227810089 1.007825032230 H -6.783560009480 -1.173926656530 -0.530477449489 1.007825032230 H -5.969000464880 -0.983732799130 1.794283612611 1.007825032230 H -6.429532541280 0.534577281670 3.830475615811 1.007825032230 H -0.164226375180 -4.887055000930 -1.334640467389 1.007825032230 H -0.424643353080 -3.328691049330 -2.157438422789 1.007825032230 H 7.835500291220 0.065938515470 -0.921769781789 1.007825032230 H 0.280221214920 0.956683993470 2.728645658111 1.007825032230 H 1.747024140320 1.359908648670 3.346038551411 1.007825032230 H -3.829421013780 -2.571106092430 0.812789978711 1.007825032230 H 4.203719336220 3.869463925070 2.550812242411 1.007825032230 H 2.764653104420 -0.607249357830 2.070992087411 1.007825032230 H 1.822963565420 3.998962029270 2.347115302211 1.007825032230 H 2.873313233220 3.428326310170 -0.458437339889 1.007825032230 H 1.423296421920 4.345057722370 -0.053103985889 1.007825032230 H 4.329642385220 2.561431424770 1.383922458011 1.007825032230 H 2.881730458020 6.017326375870 1.161775109511 1.007825032230 H 3.625277287820 2.271986037870 2.983011256011 1.007825032230 H -7.443463685380 0.441052298070 -2.589101689389 1.007825032230 H 3.201106403020 5.914044369070 -0.570468338989 1.007825032230 H 0.508864437320 2.424575945370 0.926543063511 1.007825032230 N 6.089376920620 2.042268910870 -0.153045002889 14.003074004430 N -1.862112502780 0.403485545970 2.293716133311 14.003074004430 N -1.413868267180 0.644418478470 0.045530940511 14.003074004430 N -0.709496196580 -1.263243640330 -1.167615241789 14.003074004430 N -0.614250361980 -3.293788727030 -0.100694759089 14.003074004430 N 1.978424723220 -3.681420601330 0.149257795111 14.003074004430 N -0.943697896280 0.770718955070 -2.195103718789 14.003074004430 N 3.030604932020 -0.774545035830 -0.027468195689 14.003074004430 N 2.139390467920 1.040663076670 -0.664008978489 14.003074004430 N 2.724104508920 -0.024887744330 -1.077948484789 14.003074004430 N 0.964366717620 1.701693516370 2.799565885411 14.003074004430 N -1.189761996980 -1.443402396430 1.139850386111 14.003074004430 O -4.650460781380 0.935132853670 3.030086427411 15.994914619570 O 0.553972156620 2.379030254270 -3.938430599989 15.994914619570 O -6.376361129480 0.574029773370 0.505587091411 15.994914619570 O -5.547778804180 -0.131717416630 -2.508918002589 15.994914619570 O 3.352911070220 -3.242204449630 1.878767812611 15.994914619570 Running in c1 symmetry. Rotational constants: A = 0.00301 B = 0.00144 C = 0.00130 [cm^-1] Rotational constants: A = 90.10550 B = 43.03146 C = 38.85152 [MHz] Nuclear repulsion = 7705.917246066113876 Charge = 0 Multiplicity = 1 Electrons = 378 Nalpha = 189 Nbeta = 189 ==> Algorithm <== SCF Algorithm Type is DF. DIIS enabled. MOM disabled. Fractional occupation disabled. Guess Type is SAD. Energy threshold = 1.00e-08 Density threshold = 1.00e-08 Integral threshold = 1.00e-12 ==> Primary Basis <== Basis Set: DEF2-TZVPPD Blend: DEF2-TZVPPD Number of shells: 1047 Number of basis functions: 2817 Number of Cartesian functions: 3173 Spherical Harmonics?: true Max angular momentum: 3 ==> DFT Potential <== => LibXC <= Version 5.2.3 S. Lehtola, C. Steigemann, M. J. Oliveira, and M. A. Marques, SoftwareX 7, 1 (2018) (10.1016/j.softx.2017.11.002) => Composite Functional: WB97M-D3BJ <= wB97M-V with D3(BJ) instead of VV10 dispersion A. Najib, L. Goerigk, J. Comput. Theory Chem.,14, 5725, 2018 N. Mardirossian, M. Head-Gordon, J. Chem. Phys. 144, 214110, 2016 Deriv = 1 GGA = TRUE Meta = TRUE Exchange Hybrid = TRUE MP2 Hybrid = FALSE => Exchange-Correlation Functionals <= 1.0000 wB97M-V exchange-correlation functional => Exact (HF) Exchange <= 0.8500 HF,LR [omega = 0.3000] 0.1500 HF => LibXC Density Thresholds <== XC_HYB_MGGA_XC_WB97M_V: 1.00E-13 => Molecular Quadrature <= Radial Scheme = TREUTLER Pruning Scheme = NONE Nuclear Scheme = TREUTLER Blocking Scheme = OCTREE BS radius alpha = 1 Pruning alpha = 1 Radial Points = 75 Spherical Points = 302 Total Points = 2115829 Total Blocks = 14225 Max Points = 256 Max Functions = 1539 Weights Tolerance = 1.00E-15 => Loading Basis Set <= Name: (DEF2-TZVPPD AUX) Role: JKFIT Keyword: DF_BASIS_SCF atoms 1-33 entry C line 198 file /aerosol/users/njzjz/anaconda3/envs/p4env/share/psi4/basis/def2-universal-jkfit.gbs atoms 34-89 entry H line 18 file /aerosol/users/njzjz/anaconda3/envs/p4env/share/psi4/basis/def2-universal-jkfit.gbs atoms 90-101 entry N line 258 file /aerosol/users/njzjz/anaconda3/envs/p4env/share/psi4/basis/def2-universal-jkfit.gbs atoms 102-106 entry O line 318 file /aerosol/users/njzjz/anaconda3/envs/p4env/share/psi4/basis/def2-universal-jkfit.gbs ==> Integral Setup <== PSIO_ERROR: unit = 97, errval = 12 forrtl: error (76): Abort trap signal Image PC Routine Line Source libgdma.so 00007F6F7C68B10F for__signal_handl Unknown Unknown libpthread-2.28.s 00007F6F9C906B20 Unknown Unknown Unknown libc-2.28.so 00007F6F9BDDC37F gsignal Unknown Unknown libc-2.28.so 00007F6F9BDC6DB5 abort Unknown Unknown libstdc++.so.6.0. 00007F6F8C6B535A Unknown Unknown Unknown libstdc++.so.6.0. 00007F6F8C6B53C5 Unknown Unknown Unknown libstdc++.so.6.0. 00007F6F8C6B5658 Unknown Unknown Unknown core.cpython-310- 00007F6F88B5A129 Unknown Unknown Unknown core.cpython-310- 00007F6F88B69438 _ZN3psi4PSIO2rwEm Unknown Unknown core.cpython-310- 00007F6F88B6B963 _ZN3psi4PSIO5writ Unknown Unknown core.cpython-310- 00007F6F88B566E9 Unknown Unknown Unknown libstdc++.so.6.0. 00007F6F8C6DFBF4 Unknown Unknown Unknown libpthread-2.28.s 00007F6F9C8FC14A Unknown Unknown Unknown libc-2.28.so 00007F6F9BEA1DC3 clone Unknown Unknown /var/spool/slurm/slurmd/job83930/slurm_script: line 27: 2245656 Aborted (core dumped) ```

CPU information:

``` processor : 15 vendor_id : GenuineIntel cpu family : 6 model : 45 model name : Intel(R) Xeon(R) CPU E5-2670 0 @ 2.60GHz stepping : 6 microcode : 0x621 cpu MHz : 1832.125 cache size : 20480 KB physical id : 1 siblings : 8 core id : 7 cpu cores : 8 apicid : 46 initial apicid : 46 fpu : yes fpu_exception : yes cpuid level : 13 wp : yes flags : fpu vme de pse tsc msr pae mce cx8 apic sep mtrr pge mca cmov pat pse36 clflush dts acpi mmx fxsr sse sse2 ss ht tm pbe syscall nx pdpe1gb rdtscp lm constant_tsc arch_perfmon pebs bts rep_good nopl xtopology nonstop_tsc cpuid aperfmperf pni pclmulqdq dtes64 monitor ds_cpl vmx smx est tm2 ssse3 cx16 xtpr pdcm pcid dca sse4_1 sse4_2 x2apic popcnt tsc_deadline_timer aes xsave avx lahf_lm pti ssbd ibrs ibpb stibp tpr_shadow vnmi flexpriority ept vpid xsaveopt dtherm arat pln pts md_clear flush_l1d bugs : cpu_meltdown spectre_v1 spectre_v2 spec_store_bypass l1tf mds swapgs itlb_multihit bogomips : 5194.34 clflush size : 64 cache_alignment : 64 address sizes : 46 bits physical, 48 bits virtual power management: ```
TiborGY commented 1 year ago

PSIO error 12 means write error: https://github.com/psi4/psi4/blob/58efb602ee9fdf3d774a3a6f073dbf47f2ac249b/psi4/src/psi4/libpsio/config.h#L55

Unfortunately your environment does not seem to give you the rest of the error message, so we dont actually know why a write has failed. Are you sure you have enough free disk space for the temporary files?

njzjz commented 1 year ago

The disk space for the temporary file is about 245 GB. I don't know exactly how much space it requires.

Could this kind of error raise a Python exception? I use try...except to mark failed points, and core dumped can not be caught within the Python script.

loriab commented 1 year ago

The input file would be needed to see if it can be reproduced.

njzjz commented 1 year ago

Thanks. Below is the code to call psi4. I have omitted the code for reading data and processing and storing data.

import psi4
import numpy as np

# Actually, data is read from the file
data = {
    "atom_names": ["C", "H", "N", "O"],
    "atom_types": np.array([0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3]),
    "coords": np.array([[[-3.9865585648, -1.8582760752, 0.0587669877], [0.7980528366, 0.9905034199, -3.8757579839], [0.0178106787, -3.8114712569, -1.287341706], [0.1736042014, 2.8676674012, -2.6686480772], [1.5168871292, -3.4840112517, -1.2165490831], [5.9928177782, -1.0263175637, -0.727870595], [3.017388062, -3.0468838322, 0.725796647], [-1.1257063419, 2.2024569872, -2.1982624119], [3.7523959141, -2.0176820993, -0.1459387132], [2.6259054243, -0.1622076966, 1.1085321105], [2.049646823, 1.0182428923, 0.6838607968], [1.4062565931, 2.114871008, 1.4772425059], [-0.2523568246, -4.016672408, 1.0999647374], [2.3255033754, 3.3522821937, 1.6200106665], [3.7046766018, 2.9848831541, 2.1651825241], [1.1735625484, -4.5811334275, 0.9540162067], [2.4315347043, 4.0937068732, 0.2851573364], [3.2575655511, 5.3716076028, 0.3706631766], [6.7031283129, 1.0016787128, 0.6498103659], [5.2035205069, -1.8267984714, 0.310962209], [-0.8282256535, -1.9486626369, -0.0358093116], [-0.4674373308, 0.2364862053, -3.4484642696], [7.1610170448, -0.2486147189, -0.1207126491], [-6.5341397454, -0.0411486966, 1.763398572], [-7.0404010675, -0.1062099769, -0.5364760195], [-1.0082892403, 0.037578979, -1.0747103581], [-6.0424831429, 0.8949890937, 2.8738572621], [-3.9739485022, 1.6408470604, 2.0051764222], [-6.6202866638, 0.5125099574, -1.8747825793], [-2.4878029945, 1.6996757336, 2.3582928852], [-1.4969461675, -0.143208941, 1.1102807713], [-4.1119885595, -1.0217653672, -0.7882992237], [-4.2973010873, -0.0415607786, -1.847178749], [5.6555819166, -2.8072201692, 0.4649965987], [5.201041384, -1.3183226287, 1.275497881], [-6.4356133459, 1.9033044571, 2.6781390688], [-7.5986205227, -0.2620689834, 1.9372581796], [6.3682326149, -1.7072084584, -1.4948520646], [5.3206432155, -0.3239359453, -1.229189328], [7.7228189317, -0.8981040693, 0.5527825706], [-4.3826738861, 2.6590816954, 1.9262780172], [5.9839080242, 0.7189536547, 1.4232364458], [7.5670658959, 1.4471289795, 1.1506060496], [5.2715598856, 1.6911175648, -0.6410698838], [6.7401534369, 2.3783433323, -0.8548688716], [-0.2650650476, -0.8249766121, -3.3167075961], [-0.9582567982, -4.838510552, 1.2653010477], [1.6552489511, -4.6724884576, 1.9281749784], [-0.3212780436, -3.315804755, 1.9333902078], [-1.9737534209, 2.3488013352, 1.6496947171], [1.157489124, -5.5644038395, 0.4736865667], [1.6275273384, -2.4374739884, -1.5061238113], [-2.3642992505, 2.0857513251, 3.371584088], [-2.1504701352, -0.2655020749, 2.9919815006], [2.0884535097, -4.1168389228, -1.9076390582], [-1.9375728156, 2.4525666513, -2.8929217819], [0.9661878722, 2.6551256953, -1.938751508], [1.6022040792, 0.7841724731, -3.157160468], [0.0402096324, 3.9447647083, -2.7816946538], [-1.2467265935, 0.3780860721, -4.2084641204], [1.1057392299, 0.6873032458, -4.8782205076], [-4.102496078, 1.1336002767, 1.0432065929], [-1.3867786933, 2.5179094766, -1.1894915886], [4.3013076735, 5.1465739282, 0.5694956485], [3.7453584779, -2.3026160384, -1.2011059512], [-3.5436704242, -0.2263867341, -2.6214227607], [-4.1477441305, 0.9681629588, -1.4354234002], [-6.385367684, 1.5720686653, -1.6936358888], [-8.1293707449, -0.008364784, -0.4022235369], [-6.7835818427, -1.173938393, -0.5304731763], [-5.9690222981, -0.9837445356, 1.7942878858], [-6.4295543745, 0.5345655452, 3.830479889], [-0.1642482084, -4.8870667374, -1.3346361942], [-0.4246651863, -3.3287027858, -2.1574341496], [7.835478458, 0.065926779, -0.9217655086], [0.2801993817, 0.956672257, 2.7286499313], [1.7470023071, 1.3598969122, 3.3460428246], [-3.829442847, -2.5711178289, 0.8127942519], [4.203697503, 3.8694521886, 2.5508165156], [2.7646312712, -0.6072610943, 2.0709963606], [1.8229417322, 3.9989502928, 2.3471195754], [2.8732914, 3.4283145737, -0.4584330667], [1.4232745887, 4.3450459859, -0.0530997127], [4.329620552, 2.5614196883, 1.3839267312], [2.8817086248, 6.0173146394, 1.1617793827], [3.6252554546, 2.2719743014, 2.9830155292], [-7.4434855186, 0.4410405616, -2.5890974162], [3.2010845698, 5.9140326326, -0.5704640658], [0.5088426041, 2.4245642089, 0.9265473367], [6.0893550874, 2.0422571744, -0.1530407297], [-1.862134336, 0.4034738095, 2.2937204065], [-1.4138901004, 0.644406742, 0.0455352137], [-0.7095180298, -1.2632553768, -1.1676109686], [-0.6142721952, -3.2938004635, -0.1006904859], [1.97840289, -3.6814323378, 0.1492620683], [-0.9437197295, 0.7707072186, -2.1950994456], [3.0305830988, -0.7745567723, -0.0274639225], [2.1393686347, 1.0406513402, -0.6640047053], [2.7240826757, -0.0248994808, -1.0779442116], [0.9643448844, 1.7016817799, 2.7995701586], [-1.1897838302, -1.4434141329, 1.1398546593], [-4.6504826146, 0.9351211172, 3.0300907006], [0.5539503234, 2.3790185178, -3.9384263268], [-6.3763829627, 0.5740180369, 0.5055913646], [-5.5478006374, -0.1317291531, -2.5089137294], [3.352889237, -3.2422161861, 1.8787720858]]])
}

psi4.set_memory("8 GB")
psi4.set_num_threads(4)
types = np.array(data["atom_names"])[data["atom_types"]]
buff = ["0 1"]
for tt, cc in zip(types, data["coords"][0]):
    buff.append(" ".join([tt] + [str(x) for x in cc]))
psi4.geometry("\n".join(buff))

G, wfn = psi4.gradient("WB97M-D3BJ/def2-TZVPPD", return_wfn=True)
energy = wfn.energy()
wfn.gradient().print_out()
gradient = np.array(G)

# post process energy and gradient
susilehtola commented 8 months ago

I bet this is because of the calls to exit() in libpsio, combined with the lack of flushing of the error messages to the output. I have replaced the calls to exit() in #3118, which should solve this issue.