Designing New Chemicals
Applicability: Synapse (core versions 0315+)
The goal of a chemical design is to generate molecular
structures whose physical properties satisfy a set of
constraints. Thus, to design chemicals, Synapse requires:
-
Groups: the pieces of molecular structure from
which new chemicals will be designed.
-
Limits: the minimum and maximum number of times each
group can occur in a designed structure as well as the
minimum and maximum number of groups and rings that
designed structures may contain.
-
Constraints: once molecular structures have been
created, we can estimate the candidate chemical's
physical properties. These properties are then used to
evaluate constraints.
Refractive Index Matching Fluids
Refractive index matching fluids are often used to fill the
gaps between the ends of two optical fibers. If the fluid
matches the refractive index of the fibers, there will be
very little reflective loss as the light travels from one
fiber to the next. In addition to constraints on refractive
index, it is desired that candidate matching fluids have
low freezing points, low vapor pressures and low flash
points. It is also important that its refractive index
does not change greatly as the temperature changes.
The following examples detail the steps needed to design
refractive index matching fluids that have acceptable
physical properties.
Example: Evaluating property estimation methods
Synapse combines design groups into the molecular structures of
candidate chemicals. Synapse then estimates those physical properties
of each candidate chemical needed to evaluate the design's constraints.
Thus, every chemical design requires an associated knowledge base
which is used as the source of design groups and estimation techniques.
Thus, the first step in a chemical design is develop a
knowledge base that contains the all the desired design groups
and whose estimation techniques can accuractely
predict the properties needed to evaluate design constraints.
Often the MKS Core Knowledge Base will be sufficient for
a design. However, it is possible that the knowledge base
will need to be expanded by adding new groups, new physical
property data and new estimation techniques.
At a minimum, you should evaluate your design's constraints
on existing chemicals to determine the accuracy of the
knowledge base's estimation techniques.
-
Open the MKS Core Knowledge Base document. (Create a
copy of the document
(see here)
to use for these examples.)
-
Change to the Techniques Chapter and navigate to the
RI,l,293: Vogel nD Method.
(See the
Navigation Overview documentation
for details on navigating chapters and pages.)
-
Select the Evaluate Technique command from the chapter's
Commands menu. The application will use the current
technique to estimate the refractive index of every
valued chemical in the current knowledge base.
Once the application has evaluated all chemicals having
data values, it displays the results in the Technique
Evaluation dialog.
1
|
The Evaluation Records table displays data values,
estimated values, errors, absolute errors, percent
errors and absolute percent errors for each
chemical.
|
2
|
The Statistics table displays summary statistics
for all "considered" records. Records that have
been ignored are not included in the calculation
of summary statistics.
|
3
|
The Graphs control displays several scatter plots
and histograms comparing data values and estimates.
Clicking and dragging the left mouse button to
enclose data points will select their corresponding
records in teh Evaluation Records table.
|
The Technique Evaluation dialog shows that Vogel's technique
generates estimates with an average error of ±0.006 and
a maximum absolute error of ±0.0371 for the 124
chemicals analyzed. These errors are quite low. Therefore,
Vogel's technique should give good design results without
modification.
Example: Creating a design function
Once we have verified that we can adequately estimate
the physical properties we plan to use in our design,
we can develop the design functions that use these
estimated physical properties.
-
Select the New command from the File menu. The
application will activate the Create a New
Document dialog. Select the Chemical Design
Document document type and press the dialog's
OK button.
-
The application will activate the File Dialog
prompting you for the filename of the
new document. Enter a name and press the dialog's
Save button. The application will create and open
the new design document.
-
Enter values for the document's title, subtitle
and descriptions. See documentation on the
Document Titles Section and
Document Information Section for
details.
-
Using the tabs at the top of the document, change to the
Functions Chapter by clicking the left mouse button on
Functions tab. (See the
Navigation Overview documentation
for details on navigating chapters and pages.)
-
Create a new design function entity by pressing the "+" button
in the menubar or executing the "Add New Page" command found
on the Edit menu. A new, blank Function page will be added
to the current document.
-
Click the left mouse button on the Identifier Pane's
edit control. The application will activate the edit
dialog. Enter the name "Refractive Index at 20°C" and
press the dialog's Save button.
-
Now click the left mouse button on the Function
Section's large edit control. The application will
activate the Function Code dialog. Enter the
following code:
// Default declarations
string prop, candidate;
double ri, temp;
int err;
// Default assignments
prop = "Refractive Index, Liquid at 293K";
candidate = Chemical();
temp = 273.15 + 20.0;
// Retrieve the candidate's refractive index
ri = CProp(candidate, prop, temp, 0, err);
if( err != 0 ) return FALSE;
// Assign estimate
SetResult(ri);
// Successful
return TRUE;
-
Finally, press the Code dialog's Save button.
-
Refractive index matching liquids should also
have a low vapor pressure so they to not
evaporate quickly from connections where they
are applied. We thus need to add another
design function to estimate the vapor pressure
at 20°C.
Repeat the previous steps to create another
design function. Name this function "Vapor
Pressure at 20°C [kPa]" and enter the
following code.
// Default declarations
string prop, candidate;
double pvp, temp;
int err;
// Initialize values
prop = "Vapor Pressure, Liquid - f(T)";
candidate = Chemical();
temp = 273.15 + 20.0;
// Retrieve vapor pressure value
pvp = CProp(candidate, prop, temp, 0, err);
if( err != 0 ) return FALSE;
// Convert from Pa to kPa
pvp = pvp / 1000.0;
// Assign value
SetResult(pvp);
// Successful
return TRUE;
Example: Testing a design function
Before performing a chemical design, it is important to ensure the
constraints and design functions have been entered correctly.
Synapse provides a testing mechanism in which you can run a
design function on an existing chemical entity.
-
Activate the copy of the knowledge base document we created
in the first example.
-
Change to the Chemicals Chapter and navigate to
Diethylene glycol dimethyl ether.
(See the
Navigation Overview documentation
for details on navigating chapters and pages.)
-
Select the Compute Estimates command from the Commands
menu. The application will activate the Property
Estimation Dialog.
-
Press the dialog's Start button. The application
will begin estimating all properties of the current
chemical.
-
Once all estimations have been performed, press the
dialog's Save button to store the estimated values
into the current document.
-
Scroll to the Optical Properties Section and make a
note of the estimated liquid refractive index, i.e.,
nD,l,293.
1
|
The estimation liquid refractive index
is 1.39251.
|
-
Now navigate to the "Refractive Index at 20°C" design
function in the Functions Chapter of our newly created
chemical design document.
-
Select the Test Chemical Function command from the Commands
menu.
The application will activate the Test Chemical Function
Dialog.
-
Select Diethylene glycol dimethyl ether as the chemical
candidate press the
dialog's Calculate button. Synapse will execute the design
function using the selected chemical candidate and display
the results in the dialog. In this example, the design
function successfully calculated a refractive index of
1.39251, the same value estimated in the knowledge base.
Example: Create a combinatorial chemical design
-
Activate the chemical design document we created
in a previous example.
-
Using the tabs at the top of the document, change to the
Combinatorial Designs Chapter by clicking the left mouse
button on the Combinatorials tab.
-
Create a new Combinatorial design by pressing the "+" button in the
menubar or executing the "Add New Page" command found on the
Edit menu. A new, blank page will be added to the current
document.
-
Click the left mouse button in the Identifier Pane, the large
white box at the top of the page. The application will activate
the pane's datum edit dialog.
-
Enter a name for the new design. Optionally enter a
reference and comment.
-
Finally, press the dialog's Save button. The application
will save the new design's identifier into the
current document and display the new name you just
entered.
Example: Specify the design's knowledge base document
A combinatorial chemical design will assemble groups
in all possible combinations to generate candidate molecular
structures, use estimation techniques to obtain physical
properties for each candidate
chemical and then use these physical property values to evaluate
each design constraint. The groups and estimation
techniques used in this process are retrieved from the design's
associated knowledge base.
-
Ensure the copy of the knowledge base document we created
in the first example is open.
-
Navigate to the new combinatorial chemical design we created
in the previous example.
-
Click the left mouse button on the Source Knowledge Base
Section's edit control. The application will activate the
Knowledge Base selection dialog.
-
Select the knowledge base we have been using in these
examples and press the dialog's OK button. The application
will store this selection into the current document.
Example: Specifying design parameters
Synapse generates candidate chemicals by assembling a given
set of groups in all possible combinations. The limits on
the total number of groups to use and the total number of
rings that may occur in these generated molecular
structures are
specified in the Design Parameters Section.
For this example we will use the following design
parameters:
Parameter |
Minimum |
Maximum |
# Groups |
6 |
10 |
# Rings |
0 |
0 |
-
Activate the chemical design document we created
in a previous example.
-
Using the tabs at the top of the document, change to the
Combinatorials Chapter by clicking the left mouse button on
Combinatorials tab.
-
In the Design Parameters Section click the left mouse
button on the Min # Groups field's edit control. The
application activates the Enter Minimum Groups
dialog.
-
Enter a value of 6 and, optionally, a comment. Then
press the dialog's OK button.
-
Enter the remaining design parameters.
1
|
Candidate chemicals are limited to having
between 6 and 10 groups.
|
2
|
Candidate chemicals are limited to having
between 0 and 0 ring, i.e., candidate
chemicals cannot contain rings.
|
Example: Specifying design groups
Synapse generates candidate chemicals by assembling
design groups into new molecular structures. The Design
Groups Section is used to specify these design groups as
well as the limits on their occrrence in new molecular
structures.
For example, the previous image shows that five groups
can be used to generate new molecular structures. It
also shows that every new molecular structure must have
between 0 and 2 -OH groups.
-
Click the left mouse button on the Design Groups Section's
table control. The application activates the Design Groups
Limits Dialog.
-
Click the left mouse button on the first row of the
Design Groups Limits Dialog's table control and then
press the dialog's Edit button. The application
activates the Design Group Limits Dialog.
-
Enter the name of a design group, the minimum number of
times the group must occur in a candidate structure, the
maximum number of times the group may occur in a
candidate structure and, optionally, values for the
reference and comments. Finally, press the dialog's
OK button.
-
Continue adding all the design groups and limits
shown in the image below. Once all the groups and
limits have been entered, press the dialog's Save
button. The values are stored in the current
document.
Example: Specifying substructure limits
Synapse bonds design groups togehter in all possible
combinations to form candidate molecular structures.
Some combinations of design groups may not be
desired due to chemical stability or reactivity. For
example, the design groups presented in the previous
example can be bonded for form acetals and
peroxides, e.g.,
For the design of refractive index matching fluids,
such unstable chemicals are not desired.
The Sustructure Limits Section enables you to enter
limits on the occurrence of non-design groups. Typically,
the groups you enter into this section are larger
groups which require or disallow particular combinations
of design groups.
-
Click the left mouse button on the Substructure Limits
Section's table control. The application activates the
Substructures Limits dialog.
-
Click the left mouse button on the table control's first
row and then press the dialog's Edit button. The application
activates the Substructure Limits dialog.
-
Enter "-O-CH2-O-", the name of the group to restrict, 0,
the minimum number of times the group may occur in a
candidate structure, 0, the maximum number of times the
group may occur in a
candidate structure and, optionally, values for the
reference and comments. Finally, press the dialog's
OK button.
-
Repeat the previous steps to add a restriction on the
occurrence of the "-O-O-" group. Then press the dialog's
Save button. The values are stored in the current
document.
Example: Entering design constraints
Each chemical design contains one or more constraints which viable
candidates must satisfy. Design constraints can be imposed on a
single physical property, such as density or viscosity, or on a
complex function of physical properties such as a heat transfer
coefficient calculated by a Nusselt number correlation.
Each constraint contains a function name, a minimum value, a
goal value and a maximum value. Some example constraints are
shown in the following table.
Function |
Minimum |
Goal |
Maximum |
Freezing Point [°C] |
-40 |
-10 |
0 |
Liquid Density at 20°C [kg/m3] |
850 |
1000 |
1150 |
Heat Transfer Coefficient at 20°C [W/m2 K] |
1000 |
1200 |
2000 |
We will use two constraints for our example refractive index
matchnig fluid design:
-
Refractive Index: the refractive index of viable
candidates must be close to 1.458, the refractive index of
the fused silica used in fiber optic cable.
-
Vapor Pressure: the vapor pressure of viable
candidates must be less than 0.40 kPa to reduce
evaporation.
Additional constraints on density and viscosity could
also be imposed but these two contraints will be
enough for demonstration purposes.
-
Navigate to the new combinatorial chemical design we created
in a previous example.
-
Scroll the datapane and click the left mouse button on the
Constraints Section's large table control. The application
activates the Edit Constraints dialog.
-
Click the left mouse button on the dialog table's first row.
Then press the dialog's Edit button. The application activates the
Edit Constraint dialog.
-
Enter "Refractive Index at 20°C", the name of the function we
created in a previous example into the Constraint Function
edit control. You can also press the dialog's List button for
a list of all functions present
in the current document. Enter a minimum value of 1.438, a goal
value of 1.458 and a maximum value of 1.478. Optionally, enter
values into the Reference and Comment controls. Finally, press
the dialog's OK button.
-
Repeat the previous steps to enter a second constraint on
the vapor pressure as shown in the following image.
-
Finally, press the Constraints Dialog to save the entered
constraints into the current document.
Example: Designing chemicals
-
Navigate to the new combinatorial mixture design we created
in a previous example.
-
Select the Design Candidates command from the Commands menu.
Synapse activates the Combinatorial Chemical Design dialog.
-
Press the dialog's Start button. Synapse will:
-
Generate candidate chemicals by combining all groups
in all possible ways.
-
Estimate the refractive index and vapor pressure
of each candidate chemical.
-
Check if these estimated properties satisfy
the entered constraints.
The Design dialog shows that Synapse generated 26,248
candidates, 329 of which satified our design
constraints.
-
Finally, press the dialog's Save button to store the
design results into the current document.
Tip: Use results to update constraints and substructure limits
The Design Dialog will often very good insight into how the
design inputs might be modified to produce improved results.
For example, the Vapor Pressure histogram, shown in the
Sequential Constraint Statistics section, shows that
candidates were rejected because their vapor pressure was
too low. Since, low vapor pressure chemicals are desired
in this design, this result indictes that we should modify
our vapor pressure constraint by lowering the range of
acceptable values.
Further examining the design's results, we will also find
candiates that contain the substructure:
Note that this substructure violates our desire to avoid
acetals (and hemiacetals, ketals, and hemiketals). However,
our current substructure limits do not restrict the
occurrence of this group. We should thus replace the
limit on the -O-CH2-O- group with a limit on a more
general group.
Example: Copy a template knowledge base
In the previous example, Synapse designed 329 candidate
chemicals that had refractive indices near 1.458 and
vapor pressures below 0.40 kPa. This
example shows how it is often useful to transfer
candidates to a knowledge base for further analysis.
-
Open the MKS TGER Template Knowledge Base document. The template
knowledge base contains techniques, elements, groups and
references but no chemicals and no mixtures. We
will use the template knowledge base as the starting
point for creating our compilation of design
candidates. (Note: the MKS TGER Template Knowledge Base can
be downloaded from the MKS Documents section of our
website's
Library page.)
-
Choose the 'Save a Copy' command from the File menu. Synapse
displays the Save As dialog.
-
Choose a location and enter a name for the knowledge base
copy, e.g., “RI Matching Fluids Candidates”.
-
Finally press the dialog’s Save button. The application will
display the Save Copy Progress dialog
showing the progress of the save operation. The dialog will
indicate when the document has been copied and give you the
option of opening the copied document.
-
Press the dialog's Open button.
-
Update the document's titles and descriptions. (See
documentation on the
Document Titles Section and the
Document Information Section for
editing details.)
Example: Analyze design candidates
-
Activate the chemical Design document we created and navigate
to the combinatorial chemical design we ran
previously.
-
Select the Transfer Candidates command from the Commands
menu. Synapse activates the Transfer Chemical Candidates
dialog.
-
In the Destination Knowledge Base control, select the name
of the newly copied template knowledge base. (IMPORTANT:
make sure the original MKS TGER Template Knowledge Base is
not selected.) Then press the
dialog's Select All button. All candidates
will be selected for transfer. Finally, press the
dialog's Transfer button. Synapse will create a new
chemical entity in the knowledge base for each chemical
candidate.
-
Once the transfer is complete, press the Transfer
Candidates Dialog's Done button.
-
Activate the knowledge base, change to the Chemicals
chapter and navigate to any one of the newly transferred
candidates.
-
Select the Estimate Multiple Chemicals command from
the Commands menu.
-
Synapse activates the Estimated Multiple Chemicals dialog.
Select "Refractive Index, Liquid at 293K" from the
Estimated Property control and press the dialog's
All button.
-
Press the dialog's Start button. Synapse will estimate
the liquid refractive index at 293K for every chemical
in the current knowledge base.
-
Finally, press the dialog's Save button. The application
stores the estimated properties into the current document.
Each newly designed candidate contained in the new knowledge
base can now be examined further. Refractive index can be
estimated at other temperatures. Vapor press curves can be
generated. Other properties such as viscosity and density
can also be estimated.
The main concept is that using Synapse, you generated a
new knowledge base containing refractive index matching
fluid candidates, many of which could be novel chemical
products.
Example: Examine the range of estimated property values
-
Click the right mouse button on the liquid refractive
index field's estimate control and select the
Compile Values command from the displayed commands
menu.
-
Synapse activates the Compile Property Values dialog.
Select the "Estimates Only" option from the Status
control and press the dialog's Compile button.
1
|
Select the Estimates Only option from
the Status control.
|
-
Synapse will compile all the refractive index values
estimated in the previous example and display them
in the Compiled Property Values dialog.
It is very interesting to note that the Maximum value
shown in the dialog's Value Statistics section is
below our goal value of 1.458, i.e., the refractive
index of fused silica. This might indicate that we
should perform another design with larger molecules
or additional design groups.
Related Documentation