How to construct a coarse-grain polymer model using EMC?

Hello,

I would like to simulate coarse-grained PMMA using the parameters reported in this reference paper. Systematic Method for Thermomechanically Consistent Coarse-Graining: A Universal Model for Methacrylate-Based Polymers - PubMed

Could you please advise on how to construct the polymer model in EMC using these tabulated coarse-grain parameters?

Thank you for your help.

Best regards,
Sanjib

Dear user,

The paper you are referring to has combinations of functional forms for which EMC does not have a predefined force field. You can, however, always use tabulated potentials instead, which means, that you would need to convert the needed functions to tables. EMC supports tabulated potentials. A parameter file consists of a force field definition section in which the potential type is defined,

# Define

ITEM	DEFINE

ANGLE		WARN
DENSITY		G/CC
ENERGY		KCAL/MOL
FFAPPLY		ALL
FFMODE		TABLE
FFTYPE		COARSE
LENGTH		ANGSTROM
NBONDED		1
TORSION		WARN
VERSION		1.0

ITEM	END	# DEFINE

Each parameter paragraph follows the same idea of representing parameters. For example, non-bonded parameters definitions follow

# Nonbond

ITEM	NONBOND

# type1 type2	T	t	TABLE	n	x0	dx	table

ch3	ch3	T	300	TABLE	628	2.35	0.02
187.01634698	168.4523544	151.84170742	136.96714462	123.63523077
...

where the ch3 ch3 definition at T = 300K with 628 entries, starting at r = 2.35 A, with steps of dr = 0.02. Bonded contributions are available for bonds, angles, torsions, and impropers using ITEM BOND, ITEM ANGLE, ITEM TORSION, and ITEM IMPROPER. For instance, the angle contribution would look like

# Nonbond

ITEM	NONBOND

# type1 type2	type3	T	t	TABLE	n	x0	dx	table

ch3	ch2	ch3	T	300	TABLE	140		1.74		0.01
...

for the ch3 ch2 ch3angle. Units are in radians.

Hi Pieter,

Thank you for your detailed reply with instructions.

I encounter the following error when I run the polymer.esh script with the cg_Hsu.prm file. The parameter file contains data in tabular form, separated by tabs.

Error: core/parse.c:169 ParseUnary:

   Illegal expression termination in line 1 of input

   Program aborted.

I have attached the relevant files (polymer.esh and cg_Hsu.prm) below for your reference. Could you please have a look and let me know if you find any issues with them?

There are two types of beads (A & B) in this coarse grain parameter set, and it doesn’t consider all types of angle and torsion interactions (for instance B-A-A-A is missing). Should I need to specify the missing interactions as zero in the parameter file?

Best regards,

Sanjib

#!/usr/bin/env emc.pl

# Options

ITEM		OPTIONS

project		cg
field		cg/cg_Hsu
auto		false
ntotal		1000
density		1.0
replace		true
build_dir	.
sample		volume
profile		density:true
units		real	

ITEM		END

# Groups

ITEM		GROUPS

a           *A([B])*,1,a:2
t           *B,1,a:1,1,a:2

ITEM		END

# Clusters

ITEM		CLUSTERS

poly		random,1

ITEM		END

# Polymers

ITEM		POLYMERS

poly
1		a,10,t,2

ITEM		END
#################################
#  CG_Hsu interaction parameters
#################################

# Force field definition

ITEM	DEFINE

ANGLE		WARN
CUTOFF	    15.0
DENSITY		G/CC
ENERGY		KCAL/MOL
FFAPPLY		ALL
FFMODE		TABLE
FFTYPE		COARSE
LENGTH		ANGSTROM
NBONDED		1
# INNER	    12.0
TORSION		WARN
VERSION		1.0

ITEM	END

# Masses

ITEM	MASS

# type	mass	name	ncons	charge	comment

A	85	A	2	0	Large_Bead
B	15	B	1	0	Small_Bead

ITEM	END

# Typing equivalences

ITEM	EQUIVALENCE

# type	pair	incr	bond	angle	torsion	improper

A	A	A	A	A	A	A
B	B	B	B	B	B	B

ITEM	END

# Nonbonded parameters

ITEM	NONBOND

# type1 type2	T	t	TABLE	n	x0	dx	table

