Task: Assume 46 and 78 core configurations to be f-in-valence (pyscf/pyscf#1891)

org: pyscf
repo: pyscf
pr_number: 1891
issue_title: Assume 46 and 78 core configurations to be f-in-valence
issue_body: This is a reasonable default and should fix  #1889. I added a test with `mf.get_init_guess()` to catch the error.

What should be the way to introduce another core configuration to the same element and number of electrons? 
An `f_in_core`  attribute should be added somewhere.
resolved_issues: number: 1889
title: Using ECP gives an error for Actinide elements
body: I am trying to do a HF/DFT calculation on a system involving a U atom and trying to invoke the effective core potential (ECP) 'crenbl' but that gives a numpy error:   LinAlgError('Last 2 dimensions of the array must be square'). The error goes away if ECP is turned off. The same thing happens when using other ECPs (eg. 'stuttgart') listed in the pyscf User Guide or on Basis set exchange for U and other Actinide elements. I am using pyscf v2.3.0 with numpy v1.26.0, scipy v1.11.3, h5py v3.9.0 and python v3.10.13 on Mac ARM architechture. The code is as follows:

```python
from pyscf import gto, scf, dft                                             

mol = gto.Mole()
mol.atom='''
U    -0.00000    0.00000    0.00000
H     -1.70048   -1.32378   -0.00000
H      1.70048   -1.32378   -0.00000
'''
mol.basis={'U':'crenbl','H':'6-31++g**'}                                      
mol.ecp={'U':'crenbl'}                                                        
mol.symmetry = True          
mol.build()

mf = dft.ROKS(mol)
mf.xc='b3lyp,b3lyp'
mf.max_cycle = 60
mf.conv_tol = 1e-9
mf.kernel()
``` 

Reference solution:
fix_patch: diff --git a/pyscf/gto/ecp.py b/pyscf/gto/ecp.py
index e686b88f26..d6a264318e 100644
--- a/pyscf/gto/ecp.py
+++ b/pyscf/gto/ecp.py
@@ -145,10 +145,12 @@ def core_configuration(nelec_core, atom_symbol=None):
     elements_4f = ELEMENTS[57:71]
     elements_5f = ELEMENTS[89:103]
     if atom_symbol in elements_4f:
-        for i in range(46, 60):
+        # TODO: fix La3+ and Ce4+ ECP46MWB f-in-core ECPs
+        for i in range(47, 59):
             conf_dic[i] = '4s3p2d1f'
     if atom_symbol in elements_5f:
-        for i in range(78, 92):
+        # TODO: fix Ac3+, Th4+, Pa5+ and U6+ ECP78MWB f-in-core ECPs
+        for i in range(79, 91):
             conf_dic[i] = '5s4p3d2f'
     if nelec_core not in conf_dic:
         raise RuntimeError('Core configuration for %d core electrons is not available.' % nelec_core)