CHAPTER 3

Binarity in A Uniform Cluster for Flat / Steep IMFs

3.1 Introduction

In Astronomy, there is two ways of study for the investigation over a group of objects. In the first case, a separate study of individuals in a group is encouraged. Second case deals with the distribution of the members of the group under study. Since binary system is a special case of multiple systems, in order to study about the formation of binaries it is significant to know about stellar multiplicity 1. Hence in this dissertation we study stellar multiplicity as a function of primary mass in stellar clusters with flatter / steeper IMFs than Salpeter.

3.2 Binary Formation by Capture- The Condition for Binding

Formation of a single binary system may happen in accordance with three basic theories.

Capture theory: According to the theory, two independent stars may deviate from their own path and will orbit around their common centre of mass under suitable conditions.

Fission theory: The theory suggests that, a single star in its initial state or afterwards may breaks into two to form a binary star system.

Independent Nuclei Theory: If the constitute of the nebulae condense to form two nuclei separated by appropriate distance to become a binary star system this process takes place as per the independent nuclei theory.

In this dissertation we are considering binaries formed by capture process. Capture theory was put forward by Dr. G. Johnston Stoney. He says that if two stars come closer to each other by their usual motions then one of the three things will happen 2. Either they will pass each other without interrupting the visibility of the companion, or they caught together and look like the fragments of a single star, or lastly the stars may come face to face without affecting the visibility of each other. In the third case the stars will revolve in elliptical orbits after the event. Since density is not uniform in the atmosphere of each star the resistance act as forces. This force may act tangentially and reduce the periastron distance or it may act normally and increase the distance 3. The force dominates each other in a special way and one star will fall into the other. This process results in the formation of binaries 2.

Total energy of a system is the sum of its kinetic energy and potential energy. Imagine two stars are seen to be isolated from other stars. Let the total energy of the combined system is Etotal , then

If Etotal < 0, stars are bounded or stars forms a binary system.

If Etotal > 0, stars are unbounded and excess energy should be dissipated in any way to form a binary system.

Another approach to the boundedness which is used as the binding condition in this dissertation is given below,

Let mi and mj are the masses of two stars which are expected to form a binary.

For such a system to be bounded its potential energy should be greater than its kinetic energy.

Gmimjrij > miVi-Vj2Where mi > mjG is the gravitational constant

Vi, Vj – velocities of mi mass star and mj mass star respectively.

ri , rj – radius of mi mass star and mj mass star respectively.

rij = | ri – rj |

3.3 Stellar Multiplicity as a Function of Primary Mass

In order to study the stellar multiplicity as a function of primary mass, we start with the process of random sampling from Salpeter function. Firstly we designed a program in the R language to draw 1000 random numbers between 0 and 1 from a uniform distribution and the value of unknown mass can be calculated using the equation 2.4. Salpeter’s slope also obtained from the program. Second step consists of randomly drawing velocity from a Gaussian distribution with mean zero and varying standard deviation from 0.1 to 1.3. Our aim is to study about binary systems in stellar clusters with IMF flatter / steeper than Salpeter. For this study we developed a program for a cluster having Salpeter IMF with a slope of 2.3. Then compared with clusters having steeper (slope = 3.3) and flatter IMF (slope=1.3, 1.7) than Salpeter IMF.

Next we developed a program for a cluster with IMFs 1.3, 2.3, 3.3 and for a cluster with IMFs 1.7, 2.3, 3.3.

All programs are developed in the R language and are given in the appendix.

3.4 Results

Calculation of unknown mass

Coefficients:

(Intercept) w

2.945 -1.349

W is the slope of the graph which is the same as Salpeter’s slope.

Plot 3.1: Salpeter’s slope

Binary frequency as a function of primary mass for a flat IMF

Case 1: slope =1.3

1. Sd = 0.1

“bounded”

499500

“notbounded”

0

Plot 3.2: obtained by taking standard deviation = 0.1

2. Sd = 0.2

“bounded”

499198

“notbounded”

302

Plot 3.3: obtained by taking standard deviation = 0.2