A	A	T	300	TABLE	106	4.50	0.10
15.5575	11.2288	8.0530	5.7180	3.9992	2.7337	1.8031	1.1203	0.6216	0.2599	0.0000	-0.1839	-0.3114	-0.3968	-0.4512	-0.4826	-0.4972	-0.4997	-0.4934	-0.4811	-0.4646	-0.4455	-0.4247	-0.4032	-0.3814	-0.3598	-0.3388	-0.3184	-0.2989	-0.2803	-0.2627	-0.2460	-0.2303	-0.2156	-0.2018	-0.1889	-0.1768	-0.1655	-0.1550	-0.1452	-0.1360	-0.1275	-0.1195	-0.1121	-0.1052	-0.0987	-0.0927	-0.0871	-0.0819	-0.0770	-0.0725	-0.0682	-0.0643	-0.0605	-0.0571	-0.0538	-0.0508	-0.0480	-0.0453	-0.0428	-0.0405	-0.0383	-0.0362	-0.0343	-0.0325	-0.0308	-0.0292	-0.0277	-0.0262	-0.0249	-0.0236	-0.0225	-0.0213	-0.0203	-0.0193	-0.0184	-0.0175	-0.0166	-0.0159	-0.0151	-0.0144	-0.0137	-0.0131	-0.0125	-0.0119	-0.0114	-0.0109	-0.0104	-0.0100	-0.0095	-0.0091	-0.0087	-0.0083	-0.0080	-0.0076	-0.0073	-0.0070	-0.0067	-0.0065	-0.0062	-0.0059	-0.0057	-0.0055	-0.0053	-0.0050	-0.0048	
B	B	T	300	TABLE	112	3.90	0.10
14.22838958	8.961191336	5.366189635	2.921432031	1.270873793	0.170072614	-0.549751803	-1.005711768	-1.279401518	-1.427850896	-1.491008006	-1.496884914	-1.465118168	-1.409446551	-1.339444319	-1.261739369	-1.180872985	-1.099908853	-1.020865782	-0.945025925	-0.873154702	-0.805657876	-0.742693742	-0.684253188	-0.630216687	-0.580394709	-0.534556221	-0.492448584	-0.453811295	-0.418385289	-0.385919075	-0.356172609	-0.328919574	-0.303948535	-0.281063328	-0.260082912	-0.240840891	-0.223184806	-0.206975312	-0.19208528	-0.178398895	-0.165810754	-0.154225003	-0.143554514	-0.133720112	-0.124649861	-0.116278396	-0.10854632	-0.101399644	-0.094789277	-0.08867057	-0.08300289	-0.077749246	-0.072875942	-0.068352267	-0.064150214	-0.060244226	-0.056610969	-0.053229122	-0.050079192	-0.047143347	-0.044405261	-0.041849977	-0.039463786	-0.037234111	-0.035149405	-0.033199065	-0.031373341	-0.029663268	-0.028060593	-0.026557713	-0.025147624	-0.023823864	-0.022580471	-0.021411939	-0.02031318	-0.019279491	-0.018306518	-0.017390235	-0.016526907	-0.015713078	-0.014945537	-0.01422131	-0.013537631	-0.012891933	-0.012281828	-0.011705098	-0.011159676	-0.010643637	-0.010155192	-0.009692668	-0.009254508	-0.008839259	-0.008445564	-0.008072156	-0.00771785	-0.00738154	-0.007062189	-0.006758828	-0.006470551	-0.006196506	-0.005935897	-0.005687978	-0.005452047	-0.005227449	-0.005013564	-0.004809814	-0.004615654	-0.004430572	-0.004254087	-0.004085745	-0.00392512	
A	B	T	300	TABLE	109	4.20	0.10
16.0932	11.0599	7.4776	4.9266	3.1118	1.8247	0.9169	0.2823	-0.1553	-0.4508	-0.6440	-0.7638	-0.8310	-0.8610	-0.8650	-0.8510	-0.8249	-0.7911	-0.7527	-0.7118	-0.6700	-0.6285	-0.5879	-0.5489	-0.5116	-0.4763	-0.4431	-0.4120	-0.3829	-0.3558	-0.3307	-0.3073	-0.2856	-0.2656	-0.2470	-0.2298	-0.2139	-0.1992	-0.1856	-0.1730	-0.1614	-0.1506	-0.1406	-0.1314	-0.1228	-0.1149	-0.1075	-0.1007	-0.0943	-0.0884	-0.0830	-0.0779	-0.0732	-0.0687	-0.0646	-0.0608	-0.0572	-0.0539	-0.0508	-0.0479	-0.0452	-0.0427	-0.0403	-0.0381	-0.0360	-0.0340	-0.0322	-0.0305	-0.0289	-0.0274	-0.0259	-0.0246	-0.0233	-0.0222	-0.0210	-0.0200	-0.0190	-0.0181	-0.0172	-0.0164	-0.0156	-0.0148	-0.0141	-0.0135	-0.0128	-0.0122	-0.0117	-0.0112	-0.0107	-0.0102	-0.0097	-0.0093	-0.0089	-0.0085	-0.0081	-0.0078	-0.0075	-0.0071	-0.0068	-0.0066	-0.0063	-0.0060	-0.0058	-0.0055	-0.0053	-0.0051	-0.0049	-0.0047	-0.0045	

