eff_z_Ar.h

/*
 * eff_z_Ar.h
 *
 *  Created on: Mar 26, 2020
 *      Author: Evan
 */

// Effective Z as a function of the radial distance from a point nucleus based on numerical Dirac-Fock electron densities from ELSEPA.

#ifndef EFF_Z_AR_H_
#define EFF_Z_AR_H_

double radius_ar(int index){
	double r_list[677]{
		3.07661E-5,
		3.15449E-5,
		3.23435E-5,
		3.31622E-5,
		3.40017E-5,
		3.48625E-5,
		3.5745E-5,
		3.66499E-5,
		3.75777E-5,
		3.8529E-5,
		3.95044E-5,
		4.05044E-5,
		4.15298E-5,
		4.25811E-5,
		4.36591E-5,
		4.47643E-5,
		4.58975E-5,
		4.70594E-5,
		4.82507E-5,
		4.94722E-5,
		5.07246E-5,
		5.20087E-5,
		5.33253E-5,
		5.46752E-5,
		5.60593E-5,
		5.74784E-5,
		5.89335E-5,
		6.04254E-5,
		6.1955E-5,
		6.35234E-5,
		6.51315E-5,
		6.67803E-5,
		6.84709E-5,
		7.02042E-5,
		7.19814E-5,
		7.38036E-5,
		7.5672E-5,
		7.75876E-5,
		7.95517E-5,
		8.15656E-5,
		8.36304E-5,
		8.57475E-5,
		8.79182E-5,
		9.01438E-5,
		9.24258E-5,
		9.47655E-5,
		9.71645E-5,
		9.96242E-5,
		1.02146E-4,
		1.04732E-4,
		1.07383E-4,
		1.10102E-4,
		1.12889E-4,
		1.15747E-4,
		1.18677E-4,
		1.21681E-4,
		1.24761E-4,
		1.2792E-4,
		1.31158E-4,
		1.34478E-4,
		1.37882E-4,
		1.41373E-4,
		1.44952E-4,
		1.48621E-4,
		1.52383E-4,
		1.56241E-4,
		1.60196E-4,
		1.64252E-4,
		1.68409E-4,
		1.72673E-4,
		1.77044E-4,
		1.81526E-4,
		1.86121E-4,
		1.90833E-4,
		1.95663E-4,
		2.00617E-4,
		2.05695E-4,
		2.10902E-4,
		2.16241E-4,
		2.21715E-4,
		2.27328E-4,
		2.33082E-4,
		2.38983E-4,
		2.45032E-4,
		2.51235E-4,
		2.57595E-4,
		2.64116E-4,
		2.70802E-4,
		2.77657E-4,
		2.84686E-4,
		2.91893E-4,
		2.99282E-4,
		3.06858E-4,
		3.14626E-4,
		3.2259E-4,
		3.30756E-4,
		3.39129E-4,
		3.47714E-4,
		3.56516E-4,
		3.65541E-4,
		3.74794E-4,
		3.84282E-4,
		3.9401E-4,
		4.03984E-4,
		4.1421E-4,
		4.24696E-4,
		4.35446E-4,
		4.46469E-4,
		4.57771E-4,
		4.69359E-4,
		4.81241E-4,
		4.93423E-4,
		5.05913E-4,
		5.1872E-4,
		5.3185E-4,
		5.45313E-4,
		5.59117E-4,
		5.73271E-4,
		5.87782E-4,
		6.02661E-4,
		6.17917E-4,
		6.33558E-4,
		6.49596E-4,
		6.66039E-4,
		6.82899E-4,
		7.00185E-4,
		7.17909E-4,
		7.36082E-4,
		7.54715E-4,
		7.73819E-4,
		7.93407E-4,
		8.1349E-4,
		8.34082E-4,
		8.55195E-4,
		8.76843E-4,
		8.99038E-4,
		9.21795E-4,
		9.45128E-4,
		9.69052E-4,
		9.93581E-4,
		0.00102,
		0.00104,
		0.00107,
		0.0011,
		0.00113,
		0.00115,
		0.00118,
		0.00121,
		0.00124,
		0.00128,
		0.00131,
		0.00134,
		0.00138,
		0.00141,
		0.00145,
		0.00148,
		0.00152,
		0.00156,
		0.0016,
		0.00164,
		0.00168,
		0.00172,
		0.00177,
		0.00181,
		0.00186,
		0.0019,
		0.00195,
		0.002,
		0.00205,
		0.0021,
		0.00216,
		0.00221,
		0.00227,
		0.00232,
		0.00238,
		0.00244,
		0.00251,
		0.00257,
		0.00263,
		0.0027,
		0.00277,
		0.00284,
		0.00291,
		0.00298,
		0.00306,
		0.00314,
		0.00322,
		0.0033,
		0.00338,
		0.00347,
		0.00355,
		0.00364,
		0.00374,
		0.00383,
		0.00393,
		0.00403,
		0.00413,
		0.00423,
		0.00434,
		0.00445,
		0.00456,
		0.00468,
		0.0048,
		0.00492,
		0.00504,
		0.00517,
		0.0053,
		0.00544,
		0.00557,
		0.00571,
		0.00586,
		0.00601,
		0.00616,
		0.00632,
		0.00648,
		0.00664,
		0.00681,
		0.00698,
		0.00716,
		0.00734,
		0.00752,
		0.00771,
		0.00791,
		0.00811,
		0.00831,
		0.00852,
		0.00874,
		0.00896,
		0.00919,
		0.00942,
		0.00966,
		0.0099,
		0.01015,
		0.01041,
		0.01067,
		0.01094,
		0.01122,
		0.0115,
		0.01179,
		0.01209,
		0.0124,
		0.01271,
		0.01303,
		0.01336,
		0.0137,
		0.01404,
		0.0144,
		0.01476,
		0.01514,
		0.01552,
		0.01591,
		0.01631,
		0.01673,
		0.01715,
		0.01758,
		0.01803,
		0.01848,
		0.01895,
		0.01943,
		0.01992,
		0.02042,
		0.02094,
		0.02147,
		0.02201,
		0.02257,
		0.02313,
		0.02372,
		0.02432,
		0.02493,
		0.02556,
		0.02621,
		0.02687,
		0.02755,
		0.02824,
		0.02896,
		0.02969,
		0.03044,
		0.0312,
		0.03199,
		0.0328,
		0.03363,
		0.03447,
		0.03534,
		0.03624,
		0.03715,
		0.03809,
		0.03905,
		0.04003,
		0.04104,
		0.04207,
		0.04314,
		0.04422,
		0.04534,
		0.04648,
		0.04765,
		0.04885,
		0.05008,
		0.05134,
		0.05264,
		0.05396,
		0.05532,
		0.05671,
		0.05814,
		0.0596,
		0.0611,
		0.06264,
		0.06422,
		0.06583,
		0.06749,
		0.06918,
		0.07092,
		0.0727,
		0.07453,
		0.0764,
		0.07832,
		0.08029,
		0.0823,
		0.08437,
		0.08649,
		0.08866,
		0.09088,
		0.09316,
		0.0955,
		0.09789,
		0.10035,
		0.10286,
		0.10544,
		0.10808,
		0.11078,
		0.11356,
		0.1164,
		0.11931,
		0.12229,
		0.12535,
		0.12848,
		0.1317,
		0.13498,
		0.13836,
		0.14181,
		0.14535,
		0.14897,
		0.15269,
		0.15649,
		0.16039,
		0.16439,
		0.16848,
		0.17267,
		0.17697,
		0.18137,
		0.18587,
		0.19049,
		0.19522,
		0.20007,
		0.20503,
		0.21011,
		0.21532,
		0.22065,
		0.22612,
		0.23171,
		0.23744,
		0.24331,
		0.24932,
		0.25547,
		0.26177,
		0.26823,
		0.27484,
		0.2816,
		0.28853,
		0.29563,
		0.30289,
		0.31033,
		0.31794,
		0.32574,
		0.33372,
		0.34188,
		0.35025,
		0.35881,
		0.36757,
		0.37653,
		0.38571,
		0.3951,
		0.40472,
		0.41455,
		0.42462,
		0.43492,
		0.44546,
		0.45625,
		0.46728,
		0.47857,
		0.49012,
		0.50193,
		0.51402,
		0.52638,
		0.53902,
		0.55195,
		0.56518,
		0.5787,
		0.59253,
		0.60667,
		0.62113,
		0.63592,
		0.65103,
		0.66648,
		0.68227,
		0.69842,
		0.71492,
		0.73178,
		0.74901,
		0.76662,
		0.78462,
		0.803,
		0.82179,
		0.84098,
		0.86058,
		0.8806,
		0.90105,
		0.92193,
		0.94325,
		0.96503,
		0.98726,
		1.00996,
		1.03313,
		1.05678,
		1.08092,
		1.10556,
		1.1307,
		1.15635,
		1.18253,
		1.20923,
		1.23647,
		1.26425,
		1.29259,
		1.32149,
		1.35095,
		1.381,
		1.41163,
		1.44285,
		1.47467,
		1.50711,
		1.54016,
		1.57384,
		1.60816,
		1.64311,
		1.67872,
		1.71499,
		1.75192,
		1.78953,
		1.82782,
		1.8668,
		1.90648,
		1.94686,
		1.98796,
		2.02978,
		2.07233,
		2.11561,
		2.15964,
		2.20442,
		2.24995,
		2.29625,
		2.34333,
		2.39118,
		2.43982,
		2.48925,
		2.53947,
		2.59051,
		2.64235,
		2.69501,
		2.7485,
		2.80281,
		2.85796,
		2.91395,
		2.97079,
		3.02847,
		3.08701,
		3.14641,
		3.20668,
		3.26782,
		3.32983,
		3.39272,
		3.4565,
		3.52116,
		3.58671,
		3.65315,
		3.72049,
		3.78873,
		3.85788,
		3.92793,
		3.99889,
		4.07076,
		4.14354,
		4.21724,
		4.29186,
		4.36739,
		4.44385,
		4.52122,
		4.59952,
		4.67875,
		4.75889,
		4.83997,
		4.92197,
		5.00489,
		5.08874,
		5.17352,
		5.25922,
		5.34585,
		5.4334,
		5.52187,
		5.61127,
		5.70159,
		5.79283,
		5.88499,
		5.97807,
		6.07206,
		6.16697,
		6.26279,
		6.35952,
		6.45716,
		6.5557,
		6.65515,
		6.75549,
		6.85673,
		6.95887,
		7.0619,
		7.16582,
		7.27062,
		7.37631,
		7.48288,
		7.59032,
		7.69863,
		7.80781,
		7.91785,
		8.02875,
		8.14051,
		8.25313,
		8.36659,
		8.48089,
		8.59604,
		8.71202,
		8.82883,
		8.94646,
		9.06492,
		9.1842,
		9.30429,
		9.42519,
		9.54689,
		9.66939,
		9.79269,
		9.91677,
		10.04164,
		10.16729,
		10.29371,
		10.42091,
		10.54887,
		10.67758,
		10.80706,
		10.93728,
		11.06825,
		11.19995,
		11.33239,
		11.46557,
		11.59946,
		11.73408,
		11.86941,
		12.00544,
		12.14219,
		12.27963,
		12.41777,
		12.55659,
		12.6961,
		12.83629,
		12.97715,
		13.11868,
		13.26087,
		13.40372,
		13.54723,
		13.69138,
		13.83618,
		13.98162,
		14.12768,
		14.27438,
		14.4217,
		14.56965,
		14.7182,
		14.86736,
		15.01713,
		15.1675,
		15.31846,
		15.47001,
		15.62215,
		15.77487,
		15.92816,
		16.08203,
		16.23646,
		16.39145,
		16.547,
		16.7031,
		16.85975,
		17.01695,
		17.17468,
		17.33295,
		17.49175,
		17.65108,
		17.81093,
		17.97129,
		18.13217,
		18.29356,
		18.45545,
		18.61784,
		18.78073,
		18.94411,
		19.10798,
		19.27234,
		19.43717,
		19.60248,
		19.76827,
		19.93452,
		20.10124,
		20.26841,
		20.43605,
		20.60413,
		20.77267,
		20.94165,
		21.11107,
		21.28093,
		21.45123,
		21.62196,
		21.79311,
		21.96469,
		22.13669,
		22.3091,
		22.48193,
		22.65517,
		22.82881,
		23.00286,
		23.17731,
		23.35216,
		23.52739,
		23.70302,
		23.87904,
		24.05544,
		24.23222,
		24.40938,
		24.58691,
		24.76481,
		24.94308,
		25.12172,
		25.30072,
		25.48008,
		25.65979,
		25.83986,
		26.02028,
		26.20105,
		26.38217,
		26.56362,
		26.74542,
		26.92755,
		27.11002,
		27.29281,
		27.47594,
		27.65939,
		27.84317,
		28.02727,
		28.21168,
		28.39642,
		28.58146,
		28.76682,
		28.95249,
		29.13846,
		29.32473,
		29.51131,
		29.69819,
		29.88536,
		30.07283,
		30.26059,
		30.44863,
		30.63697,
		30.82559};
	return r_list[index];
	}