Sd = 0.3

“bounded”

493810

“notbounded”

5690

Plot 3.4: Obtained by taking standard deviation = 0.3

Sd = 0.4

“bounded”

482970

“notbounded”

16530

Plot 3.5: Obtained by taking standard deviation = 0.4

Sd = 0.5

“bounded”

467039

“notbounded”

32461

Plot 3.6: Obtained by taking standard deviation = 0.5

Sd = 0.6

“bounded”

451460

“notbounded”

48040

Plot 3.7: Obtained by taking standard deviation = 0.6

Sd = 0.7

“bounded”

430286

“notbounded”

69214

Plot 3.8: Obtained by taking standard deviation = 0.7

Sd = 0.8

“bounded”

413555

“notbounded”

85945

Plot 3.9: Obtained by taking standard deviation = 0.8

Sd = 0.9

“bounded”

395251

“notbounded”

104249

Plot 3.10: Obtained by taking standard deviation = 0.9

Sd = 1

“bounded”

380665

“notbounded”

118835

Plot 3.11: Obtained by taking standard deviation = 1

Sd = 1.1

“bounded”

367988

“notbounded”

131512

Plot 3.12: Obtained by taking standard deviation = 1.1

Sd = 1.2

“bounded”

350388

“notbounded”

149112

Plot 3.13: Obtained by taking standard deviation = 1.2

Sd = 1.3

“bounded”

345426

“notbounded”

154074

Plot 3.14: Obtained by taking standard deviation = 1.3

Case 2: Slope = 1.7

Sd = 0.1

“bounded”

499500

“notbounded”

0

Plot 3.15: Obtained by taking standard deviation = 0.1

2. Sd = 0.2

“bounded”

499198

“notbounded”

302

Plot 3.16: Obtained by taking standard deviation = 0.2

3. Sd = 0.3

“bounded”

493810

“notbounded”

5690

Plot 3.17: Obtained by taking standard deviation = 0.3

Sd = 0.4

“bounded”

482970

“notbounded”

16530

Plot 3.18: Obtained by taking standard deviation = 0.4

Sd = 0.5

“bounded”

467039

“notbounded”

32461

Plot 3.19: Obtained by taking standard deviation = 0.5

Sd = 0.6

“bounded”

451460

“notbounded”

48040

Plot 3.20: Obtained by taking standard deviation = 0.6

Sd = 0.7

“bounded”

430286

“notbounded”

69214

Plot 3.21: Obtained by taking standard deviation = 0.7

Sd = 0.8

“bounded”

413555

“notbounded”

85945

Plot 3.22: Obtained by taking standard deviation = 0.8

Sd = 0.9

“bounded”

395251

“notbounded”

104249

Plot 3.23: Obtained by taking standard deviation = 0.9

Sd = 1

“bounded”

380665

“notbounded”

118835

Plot 3.24: Obtained by taking standard deviation = 1

Sd = 1.1

“bounded”

367988

“notbounded”

131512

Plot 3.25: Obtained by taking standard deviation = 1.1

Sd = 1.2

“bounded”

350388

“notbounded”

149112

Plot 3.26: Obtained by taking standard deviation = 1.2

Sd = 1.3

“bounded”

345426

“notbounded”

154074

Plot 3.27: Obtained by taking standard deviation = 1.3

Binary frequency as a function of primary mass for steep IMF

Slope = 3.3

Sd = 0.1

“bounded”

499499

“notbounded”

1

Plot 3.28: Obtained by taking standard deviation = 0.1

Sd = 0.2

“bounded”

498969

“notbounded”

531

Plot 3.29: Obtained by taking standard deviation = 0.2

Sd = 0.3

“bounded”

489340

“notbounded”

10160

Plot 3.30: Obtained by taking standard deviation = 0.3

Sd = 0.4

“bounded”

472884

“notbounded”

26616

Plot 3.31: Obtained by taking standard deviation = 0.4

Sd = 0.5

“bounded”

444376

“notbounded”

55124

Plot 3.32: Obtained by taking standard deviation = 0.5