ITEM	END

# Bond parameters

ITEM	BOND

# type1 type2	T	t	TABLE	n	x0	dx	table
A	A	T	300	TABLE	76	2.20	0.02
56.723625	53.678625	50.717625	47.840625	45.047625	42.338625	39.713625	37.172625	34.715625	32.342625	30.053625	27.848625	25.727625	23.690625	21.737625	19.868625	18.083625	16.382625	14.765625	13.232625	11.783625	10.418625	9.137625	7.940625	6.827625	5.798625	4.853625	3.992625	3.215625	2.522625	1.913625	1.388625	0.947625	0.590625	0.317625	0.128625	0.023625	0.002625	0.065625	0.212625	0.443625	0.758625	1.157625	1.640625	2.207625	2.858625	3.593625	4.412625	5.315625	6.302625	7.373625	8.528625	9.767625	11.090625	12.497625	13.988625	15.563625	17.222625	18.965625	20.792625	22.703625	24.698625	26.777625	28.940625	31.187625	33.518625	35.933625	38.432625	41.015625	43.682625	46.433625	49.268625	52.187625	55.190625	58.277625	61.448625	
A	B	T	300	TABLE	76	3.00	0.02
17.25794504	16.22477384	15.22349064	14.25409544	13.31658824	12.41096904	11.53723784	10.69539464	9.88543944	9.10737224	8.36119304	7.64690184	6.96449864	6.31398344	5.69535624	5.10861704	4.55376584	4.03080264	3.53972744	3.08054024	2.65324104	2.25782984	1.89430664	1.56267144	1.26292424	0.99506504	0.75909384	0.55501064	0.38281544	0.24250824	0.13408904	0.05755784	0.01291464	0.00015944	0.01929224	0.07031304	0.15322184	0.26801864	0.41470344	0.59327624	0.80373704	1.04608584	1.32032264	1.62644744	1.96446024	2.33436104	2.73614984	3.16982664	3.63539144	4.13284424	4.66218504	5.22341384	5.81653064	6.44153544	7.09842824	7.78720904	8.50787784	9.26043464	10.04487944	10.86121224	11.70943304	12.58954184	13.50153864	14.44542344	15.42119624	16.42885704	17.46840584	18.53984264	19.64316744	20.77838024	21.94548104	23.14446984	24.37534664	25.63811144	26.93276424	28.25930504	

ITEM	END

# Angle parameters

ITEM	ANGLE