double effz_ar(int index){
	double z_list[677]{
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		18,
		17.99999,
		17.99999,
		17.99999,
		17.99999,
		17.99999,
		17.99999,
		17.99999,
		17.99999,
		17.99999,
		17.99999,
		17.99999,
		17.99999,
		17.99999,
		17.99999,
		17.99999,
		17.99998,
		17.99998,
		17.99998,
		17.99998,
		17.99998,
		17.99998,
		17.99998,
		17.99997,
		17.99997,
		17.99997,
		17.99997,
		17.99997,
		17.99996,
		17.99996,
		17.99996,
		17.99995,
		17.99995,
		17.99995,
		17.99994,
		17.99994,
		17.99993,
		17.99993,
		17.99992,
		17.99992,
		17.99991,
		17.9999,
		17.9999,
		17.99989,
		17.99988,
		17.99987,
		17.99986,
		17.99985,
		17.99984,
		17.99982,
		17.99981,
		17.9998,
		17.99978,
		17.99977,
		17.99975,
		17.99973,
		17.99971,
		17.99969,
		17.99966,
		17.99964,
		17.99961,
		17.99958,
		17.99955,
		17.99952,
		17.99948,
		17.99944,
		17.9994,
		17.99935,
		17.9993,
		17.99925,
		17.9992,
		17.99914,
		17.99907,
		17.999,
		17.99893,
		17.99885,
		17.99876,
		17.99867,
		17.99857,
		17.99847,
		17.99835,
		17.99823,
		17.9981,
		17.99796,
		17.99781,
		17.99765,
		17.99748,
		17.99729,
		17.99709,
		17.99688,
		17.99665,
		17.99641,
		17.99615,
		17.99587,
		17.99557,
		17.99525,
		17.9949,
		17.99454,
		17.99414,
		17.99372,
		17.99327,
		17.99279,
		17.99227,
		17.99172,
		17.99113,
		17.9905,
		17.98983,
		17.98911,
		17.98834,
		17.98752,
		17.98664,
		17.98571,
		17.98471,
		17.98364,
		17.9825,
		17.98129,
		17.98,
		17.97862,
		17.97715,
		17.97559,
		17.97392,
		17.97214,
		17.97026,
		17.96824,
		17.96611,
		17.96383,
		17.96141,
		17.95884,
		17.95611,
		17.95321,
		17.95014,
		17.94687,
		17.9434,
		17.93973,
		17.93583,
		17.9317,
		17.92733,
		17.9227,
		17.91781,
		17.91262,
		17.90715,
		17.90136,
		17.89524,
		17.88879,
		17.88197,
		17.87479,
		17.86721,
		17.85924,
		17.85083,
		17.84199,
		17.83269,
		17.82292,
		17.81265,
		17.80187,
		17.79055,
		17.77869,
		17.76627,
		17.75326,
		17.73964,
		17.72541,
		17.71053,
		17.695,
		17.6788,
		17.66191,
		17.64432,
		17.62601,
		17.60697,
		17.58718,
		17.56664,
		17.54533,
		17.52325,
		17.50039,
		17.47673,
		17.45229,
		17.42706,
		17.40102,
		17.3742,
		17.34658,
		17.31818,
		17.28901,
		17.25906,
		17.22836,
		17.19692,
		17.16475,
		17.13187,
		17.0983,
		17.06406,
		17.02918,
		16.99368,
		16.95759,
		16.92093,
		16.88372,
		16.84601,
		16.80781,
		16.76914,
		16.73005,
		16.69055,
		16.65066,
		16.6104,
		16.56979,
		16.52885,
		16.48757,
		16.44598,
		16.40405,
		16.36179,
		16.31917,
		16.27618,
		16.23278,
		16.18893,
		16.14458,
		16.09968,
		16.05415,
		16.00791,
		15.96087,
		15.91292,
		15.86396,
		15.81387,
		15.76249,
		15.7097,
		15.65534,
		15.59923,
		15.54122,
		15.48113,
		15.41876,
		15.35393,
		15.28646,
		15.21614,
		15.14279,
		15.06622,
		14.98624,
		14.90268,
		14.81536,
		14.72414,
		14.62887,
		14.52942,
		14.4257,
		14.3176,
		14.20508,
		14.0881,
		13.96664,
		13.84073,
		13.71041,
		13.57578,
		13.43695,
		13.29407,
		13.14732,
		12.99692,
		12.84313,
		12.68623,
		12.52654,
		12.3644,
		12.20019,
		12.0343,
		11.86715,
		11.69919,
		11.53086,
		11.36262,
		11.19494,
		11.0283,
		10.86315,
		10.69994,
		10.53912,
		10.3811,
		10.22628,
		10.07502,
		9.92764,
		9.78444,
		9.64565,
		9.51147,
		9.38205,
		9.25748,
		9.13778,
		9.02295,
		8.91289,
		8.80748,
		8.70652,
		8.60976,
		8.5169,
		8.4276,
		8.34146,
		8.25805,
		8.1769,
		8.09752,
		8.01938,
		7.94195,
		7.86468,
		7.78703,
		7.70846,
		7.62843,
		7.54644,
		7.46198,
		7.37461,
		7.28391,
		7.18949,
		7.09102,
		6.98821,
		6.88085,
		6.76875,
		6.65178,
		6.5299,
		6.40309,
		6.2714,
		6.13495,
		5.99389,
		5.84844,
		5.69884,
		5.5454,
		5.38847,
		5.2284,
		5.06561,
		4.90052,
		4.73359,
		4.56526,
		4.39602,
		4.22635,
		4.0567,
		3.88757,
		3.71941,
		3.55266,
		3.38777,
		3.22515,
		3.06518,
		2.90823,
		2.75464,
		2.60472,
		2.45874,
		2.31697,
		2.1796,
		2.04683,
		1.91881,
		1.79567,
		1.67749,
		1.56434,
		1.45626,
		1.35326,
		1.25531,
		1.16239,
		1.07443,
		0.99134,
		0.91304,
		0.83941,
		0.77032,
		0.70564,
		0.64521,
		0.58887,
		0.53646,
		0.48781,
		0.44275,
		0.40111,
		0.36269,
		0.32734,
		0.29486,
		0.2651,
		0.23789,
		0.21305,
		0.19043,
		0.16987,
		0.15123,
		0.13436,
		0.11912,
		0.1054,
		0.09305,
		0.08198,
		0.07206,
		0.0632,
		0.0553,
		0.04827,
		0.04203,
		0.03651,
		0.03162,
		0.02731,
		0.02352,
		0.02019,
		0.01727,
		0.01472,
		0.01249,
		0.01055,
		0.00886,
		0.0074,
		0.00614,
		0.00505,
		0.00411,
		0.0033,
		0.00261,
		0.00202,
		0.00151,
		0.00108,
		7.11958E-4,
		4.01179E-4,
		1.38132E-4,
		-8.40299E-5,
		-2.71256E-4,
		-4.28697E-4,
		-5.60805E-4,
		-6.71418E-4,
		-7.63833E-4,
		-8.40878E-4,
		-9.04972E-4,
		-9.58177E-4,
		-0.001,
		-0.00104,
		-0.00107,
		-0.00109,
		-0.00111,
		-0.00113,
		-0.00114,
		-0.00116,
		-0.00116,
		-0.00117,
		-0.00118,
		-0.00118,
		-0.00119,
		-0.00119,
		-0.00119,
		-0.00119,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012,
		-0.0012};
	return z_list[index];
}


double get_z_ar(double r){
	 double Z;
	 int index2 = 0;
	 int pass = 0;
	 if (r > 5.97807){
		 Z = 0;
	 }
	 else{
		 while (pass == 0){
		 			 if(r > radius_ar(index2)){
		 				 index2 = index2 + 10;
		 					}
		 			 else{
		 				 Z = effz_ar(index2);
		 				 pass = 1;
		 					}
		 				}
	 }
	 return Z;
}





#endif /* EFF_Z_AR_H_ */