Sd = 0.6

“bounded”

423807

“notbounded”

75693

Plot 3.33: Obtained by taking standard deviation = 0.6

Sd = 0.7

“bounded”

395298

“notbounded”

104202

Plot 3.34: Obtained by taking standard deviation = 0.7

Sd = 0.8

“bounded”

375378

“notbounded”

124122

Plot 3.35: Obtained by taking standard deviation = 0.8

Sd = 0.9

“bounded”

355984

“notbounded”

143516

Plot 3.36: Obtained by taking standard deviation = 0.9

Sd = 1

“bounded”

333765

“notbounded”

165735

Plot 3.37: Obtained by taking standard deviation = 1

Sd = 1.1

“bounded”

316659

“notbounded”

182841

Plot 3.38: Obtained by taking standard deviation = 1.1

Sd = 1.2

“bounded”

301037

“notbounded”

198463

Plot 3.39: Obtained by taking standard deviation = 1.2

Sd = 1.3

“bounded”

287555

“notbounded”

211945

Plot 3.40: Obtained by taking standard deviation = 1.3

Binary frequency as a function of stellar mass for a cluster with IMFs

1.3, 2.3, 3.3.

Case1: In these graphs normalization is done by diving each and every y value by sum of all y values.

1.

Plot 3.41: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 0.1

2.

Plot 3.42: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 0.2

3.

Plot 3.43: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 0.3

4.

Plot 3.44: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 0.4

5.

Plot 3.45: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 0.5

6.

Plot 3.46: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 0.6

7.

Plot 3.47: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 0.7

8.

Plot 3.48: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 0.8

9.

Plot 3.49: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 0.9

10.

Plot 3.50: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 1

11.

Plot 3.51: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 1.1

12

Plot 3.52: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 1.2

13.

Plot 3.53: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 1.3

Case 2: In these graphs normalization is done by dividing each and every y value by ym, where ym is the value of y when x = log0

1.

Plot 3.54: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 0.1

2.

Plot 3.55: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 0.2

3.

Plot 3.56: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 0.3

4.

Plot 3.57: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 0.4

5.

Plot 3.58: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 0.5

6.

Plot 3.59: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 0.6

7.

Plot 3.60: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 0.7

8.

Plot 3.61: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 0.8

9.

Plot 3.62: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 0.9

10.

Plot 3.63: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 1

11.

Plot 3.64: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 1.1

12.

Plot 3.65: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 1.2

13.

Plot 3.66: A cluster with IMFs 1.3, 2.3, 3.3 having standard deviation = 1.3

Binary frequency as a function of stellar mass for a cluster with IMFs

1.7, 2.3, 3.3.

Case1: In these graphs normalization is done by diving each and every y value by sum of all y values.

1.

Plot 3.67: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 0.1

2.

Plot 3.68: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 0.2

3.

Plot 3.69: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 0.3

4.

Plot 3.70: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 0.4

5.

Plot 3.71: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 0.5

6.

Plot 3.72: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 0.6

7.

Plot 3.73: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 0.7

8.

Plot 3.74: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 0.8

9.

Plot 3.75: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 0.9

10.

Plot 3.76: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 1

11.

Plot 3.77: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 1.1

12.

Plot 3.78: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 1.2

13.

Plot 3.79: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 1.3

Case 2: In these graphs normalization is done by dividing each and every y value by ym, where ym is the value of y when x = log0 1.

1.

Plot 3.80: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 0.1

2.

Plot 3.81: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 0.2

3.

Plot 3.82: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 0.3

4.

Plot 3.83: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 0.4

5.

Plot 3.84: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 0.5

6.

Plot 3.85: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 0.6

7.

Plot 3.89: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 0.7

8.

Plot 3.90: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 0.8

9.

Plot 3.91: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 0.9

10.

Plot 3.92: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 1

11.

Plot 3.93: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 1.1

12.

Plot 3.94: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 1.2

13.

Plot 3.95: A cluster with IMFs 1.7, 2.3, 3.3 having standard deviation = 1.3

Special cases