# type1 type2 type3	T	t	TABLE	n	x0	dx	table
A	A	A	T	300	TABLE	316	0.00	0.01
98.445916	97.536856	96.632112	95.731683	94.83557	93.943771	93.056287	92.173119	91.294265	90.419727	89.549504	88.683596	87.822003	86.964725	86.111762	85.263114	84.418782	83.578764	82.743062	81.911674	81.084602	80.261845	79.443403	78.629276	77.819464	77.013968	76.212786	75.41592	74.623368	73.835132	73.051211	72.271605	71.496314	70.725338	69.958677	69.196331	68.438301	67.684585	66.935185	66.190099	65.449329	64.712874	63.980734	63.252909	62.529399	61.810205	61.095325	60.384761	59.678511	58.976577	58.278958	57.585654	56.896665	56.211991	55.531632	54.855588	54.18386	53.516446	52.853348	52.194564	51.540096	50.889943	50.244105	49.602582	48.965374	48.332482	47.703904	47.079642	46.459694	45.844062	45.232745	44.625743	44.023056	43.424684	42.830627	42.240885	41.655459	41.074347	40.497551	39.92507	39.356903	38.793052	38.233516	37.678295	37.12739	36.580799	36.038523	35.500563	34.966918	34.437587	33.912572	33.391872	32.875487	32.363417	31.855662	31.352223	30.853098	30.358289	29.867794	29.381615	28.899751	28.422202	27.948968	27.480049	27.015445	26.555156	26.099183	25.647524	25.200181	24.757153	24.31844	23.884041	23.453958	23.028191	22.606738	22.1896	21.776778	21.36827	20.964078	20.564201	20.168638	19.777391	19.390459	19.007842	18.629541	18.255554	17.885882	17.520526	17.159485	16.802758	16.450347	16.102251	15.75847	15.419004	15.083854	14.753018	14.426497	14.104292	13.786401	13.472826	13.163566	12.858621	12.557991	12.261676	11.969676	11.681992	11.398622	11.119568	10.844828	10.574404	10.308295	10.046501	9.789022	9.535858	9.287009	9.042476	8.802257	8.566354	8.334765	8.107492	7.884534	7.665891	7.451563	7.24155	7.035852	6.83447	6.637402	6.44465	6.256212	6.07209	5.892283	5.716791	5.545614	5.378752	5.216205	5.057974	4.904057	4.754456	4.609169	4.468198	4.331542	4.199201	4.071175	3.947464	3.828068	3.712988	3.602222	3.495772	3.393637	3.295816	3.202311	3.113121	3.028246	2.947686	2.871442	2.799512	2.731897	2.668598	2.609614	2.554944	2.50459	2.458551	2.416827	2.379418	2.346325	2.317546	2.293083	2.272934	2.257101	2.245583	2.238379	2.235491	2.236918	2.242661	2.252718	2.26709	2.285778	2.30878	2.336098	2.367731	2.403679	2.443942	2.48852	2.537413	2.590621	2.648145	2.709983	2.776137	2.846605	2.921389	3.000488	3.083902	3.171631	3.263675	3.360034	3.460708	3.565697	3.675001	3.788618	3.906546	4.028782	4.155313	4.286119	4.421147	4.560291	4.703314	4.849725	4.99849	5.147508	5.292633	5.426159	5.535192	5.601955	5.609437	5.551601	5.43825	5.288252	5.1196	4.944898	4.771654	4.603986	4.444062	4.293006	4.1514	4.019544	3.897595	3.785636	3.68371	3.59184	3.510038	3.438311	3.376663	3.325096	3.283611	3.252209	3.23089	3.219654	3.218502	3.227433	3.246448	3.275546	3.314727	3.363993	3.423341	3.492774	3.57229	3.661889	3.761572	3.871339	3.991189	4.121123	4.261141	4.411242	4.571426	4.741695	4.922046	5.112482	5.313	5.523603	5.744289	5.975059	6.215912	6.466849	6.727869	6.998973	7.280161	7.571432	7.872786	8.184225	8.505746	8.837352	9.179041	9.530813	9.89267	10.264609	10.646633	
A	A	B	T	300	TABLE	316	0.00	0.01
154.950973	152.06899	149.226545	146.423284	143.658853	140.932898	138.24507	135.595018	132.982395	130.406855	127.868053	125.365646	122.899293	120.468654	118.07339	115.713165	113.387643	111.096492	108.839379	106.615974	104.425948	102.268974	100.144727	98.052881	95.993115	93.965108	91.968541	90.003096	88.068456	86.164307	84.290337	82.446234	80.631687	78.84639	77.090035	75.362316	73.662932	71.991579	70.347958	68.731769	67.142716	65.580503	64.044836	62.535423	61.051972	59.594195	58.161804	56.754513	55.372037	54.014093	52.680401	51.37068	50.084652	48.82204	47.58257	46.365969	45.171963	44.000284	42.850662	41.72283	40.616524	39.531478	38.46743	37.42412	36.401289	35.398679	34.416034	33.453099	32.509622	31.585351	30.680037	29.793431	28.925287	28.075361	27.243408	26.429187	25.632458	24.852982	24.090521	23.344842	22.615709	21.902891	21.206156	20.525276	19.860023	19.210171	18.575496	17.955775	17.350786	16.76031	16.18413	15.622027	15.073789	14.539201	14.018051	13.510129	13.015228	12.533139	12.063658	11.60658	11.161704	10.728829	10.307755	9.898285	9.500224	9.113376	8.73755	8.372554	8.018198	7.674294	7.340657	7.017101	6.703442	6.3995	6.105094	5.820046	5.544178	5.277316	5.019285	4.769914	4.529031	4.296469	4.072058	3.855634	3.647032	3.44609	3.252646	3.066541	2.887617	2.715718	2.550688	2.392375	2.240626	2.095293	1.956226	1.823279	1.696306	1.575163	1.459709	1.349803	1.245306	1.14608	1.051991	0.962902	0.878683	0.799201	0.724328	0.653935	0.587896	0.526087	0.468383	0.414665	0.364811	0.318704	0.276226	0.237263	0.2017	0.169425	0.140329	0.114301	0.091236	0.071026	0.053569	0.03876	0.0265	0.016689	0.009229	0.004024	0.000978	0	0.000997	0.00388	0.00856	0.014951	0.022966	0.032524	0.043541	0.055937	0.069634	0.084553	0.10062	0.117761	0.135902	0.154972	0.174904	0.195628	0.217078	0.23919	0.261901	0.285149	0.308875	0.33302	0.357527	0.382342	0.407411	0.432681	0.458102	0.483627	0.509206	0.534796	0.560351	0.585829	0.61119	0.636394	0.661403	0.686181	0.710694	0.734908	0.758793	0.782318	0.805455	0.828177	0.85046	0.87228	0.893615	0.914445	0.93475	0.954514	0.973722	0.992359	1.010412	1.027871	1.044727	1.060972	1.0766	1.091606	1.105988	1.119744	1.132874	1.14538	1.157266	1.168536	1.179198	1.189259	1.198728	1.207618	1.215942	1.223712	1.230947	1.237663	1.243879	1.249617	1.254898	1.259747	1.26419	1.268252	1.271964	1.275355	1.278458	1.281305	1.283931	1.286374	1.288672	1.290864	1.292991	1.295097	1.297227	1.299425	1.30174	1.304222	1.306921	1.309889	1.313181	1.316851	1.320959	1.325561	1.330718	1.336493	1.342949	1.35015	1.358165	1.36706	1.376905	1.387773	1.399736	1.412868	1.427247	1.442949	1.460054	1.478642	1.498797	1.520603	1.544144	1.569509	1.596785	1.626065	1.657438	1.691	1.726845	1.76507	1.805772	1.849053	1.895014	1.943757	1.995388	2.050012	2.107737	2.168673	2.23293	2.300621	2.371861	2.446765	2.52545	2.608035	2.694641	2.78539	2.880404	2.979811	3.083735	3.192306	3.305654	3.423909	3.547206	3.675679	3.809463	3.948698	

