# smint

**Repository Path**: luyi07/smint

## Basic Information

- **Project Name**: smint
- **Description**: shell model interactions  iso/upn/xpn --(Pandya)-->  upn/xpn
- **Primary Language**: C++
- **License**: MIT
- **Default Branch**: master
- **Homepage**: None
- **GVP Project**: No

## Statistics

- **Stars**: 0
- **Forks**: 0
- **Created**: 2022-01-11
- **Last Updated**: 2025-09-17

## Categories & Tags

**Categories**: Uncategorized

**Tags**: None

## README

# smint

## 1. compile

make

## 2. example

iso2upn: under smint/example/

```
cp ../smint.x .

bash test.sh
```

jpq2xpn: under smint/example/

```
cp ../smint.x .

bash test-jpq2xpn.sh sd.sp sd.jpq 2use.xpn
```

I have tested jpq2xpn with usdb interactions (after iso2jpq). It is correct. But such jpq files do not have QpQp or QnQn terms, only have QpQn terms, so QpQp/QnQn have not been checked yet, strictly speaking.

## 3. new format "oneupn"

We can consider protons and neutrons as one kind of nucleons in different single particle orbits. Then we use upn format convention, and that is "oneupn".
We take sd shell as an example.

### Edit a new single particle file

Firstly we edit a new sp file as example/sd.one.upn.sp, to define the space, such that there are 12 single particle orbits, i.e., 6 "proton orbits" and 6 "neutron orbits". Let 3 "proton orbits" correspond to the real neutron orbits. To differentiate them, we add 1 to the quantum number n of 3 neutron orbits.

```
1 1 2 3
2 1 2 5
3 2 0 1
4 2 2 3
5 2 2 5
6 3 0 1
```

In total, example/sd.one.upn.sp is

```
! sd shell
pn
16 8
12
2 6 6
1 1 2 3
2 1 2 5
3 2 0 1
4 2 2 3
5 2 2 5
6 3 0 1
1 1 2 3
2 1 2 5
3 2 0 1
4 2 2 3
5 2 2 5
6 3 0 1
```

### Run the code to get oneupn-formatted interactions

Under ``example/``, with sd.sp, usdb.int,

```
bash test-oneupn.sh
```

sd.one.upn.int will be generated, with 158 2-body proton-proton interaction elements, 0 neutron-neutron and 0 proton-neutron interaction elements.
If we feed such an interaction to bigstick, and compute Np = 2, Nn = 0 case in sd shell, then every isospin-T state of F18 will appear 2T+1 times.