Case 1: slope = 1.3

Plot 3.96: Obtained by dividing each bin number by 499500 for all standard deviation from 0.1 to 1.3.

Case 2: Slope = 1.7

Plot 3.97: Obtained by dividing each bin number by 499500, for all standard deviation from 0.1 to 1.3.

Case 3: Slope = 2.3

Plot 3.98: Obtained by dividing each bin number by 499500, for all standard deviation from 0.1 to 1.3.Case 4: Slope = 3.3

Plot 3.99: Obtained by dividing each bin number by 499500, for all standard deviation from 0.1 to 1.3.

APPENDIX

R Program to find the unknown mass

r=runif(1000,min=0,max=1)

m_up=0.4

m_low=10

beta;-2.35

bt;-(-beta+1)

n;-1/btm;- r*((m_up^bt)-(m_low^bt))+(m_low^bt)

mass;- m^(1/bt)

j;-log(mass)

h;-hist(x;-j,xlab = ‘log(mass)’,main = “Histogram”,

breaks=100)

y;-log(h$counts)

w;-h$midsis.na(y);-sapply(y,is.infinite)

print(y,na.rm=TRUE)

yis.na(y) ;- 0

plot(w,y,xlab=”log(mass)”,ylab=”Relative Frequency”)

fit;-lm(y~w)

fitBinary frequency as a function of primary mass for Salpeter IMF

p=1000

r=runif(p,min=0,max=1)

m_up=0.4

m_low=10

beta;-2.3

bt;-(-beta+1)

n;-1/btm;- r*((m_up^bt)-(m_low^bt))+(m_low^bt)

mass;- m^(1/bt)

md;-vector()

bnc;-vector()

bncn;-vector()

for(s in 1:13){

sd=s*0.1

v=rnorm(p,0,sd)

bounded=0

notbounded=0

pm;-vector()

for(i in 1:p){

m1=getElement(mass,i)

r1=getElement(r,i)

v1=getElement(v,i)

for(j in 2:p){

if(i;j){

m2=getElement(mass,j)

r2=getElement(r,j)

v2=getElement(v,j)

R=abs(r1-r2)

G=1

PE=(G*m1*m2)/R

V=abs(v1-v2)

if(m1;m2){

M=m1

}

else{ M=m2

}

KE=M*(V^2)

if(PE;KE){

if(M==m1){M1=M

M2=m2

}

else{M1=M

M2=m1}

df;-data.frame(M1,M2)

cn;-log(M1)

pm;-c(pm,cn)

bounded=bounded+1

}

else{ notbounded=notbounded+1

}

}

else{j=j+1

}

}

}

print(“bounded”)

print(bounded)

print(“notbounded”)

print(notbounded)

h;-hist(x;-pm,xlab = “log(primary mass)”,main = “Histogram of Primary Mass”)

y;-h$counts yn;-y/sum(y)

w;-h$mids plot(w,yn,xlab=”bin midvalue”,ylab=”bin number”)

bnc;-c(bnc,”bin numbers are respectively”,y)

bnc1;-data.frame(bnc)

bncn;-c(bncn,”bin numbers(normalized)”,yn)

bncn1;-data.frame(bncn)

md;-c(md,”binmid values are respectively”,w)

mdl;-data.frame(md)

}