ITEM	END

# Torsion parameters

ITEM	TORSION

# type1 type2 type3	type4	T	t	TABLE	n	x0	dx	table
A	A	A	A	T	300	TABLE	629	-3.14	0.01
3.51952	3.519539	3.519597	3.519693	3.519827	3.520001	3.520213	3.520465	3.520757	3.521089	3.521462	3.521875	3.522331	3.522829	3.52337	3.523955	3.524585	3.525261	3.525983	3.526752	3.52757	3.528437	3.529354	3.530324	3.531346	3.532421	3.533552	3.53474	3.535985	3.537288	3.538653	3.540079	3.541568	3.543121	3.54474	3.546427	3.548182	3.550007	3.551904	3.553874	3.555919	3.55804	3.560238	3.562515	3.564873	3.567313	3.569836	3.572443	3.575137	3.577919	3.58079	3.58375	3.586803	3.589948	3.593188	3.596523	3.599955	3.603484	3.607112	3.610841	3.61467	3.618601	3.622636	3.626774	3.631017	3.635365	3.63982	3.644381	3.64905	3.653827	3.658713	3.663707	3.668811	3.674025	3.679348	3.684781	3.690325	3.695978	3.701742	3.707616	3.713599	3.719692	3.725894	3.732205	3.738623	3.74515	3.751783	3.758523	3.765368	3.772318	3.779371	3.786527	3.793784	3.801141	3.808597	3.816151	3.823801	3.831546	3.839384	3.847313	3.855332	3.863439	3.871632	3.879909	3.888269	3.896708	3.905226	3.913819	3.922486	3.931225	3.940032	3.948905	3.957843	3.966842	3.9759	3.985014	3.994182	4.003401	4.012668	4.021979	4.031334	4.040728	4.050158	4.059622	4.069117	4.078639	4.088186	4.097755	4.107342	4.116944	4.126559	4.136183	4.145812	4.155445	4.165077	4.174706	4.184328	4.193941	4.20354	4.213124	4.222688	4.23223	4.241746	4.251234	4.260691	4.270113	4.279497	4.28884	4.29814	4.307394	4.316598	4.32575	4.334846	4.343885	4.352862	4.361777	4.370625	4.379405	4.388113	4.396748	4.405307	4.413786	4.422185	4.430501	4.438731	4.446873	4.454925	4.462886	4.470753	4.478524	4.486197	4.493771	4.501243	4.508613	4.515879	4.523038	4.53009	4.537034	4.543867	4.55059	4.5572	4.563696	4.570078	4.576345	4.582495	4.588529	4.594445	4.600242	4.605921	4.61148	4.616919	4.622239	4.627438	4.632517	4.637474	4.642312	4.647028	4.651624	4.6561	4.660456	4.664692	4.668808	4.672805	4.676684	4.680445	4.684088	4.687615	4.691027	4.694323	4.697505	4.700574	4.70353	4.706376	4.709111	4.711737	4.714256	4.716668	4.718975	4.721178	4.723278	4.725277	4.727177	4.728979	4.730684	4.732294	4.733811	4.735236	4.736571	4.737818	4.738978	4.740054	4.741046	4.741957	4.742789	4.743543	4.744222	4.744827	4.74536	4.745823	4.746218	4.746547	4.746812	4.747015	4.747158	4.747242	4.747271	4.747245	4.747167	4.747038	4.746862	4.746639	4.746371	4.746062	4.745712	4.745323	4.744898	4.744438	4.743945	4.743422	4.74287	4.74229	4.741686	4.741058	4.740409	4.739739	4.739052	4.738349	4.737631	4.7369	4.736158	4.735407	4.734647	4.733882	4.733112	4.732338	4.731564	4.730789	4.730015	4.729244	4.728478	4.727716	4.726962	4.726216	4.725479	4.724753	4.724039	4.723338	4.722651	4.721979	4.721324	4.720686	4.720066	4.719465	4.718884	4.718324	4.717787	4.717271	4.716779	4.716311	4.715868	4.71545	4.715058	4.714693	4.714354	4.714043	4.71376	4.713506	4.71328	4.713082	4.712915	4.712776	4.712667	4.712588	4.712539	4.71252	4.712531	4.712572	4.712643	4.712744	4.712874	4.713034	4.713223	4.713441	4.713688	4.713964	4.714267	4.714598	4.714956	4.715341	4.715752	4.716188	4.716649	4.717135	4.717644	4.718176	4.71873	4.719305	4.7199	4.720515	4.721148	4.721799	4.722467	4.72315	4.723847	4.724558	4.725281	4.726015	4.726758	4.72751	4.72827	4.729035	4.729805	4.730578	4.731353	4.732128	4.732901	4.733673	4.734439	4.735201	4.735954	4.736699	4.737433	4.738154	4.738862	4.739554	4.740228	4.740883	4.741517	4.742128	4.742715	4.743274	4.743806	4.744307	4.744776	4.745211	4.745609	4.74597	4.746291	4.74657	4.746806	4.746995	4.747137	4.747229	4.747269	4.747256	4.747187	4.74706	4.746874	4.746626	4.746314	4.745937	4.745493	4.744979	4.744394	4.743736	4.743002	4.742192	4.741302	4.740332	4.739279	4.738143	4.73692	4.735609	4.734208	4.732716	4.731132	4.729453	4.727677	4.725804	4.723832	4.72176	4.719585	4.717306	4.714923	4.712433	4.709837	4.707131	4.704316	4.70139	4.698351	4.6952	4.691935	4.688556	4.68506	4.681448	4.67772	4.673873	4.669908	4.665824	4.661621	4.657298	4.652855	4.648292	4.643608	4.638803	4.633878	4.628833	4.623666	4.61838	4.612973	4.607446	4.6018	4.596035	4.590151	4.58415	4.578031	4.571796	4.565445	4.55898	4.552401	4.545709	4.538905	4.531992	4.524969	4.517838	4.510602	4.50326	4.495815	4.488269	4.480622	4.472878	4.465037	4.457102	4.449074	4.440956	4.43275	4.424457	4.416081	4.407623	4.399086	4.390472	4.381783	4.373022	4.364193	4.355296	4.346335	4.337313	4.328232	4.319095	4.309905	4.300664	4.291377	4.282045	4.272671	4.263259	4.253812	4.244332	4.234823	4.225288	4.215729	4.206151	4.196556	4.186946	4.177326	4.167699	4.158067	4.148435	4.138804	4.129178	4.119561	4.109955	4.100363	4.090789	4.081236	4.071707	4.062204	4.052732	4.043292	4.033888	4.024522	4.015198	4.005919	3.996687	3.987505	3.978376	3.969303	3.960287	3.951333	3.942441	3.933616	3.924858	3.916172	3.907558	3.89902	3.890559	3.882177	3.873877	3.865661	3.85753	3.849487	3.841533	3.833671	3.8259	3.818224	3.810644	3.803161	3.795777	3.788492	3.781309	3.774228	3.76725	3.760376	3.753608	3.746945	3.74039	3.733941	3.727601	3.72137	3.715247	3.709234	3.70333	3.697537	3.691853	3.68628	3.680816	3.675463	3.67022	3.665086	3.660062	3.655147	3.65034	3.645642	3.641051	3.636567	3.63219	3.627919	3.623752	3.61969	3.61573	3.611873	3.608117	3.604462	3.600906	3.597448	3.594086	3.590821	3.58765	3.584572	3.581587	3.578692	3.575886	3.573168	3.570537	3.567991	3.565529	3.563149	3.56085	3.55863	3.556489	3.554423	3.552433	3.550516	3.548672	3.546898	3.545193	3.543555	3.541984	3.540478	3.539035	3.537654	3.536334	3.535073	3.53387	3.532724	3.531633	3.530597	3.529613	3.528682	3.527801	3.52697	3.526187	3.525453	3.524765	3.524123	3.523525	3.522972	3.522462	3.521995	3.52157	3.521186	3.520843	3.520541	3.520278	3.520055	3.519871	3.519726	3.519619	3.519551	3.519521	
B	A	A	B	T	300	TABLE	629	-3.14	0.01
4.1066	4.106618	4.106672	4.106763	4.10689	4.107055	4.107257	4.107498	4.107778	4.108099	4.108461	4.108866	4.109316	4.109811	4.110353	4.110945	4.111587	4.112282	4.113032	4.113839	4.114704	4.115631	4.116622	4.117679	4.118804	4.12	4.12127	4.122615	4.12404	4.125545	4.127135	4.12881	4.130576	4.132432	4.134384	4.136432	4.138579	4.140829	4.143183	4.145643	4.148213	4.150894	4.153688	4.156598	4.159625	4.162771	4.166038	4.169428	4.172941	4.17658	4.180345	4.184237	4.188258	4.192407	4.196685	4.201094	4.205632	4.210299	4.215096	4.220022	4.225076	4.230257	4.235565	4.240997	4.246552	4.252229	4.258026	4.26394	4.269968	4.276108	4.282357	4.288712	4.29517	4.301726	4.308377	4.315118	4.321946	4.328855	4.335842	4.342899	4.350024	4.357209	4.36445	4.371739	4.379072	4.386442	4.393841	4.401264	4.408704	4.416153	4.423605	4.431051	4.438484	4.445897	4.453281	4.460628	4.467931	4.475182	4.482371	4.48949	4.496532	4.503486	4.510346	4.517101	4.523744	4.530265	4.536656	4.542907	4.549011	4.554959	4.560741	4.566349	4.571775	4.577009	4.582044	4.586871	4.591482	4.595868	4.600022	4.603936	4.607601	4.611011	4.614158	4.617035	4.619634	4.621949	4.623973	4.625701	4.627125	4.62824	4.62904	4.629521	4.629676	4.629501	4.628992	4.628144	4.626954	4.625417	4.623532	4.621294	4.618701	4.615752	4.612444	4.608775	4.604745	4.600353	4.595597	4.590479	4.584999	4.579157	4.572954	4.566392	4.559472	4.552197	4.54457	4.536593	4.528269	4.519602	4.510597	4.501257	4.491588	4.481595	4.471283	4.460657	4.449725	4.438492	4.426966	4.415154	4.403063	4.3907	4.378075	4.365195	4.352069	4.338706	4.325115	4.311306	4.297287	4.28307	4.268664	4.25408	4.239328	4.224418	4.209362	4.194171	4.178855	4.163426	4.147895	4.132274	4.116575	4.100808	4.084986	4.069121	4.053223	4.037305	4.021379	4.005457	3.989549	3.973668	3.957826	3.942033	3.926302	3.910643	3.895068	3.879588	3.864214	3.848957	3.833827	3.818834	3.80399	3.789303	3.774784	3.760443	3.746288	3.732329	3.718574	3.705032	3.691711	3.678619	3.665764	3.653153	3.640793	3.62869	3.616852	3.605283	3.593989	3.582976	3.572249	3.561812	3.551668	3.541823	3.532278	3.523038	3.514105	3.50548	3.497166	3.489165	3.481476	3.474101	3.46704	3.460293	3.45386	3.447738	3.441928	3.436428	3.431234	3.426346	3.42176	3.417472	3.413481	3.409781	3.406369	3.40324	3.40039	3.397813	3.395503	3.393456	3.391666	3.390125	3.388828	3.387767	3.386937	3.386329	3.385937	3.385752	3.385767	3.385975	3.386366	3.386933	3.387667	3.38856	3.389603	3.390787	3.392105	3.393546	3.395103	3.396765	3.398526	3.400374	3.402302	3.404301	3.406362	3.408476	3.410634	3.412828	3.41505	3.41729	3.419542	3.421795	3.424044	3.426279	3.428494	3.43068	3.432831	3.43494	3.436999	3.439002	3.440944	3.442817	3.444616	3.446336	3.447972	3.449517	3.450969	3.452322	3.453573	3.454717	3.455752	3.456675	3.457482	3.458171	3.458741	3.45919	3.459517	3.45972	3.459799	3.459754	3.459584	3.459291	3.458876	3.458338	3.457681	3.456906	3.456015	3.45501	3.453895	3.452673	3.451347	3.449922	3.448402	3.44679	3.445093	3.443314	3.441461	3.439537	3.43755	3.435505	3.43341	3.43127	3.429092	3.426885	3.424654	3.422408	3.420155	3.417902	3.415658	3.413431	3.411228	3.409059	3.406932	3.404856	3.40284	3.400892	3.39902	3.397235	3.395545	3.393959	3.392485	3.391133	3.389912	3.388829	3.387894	3.387116	3.386503	3.386063	3.385805	3.385737	3.385866	3.386201	3.38675	3.387519	3.388516	3.389748	3.391222	3.392944	3.39492	3.397158	3.399661	3.402437	3.405489	3.408824	3.412445	3.416357	3.420563	3.425067	3.429873	3.434983	3.4404	3.446126	3.452162	3.458511	3.465172	3.472148	3.479437	3.48704	3.494957	3.503186	3.511726	3.520576	3.529733	3.539194	3.548958	3.559021	3.569379	3.580027	3.590963	3.602181	3.613675	3.625441	3.637472	3.649763	3.662306	3.675096	3.688124	3.701383	3.714865	3.728563	3.742468	3.75657	3.770862	3.785333	3.799975	3.814778	3.829731	3.844825	3.860048	3.875392	3.890844	3.906394	3.922031	3.937744	3.953521	3.969351	3.985223	4.001124	4.017044	4.03297	4.048891	4.064795	4.080671	4.096506	4.112288	4.128007	4.143651	4.159207	4.174665	4.190013	4.205239	4.220333	4.235284	4.25008	4.264711	4.279167	4.293436	4.307509	4.321377	4.335028	4.348454	4.361645	4.374593	4.387289	4.399724	4.411889	4.423779	4.435383	4.446696	4.457711	4.468421	4.478819	4.4889	4.498657	4.508087	4.517184	4.525943	4.534361	4.542433	4.550156	4.557527	4.564543	4.571203	4.577504	4.583444	4.589023	4.59424	4.599094	4.603585	4.607714	4.611481	4.614887	4.617934	4.620623	4.622958	4.624939	4.62657	4.627854	4.628795	4.629396	4.629661	4.629595	4.629203	4.628489	4.627459	4.626119	4.624473	4.622529	4.620293	4.61777	4.614968	4.611894	4.608555	4.604959	4.601112	4.597023	4.592699	4.588148	4.583379	4.5784	4.573219	4.567845	4.562285	4.55655	4.550647	4.544584	4.538372	4.532018	4.525532	4.518921	4.512196	4.505364	4.498434	4.491415	4.484316	4.477145	4.469911	4.462621	4.455285	4.44791	4.440504	4.433076	4.425633	4.418182	4.410732	4.403289	4.39586	4.388453	4.381075	4.373732	4.36643	4.359175	4.351974	4.344833	4.337756	4.33075	4.323819	4.316969	4.310203	4.303527	4.296945	4.29046	4.284077	4.277799	4.271629	4.265569	4.259624	4.253796	4.248086	4.242497	4.237031	4.23169	4.226474	4.221385	4.216424	4.211592	4.206889	4.202316	4.197873	4.193559	4.189374	4.185319	4.181392	4.177592	4.173919	4.170372	4.166949	4.163649	4.16047	4.15741	4.154469	4.151643	4.148932	4.146332	4.143842	4.141459	4.139182	4.137007	4.134932	4.132954	4.131072	4.129282	4.127582	4.125969	4.124441	4.122995	4.121628	4.120338	4.119123	4.117978	4.116903	4.115895	4.114951	4.114068	4.113246	4.112481	4.111771	4.111115	4.110509	4.109954	4.109446	4.108984	4.108567	4.108193	4.107861	4.10757	4.107319	4.107106	4.106931	4.106794	4.106693	4.106629	4.106601	

ITEM	END