cosmoxml

The purpose of cosmoxml (xcosmo in some versions) is to create an xml file which describes the system to be studied with the shell model. You do not need to use cosmoxml if you create xmls yourself. The cosmoxml reads database in ./interactions/int.xml, however you can have it read any other database file by using the filename as an argument e.g. "cosmoxml filename". Note if path is not setup properly and cosmoxml sees no interactions then you can also point the program to the int.xml file or copy int.xml to your local directory.

Example below shows how to create a job to study 40Ar using sdpf-nr interaction with rejection where protons are in sd and neutrons are in fp.  As a result of this you will get my40Ar.xml file.


> cosmoxml
	Reading int database from
	/Users/volya/work/CoSMo/interactions/int.xml
	Enter nucleus name: 40Ar //Here we type nucleus of interest
	22 18
	number of spaces 40

	...... all available model spaces are offered....

	you can use modelspace sdpf   16O - 80Zr

	.......

	select model space sdpf //Here we select the model space to use
	Available Hamiltonians

	.....all available interactions offered ......

	--------sdpf-nr--------
	PRC 63 (2001) 44316
	SDPF-NR INTERACTION
	.......
	select Hamiltonian sdpf-nr //Here we select interaction
	name this the system as my40Ar //Here we give a name to our file
	Add rejection (y/n)[n]<>y //Here we create first condition restricting 12 neutrons to sd shell
	Orbital: 0d5 type= pninclude (y/n)[n]y
	neutron, proton, or both (n, p, or pn): [pn]n
	Orbital: 0d3 type= pninclude (y/n)[n]y
	neutron, proton, or both (n, p, or pn): [pn]n
	Orbital: 1s1 type= pninclude (y/n)[n]y
	neutron, proton, or both (n, p, or pn): [pn]n
	Orbital: 0f7 type= pninclude (y/n)[n]n
	Orbital: 0f5 type= pninclude (y/n)[n]n
	Orbital: 1p3 type= pninclude (y/n)[n]n
	Orbital: 1p1 type= pninclude (y/n)[n]n
	Minimum N on selected orbits: [0]12
	Maximum N on selected orbits: [12]12
	Add rejection (y/n)[n]y //Here we create second condition requiring zero protons on fp shell
	Orbital: 0d5 type= pninclude (y/n)[n]n
	Orbital: 0d3 type= pninclude (y/n)[n]n
	Orbital: 1s1 type= pninclude (y/n)[n]n
	Orbital: 0f7 type= pninclude (y/n)[n]y
	neutron, proton, or both (n, p, or pn): [pn]p
	Orbital: 0f5 type= pninclude (y/n)[n]y
	neutron, proton, or both (n, p, or pn): [pn]p
	Orbital: 1p3 type= pninclude (y/n)[n]y
	neutron, proton, or both (n, p, or pn): [pn]p
	Orbital: 1p1 type= pninclude (y/n)[n]y
	neutron, proton, or both (n, p, or pn): [pn]p
	Minimum N on selected orbits: [0]0
	Maximum N on selected orbits: [20]0
	Add rejection (y/n)[n]n

Published on  December 1st, 2023