write(md,file=”C://Users//pnvko//Downloads//mid2.3″)

write(bnc,file=”C://Users//pnvko//Downloads//bin2.3″)

write(bncn,file=”C://Users//pnvko//Downloads//nrbin2.3″)

Binary frequency as a function of primary mass for a flat IMF

Case 1: Slope = 1.3

p=1000

r=runif(p,min=0,max=1)

m_up=0.4

m_low=10

beta;-1.3

bt;-(-beta+1)

n;-1/btm;- r*((m_up^bt)-(m_low^bt))+(m_low^bt)

mass;- m^(1/bt)

md;-vector()

bnc;-vector()

bncn;-vector()

for(s in 1:13){

sd=s*0.1

v=rnorm(p,0,sd)

bounded=0

notbounded=0

pm;-vector()

for(i in 1:p){

m1=getElement(mass,i)

r1=getElement(r,i)

v1=getElement(v,i)

for(j in 2:p){

if(i;j){

m2=getElement(mass,j)

r2=getElement(r,j)

v2=getElement(v,j)

R=abs(r1-r2)

G=1

PE=(G*m1*m2)/R

V=abs(v1-v2)

if(m1;m2){

M=m1

}

else{ M=m2

}

KE=M*(V^2)

if(PE;KE){

if(M==m1){M1=M

M2=m2

}

else{M1=M

M2=m1}

df;-data.frame(M1,M2)

cn;-log(M1)

pm;-c(pm,cn)

bounded=bounded+1

}

else{ notbounded=notbounded+1

}

}

else{j=j+1

}

}

}

print(“bounded”)

print(bounded)

print(“notbounded”)

print(notbounded)

h;-hist(x;-pm,xlab = “log(primary mass)”,main = “Histogram of Primary Mass

“)

y;-h$counts yn;-y/sum(y)

w;-h$mids plot(w,yn,xlab=”bin midvalue”,ylab=”bin number”)

bnc;-c(bnc,”bin numbers are respectively”,y)

bnc1;-data.frame(bnc)

bncn;-c(bncn,”bin numbers(normalized)”,yn)

bncn1;-data.frame(bncn)

md;-c(md,”binmid values are respectively”,w)

mdl;-data.frame(md)

}

write(md,file=”C://Users//pnvko//Downloads//mid1.3″)

write(bnc,file=”C://Users//pnvko//Downloads//bin1.3″)

write(bncn,file=”C://Users//pnvko//Downloads//nrbin1.3″)

Case 2: Slope 1.7

p=1000

r=runif(p,min=0,max=1)

m_up=0.4

m_low=10

beta;-1.7

bt;-(-beta+1)

n;-1/btm;- r*((m_up^bt)-(m_low^bt))+(m_low^bt)

mass;- m^(1/bt)

md;-vector()

bnc;-vector()

bncn;-vector()

for(s in 1:13){

sd=s*0.1

v=rnorm(p,0,sd)

bounded=0

notbounded=0

pm;-vector()

for(i in 1:p){

m1=getElement(mass,i)

r1=getElement(r,i)

v1=getElement(v,i)

for(j in 2:p){

if(i;j){

m2=getElement(mass,j)

r2=getElement(r,j)

v2=getElement(v,j)

R=abs(r1-r2)

G=1

PE=(G*m1*m2)/R

V=abs(v1-v2)

if(m1;m2){

M=m1

}

else{ M=m2

}

KE=M*(V^2)

if(PE;KE){

if(M==m1){M1=M

M2=m2

}

else{M1=M

M2=m1}

df;-data.frame(M1,M2)

cn;-log(M1)

pm;-c(pm,cn)

bounded=bounded+1

}

else{ notbounded=notbounded+1

}

}

else{j=j+1

}

}

}

print(“bounded”)

print(bounded)

print(“notbounded”)

print(notbounded)

h;-hist(x;-pm,xlab = “log(primary mass)”,main = “Histogram of Primary Mass

“)

y;-h$counts yn;-y/sum(y)

w;-h$mids plot(w,yn,xlab=”bin midvalue”,ylab=”bin number”)

bnc;-c(bnc,”bin numbers are respectively”,y)

bnc1;-data.frame(bnc)

bncn;-c(bncn,”bin numbers(normalized)”,yn)

bncn1;-data.frame(bncn)

md;-c(md,”binmid values are respectively”,w)

mdl;-data.frame(md)

}

write(md,file=”C://Users//pnvko//Downloads//mid1.7″)

write(bnc,file=”C://Users//pnvko//Downloads//bin1.7″)

write(bncn,file=”C://Users//pnvko//Downloads//nrbin1.7″)

Binary frequency as a function of primary mass for steep IMF

Slope = 3.3

p=1000

r=runif(p,min=0,max=1)

m_up=0.4

m_low=10

beta;-3.3

bt;-(-beta+1)

n;-1/btm;- r*((m_up^bt)-(m_low^bt))+(m_low^bt)

mass;- m^(1/bt)

md;-vector()

bnc;-vector()

bncn;-vector()

for(s in 1:13){

sd=s*0.1

v=rnorm(p,0,sd)

bounded=0

notbounded=0

pm;-vector()

for(i in 1:p){

m1=getElement(mass,i)

r1=getElement(r,i)

v1=getElement(v,i)

for(j in 2:p){

if(i;j){

m2=getElement(mass,j)

r2=getElement(r,j)

v2=getElement(v,j)

R=abs(r1-r2)

G=1

PE=(G*m1*m2)/R

V=abs(v1-v2)

if(m1;m2){

M=m1

}

else{ M=m2

}

KE=M*(V^2)

if(PE;KE){

if(M==m1){M1=M

M2=m2

}

else{M1=M

M2=m1}

df;-data.frame(M1,M2)

cn;-log(M1)

pm;-c(pm,cn)

bounded=bounded+1

}

else{#print(“not bounded”)

notbounded=notbounded+1

}

}

else{j=j+1

}

}

}

print(“bounded”)

print(bounded)

print(“notbounded”)

print(notbounded)

h;-hist(x;-pm,xlab = “log(primary mass)”,main = “Histogram of Primary Mass

“)

y;-h$counts yn;-y/sum(y)

w;-h$mids plot(w,yn,xlab=”bin midvalue”,ylab=”bin number”)

bnc;-c(bnc,”bin numbers are respectively”,y)

bnc1;-data.frame(bnc)

bncn;-c(bncn,”bin numbers(normalized)”,yn)

bncn1;-data.frame(bncn)

md;-c(md,”binmid values are respectively”,w)

mdl;-data.frame(md)

}

write(md,file=”C://Users//pnvko//Downloads//mid3.3″)

write(bnc,file=”C://Users//pnvko//Downloads//bin3.3″)

write(bncn,file=”C://Users//pnvko//Downloads//nrbin3.3″)

4. Binary frequency as a function of stellar mass for a cluster with IMFs

1.3, 2.3, 3.3.Case 1

xdata;-scan(file=”C://Users//pnvko//Downloads//midstd”)

y1;-scan(file=”C://Users//pnvko//Downloads//std1.3//std0.1-1.3″) # changing path by modifying the file name from std0.1-1.3 up to std1.3-1.3 different graphs can be obtained.

y2;-scan(file=”C://Users//pnvko//Downloads//std2.3//std0.1-2.3″)

y3;-scan(file=”C://Users//pnvko//Downloads//std3.3//std0.1-3.3″)

plot(xdata,y1,type=”o”,col=”red”,main=”sd=0.1″,xlim=c(-0.9,1.7),ylim=c(0.00,0.14),xlab=”Bin midvalue”,ylab=”Bin number”)

points(xdata,y2,type=”o”,col=”cyan”,pch=”+”)

points(xdata,y3,type=”o”,col=”blue”,pch=”*”)

Case 2

xdata;-scan(file=”C://Users//pnvko//Downloads//midstd”)

y1;-scan(file=”C://Users//pnvko//Downloads//std1.3//log0.1″) # changing path by modifying the file name from log0.1 up to log1.3 different graphs can be obtained.

y2;-scan(file=”C://Users//pnvko//Downloads//std2.3//log0.1″)

y3;-scan(file=”C://Users//pnvko//Downloads//std3.3//log0.1″)

plot(xdata,y1,type=”o”,col=”red”,main=”sd=0.1″,xlim=c(-0.9,1.7),ylim=c(0.00,2.5),xlab=”Mass”,ylab=”Relative frequency”)

points(xdata,y2,type=”o”,col=”cyan”,pch=”+”)

points(xdata,y3,type=”o”,col=”blue”,pch=”*”)

5. Binary frequency as a function of stellar mass for a cluster with IMFs

1.7, 2.3, 3.3.

Case 1:

xdata;-scan(file=”C://Users//pnvko//Downloads//midstd”)

y1;-scan(file=”C://Users//pnvko//Downloads//std1.7//std1.3-1.7″) # changing path by modifying the file name from std0.1-1.3 up to std1.3-1.3 different graphs can be obtained.

y2;-scan(file=”C://Users//pnvko//Downloads//std2.3//std1.3-2.3″)

y3;-scan(file=”C://Users//pnvko//Downloads//std3.3//std1.3-3.3″)

plot(xdata,y1,type=”o”,col=”red”,main=”sd=1.3″,xlim=c(-0.9,1.7),ylim=c(0.00,0.14),xlab=”Bin midvalue”,ylab=”Bin number”)

points(xdata,y2,type=”o”,col=”cyan”,pch=”+”)

points(xdata,y3,type=”o”,col=”blue”,pch=”*”)

Case 2:

xdata;-scan(file=”C://Users//pnvko//Downloads//midstd”)

y1;-scan(file=”C://Users//pnvko//Downloads//std1.7//log0.1″)

y2;-scan(file=”C://Users//pnvko//Downloads//std2.3//log0.1″)

y3;-scan(file=”C://Users//pnvko//Downloads//std3.3//log0.1″)

plot(xdata,y1,type=”o”,col=”red”,main=”sd=0.1″,xlim=c(-0.9,1.7),ylim=c(0.00,2.5),xlab=”Mass”,ylab=”Relative frequency”)

points(xdata,y2,type=”o”,col=”cyan”,pch=”+”)

points(xdata,y3,type=”o”,col=”blue”,pch=”*”)

Special Cases

Case 2: Slope = 1.7

xdata;-scan(file=”C://Users//pnvko//Downloads//midstd”)

y1;-scan(file=”C://Users//pnvko//Downloads//std1.7//nbin0.1″)

y2;-scan(file=”C://Users//pnvko//Downloads//std1.7//nbin0.2″)

y3;-scan(file=”C://Users//pnvko//Downloads//std1.7//nbin0.3″)

y4;-scan(file=”C://Users//pnvko//Downloads//std1.7//nbin0.4″)

y5;-scan(file=”C://Users//pnvko//Downloads//std1.7//nbin0.5″)

y6;-scan(file=”C://Users//pnvko//Downloads//std1.7//nbin0.6″)

y7;-scan(file=”C://Users//pnvko//Downloads//std1.7//nbin0.7″)

y8;-scan(file=”C://Users//pnvko//Downloads//std1.7//nbin0.8″)

y9;-scan(file=”C://Users//pnvko//Downloads//std1.7//nbin0.9″)

y10;-scan(file=”C://Users//pnvko//Downloads//std1.7//nbin1″)

y11;-scan(file=”C://Users//pnvko//Downloads//std1.7//nbin1.1″)

y12;-scan(file=”C://Users//pnvko//Downloads//std1.7//nbin1.2″)

y13;-scan(file=”C://Users//pnvko//Downloads//std1.7//nbin1.3″)

plot(xdata,y1,type=”o”,main=”Beta=1.7″,xlim=c(-0.9,2.3),xlab=”Mass”,ylab=”Relative Frequency”)

points(xdata,y2,type=”o”)

points(xdata,y3,type=”o”)

points(xdata,y4,type=”o”)

points(xdata,y5,type=”o”)

points(xdata,y6,type=”o”)

points(xdata,y7,type=”o”)

points(xdata,y8,type=”o”)

points(xdata,y9,type=”o”)

points(xdata,y10,type=”o”)

points(xdata,y11,type=”o”)

points(xdata,y12,type=”o”)

points(xdata,y13,type=”o”)

Note: By changing the path properly case 1, case 3, case 4 can be obtained from the same program

References

1 Alan H. Batten, Binary and Multiple Systems of Stars, Pergamon Press Ltd, 1973

2 Robert Grant Aitken, The Binary Stars, Douglas C. McMurtrie, 1918

3 G. Johnstone Stone, On the Physical Constitution of the Sun and Stars, The Royal Society, 16 (1867), 25-34